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THE 1600 BAND SYSTEM AMMONIA! 


ABSTRACT 


The progression ammonia bands which extends from 1689 1400 has 
been photographed absorption high resolution. Six bands have been 
analyzed and found the perpendicular type. The analysis shows that the 
molecule planar the excited state and that vibrational levels observed the 
progression are those the out-of-plane vibration. The excited electronic state 
the type. addition the normal Coriolis interaction the degenerate 
levels, second effect has been observed which behaves like the Coriolis inter- 
action recently described ‘giant doubling’ Garing, Nielsen, and 
Rao. clear evidence has been found for any distortion the degenerate state 
from symmetry. 


INTRODUCTION 


The ultraviolet absorption spectrum ammonia has been the subject 
numerous investigations, the most extensive which have been those 
Duncan (1935), Duncan and Harrison (1936), and Walsh and Warsop (1960). 
Duncan (1935) has reported that there are four band systems between 2200 
and 1150 The first these, which has its origin near 2168 diffuse and 
its rotational structure can not resolved (the corresponding bands 
show some structure (Benedict 1935)). The second system, which lies 
between 1637 and 1400 much sharper and the rotational structure can 
resolved. The band systems shorter wavelengths have been discussed 
Walsh and Warsop (1960) but rather little known about their structures. 

The 1600 band system consists single progression bands with 
upper state frequency about 900 The first known member the 
progression, 1663 weak but the intensity increases rapidly along the 
progression and reaches its maximum the seventh eighth member. The 
rotational structures few the bands have been resolved Duncan and 
Harrison and they have given partial rotational analysis these bands. 
Unfortunately, their spectrograph lacked the resolving power necessary 
show much the detail the bands and will shown later that their 
analysis not correct. 
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Some the characteristics the upper state the 1600 band system 
can deduced from its vibrational structure. The fact that long progression 
bands observed indicates that there large change the dimensions 
the molecule during the transition. Though there proof, seems most 
reasonable assign the 900 frequency the upper state vibration 
which corresponds the (~950 vibration the ground state. Since 
the vibration corresponds normal mode which the height the 
ammonia pyramid changes, has been argued that the transition one 
which the height the pyramid undergoes large change and probably one 
which the upper state planar. 

The present work consists rotational analysis number the bands 
the 1600 system. Such analysis interest for number reasons. 
till now, band arising from electronic transition symmetric top 
has been analyzed detail. The analysis should give the geometry the 
molecule and the symmetry the electronic wave function the upper state. 
Also there exists the possibility that unusual types interaction between 
electronic and vibrational motion, such the effect, may 
discovered. For these reasons have undertaken rotational analysis the 
bands and have succeeded analyzing sufficient number show the 
essential characteristics the band system. 


EXPERIMENTAL 


The absorption spectrum ammonia from 1670 1400 was photo- 
graphed using 10-meter concave grating spectrograph. The spectrograph 
the normal incidence type and the seventh and eighth order spectra which 
were used this work had dispersions 0.21 A/mm and 0.18 A/mm respec- 
tively. Third and fourth order iron lines from hollow cathode source were 
used wavelength standards. Using Lyman discharge provide the 
continuum, spectra could obtained plates with exposure times 
which varied between minutes and hours. 

The absorption cells used this work varied length from cm. 
The ammonia gas pressure was varied from about micron 
cover the wide range intensities lines within the bands and the large 
difference the intensities the various bands. would have been desirable 
use even higher pressures obtain the weak bands near 1650 but was 
found that continuous absorption blocked all transmission the pressure 
was increased beyond about mm. Spectra the gas were also taken 
various temperatures ranging from —110°C +300°C. Finally, spectra 
the isotopically substituted molecule ND; were obtained but, these 
spectra have not been analyzed, they will not discussed further. 

All the lines the bands are broad and this line width and not the 
spectrograph that limits the resolution. The diffuseness the lines increases 
slowly towards shorter wavelengths. This diffuseness presumably caused 
predissociation and, with many other examples predissociation, 
much less noticeable the deuterium substituted compound. 
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III. THEORETICAL CONSIDERATIONS 


Before giving the details the analysis the most important results will 
given here and the evidence support them presented later. The most 
significant results are follows: 

(1) The bands are the perpendicular type, the upper state being 
degenerate electronic state species 

(2) The molecule planar the excited state. 

(3) The bands involve progression vibrational levels the out-of-plane 
vibration. 

(4) The upper state rotational levels show strong Coriolis-type perturbation 
similar that recently called giant doubling Garing, Nielsen, and 
Rao (1959). 

Now shall consider the characteristics band system arising from 
states having the properties listed above. its ground electronic state 
ammonia has pyramidal structure and C;, symmetry. The symmetry the 
ground state electronic wave function wish therefore consider 
transition between this lower state and upper state which the molecule 
planar (Ds, symmetry) and the electronic wave function has symmetry. 

Figure shows the vibrational energy levels the two states and the 
transitions allowed between them. the ground electronic state only the 
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Fic. vibrational energy levels ammonia which give rise the 1600 bands. The 
lower state levels are those caused the inversion NH; and the upper state levels are the 
out-of-plane vibrational levels planar ammonia molecule. The vibrational, electronic, and 
vibronic symmetries the levels are given the right-hand side. 
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two lowest levels, which are separated about are shown. This 
small separation results from the fact that the potential barrier preventing the 
nitrogen atom from passing through the plane low (Herzberg 1945, 
220). the two energy levels, one has wave function which symmetric 
with respect reflection plane that perpendicular the axis and passes 
through the center mass the molecule and the other has wave function 
antisymmetric with respect this plane. These levels are therefore usually 
referred the symmetric and the antisymmetric levels ammonia. the 
height the ammonia pyramid were reduced, thereby reducing the barrier 
the reflection, the symmetric and antisymmetric levels would move apart 
until the limit, where the ammonia molecule became planar, these two 
levels would become and vibrational levels the molecule, which now 
would possess D3, symmetry. For all barrier heights, however, the symmetric 
and antisymmetric levels have the characteristics the and levels and 
are therefore labelled such even the ground state ammonia 

Fig. the only vibrational levels shown the upper electronic state are 
those associated with the out-of-plane vibration. shall refer this vibra- 
tion the vibration. Since the transition involves large change the 
height the ammonia pyramid, from the Franck—Condon principle follows 
that should expect find long progression the vibration the 
absorption spectrum. Other vibrational levels the upper state are 
interest present and are not shown Fig. The vibration has 
symmetry and the higher vibrational levels alternate with and sym- 
metry shown the right-hand side Fig. important for subsequent 
discussion bear mind that the progression the only progression 
vibrational levels the planar ammonia molecule which alternates the 
sense that the wave functions are alternately symmetric and antisymmetric 
with respect the plane the molecule. The vibronic states (electronic 
shown Fig. 

The selection rules for the states shown Fig. are 
The allowed transitions are shown the figure. One the most significant 
facts that the symmetric level the lower state combines only with every 
other excited state vibrational level and the antisymmetric level combines 
with those between. Since the separation between the symmetric and 
antisymmetric levels the lower state only 0.66 this effect will not 
noticeable the positions the bands unless they are measured very 
precisely. This alternation does, however, have marked effect the rotational 
structures the bands. 

few the rotational levels state and the antisymmetric 
state the lower electronic state are shown Fig. The selection rules for 
the transition between these are +1, +1. Figure essentially 
the same that given Herzberg (1945, 408) where full description 
perpendicular-type transition given. The figure has been modified, how- 
ever, that the levels are those oblate rather than prolate top. The 
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Fic. few the rotational energy levels degenerate and non-degenerate state 
The lower set levels are those the antisymmetric ground state ammonia and 
the upper set are those vibronic state. The and rotational species have been 
separated for but not for other values where they are degenerate and are marked 
The allowed transitions (AK +1) are shown the arrows. 


rotational levels have been drawn dotted lines indicate that 
where the nuclear spin the three identical nuclei these levels not 
occur. This particularly important for the levels since here all odd 
levels the antisymmetric state are levels and therefore are missing. 
Furthermore, the state the and rotational levels with 
have been separated appreciable amount; effect which shall show 
later does indeed occur. Figure refers transition. The alter- 
nate bands given Fig. which are transitions, have very nearly the 
same structure, the one important difference being the fact that the and 
subband, the bands have lines with even values only while the 
bands have lines with odd values. 

order show that the observed band system does indeed result from the 
transition described above shall have show that has the following 
characteristics. 

(1) The bands are the perpendicular type. This will once establish the 
fact that the upper state vibronic levels are degenerate but somewhat more 
detailed argument will necessary show that the upper electronic state 
degenerate. 

(2) The bands the progression have for their lower states alternately the 
symmetric and the antisymmetric and vibrational levels the lower 
electronic state. This alternation will show that the upper state vibration 
the vibration. Proof that this alternation exists can provided the 
subbands which will have alternately only even only odd 
lines present. 
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(3) The level the upper state combines with the lower state 
level. This will establish the fact that the upper state rather than 

the following discussion the band system will shown that the 
three points listed above can established. addition the rotational con- 
stants the molecule its upper electronic state have been determined and 
number unusual characteristics the spectrum have been noted. 


IV. RESULTS 


Absorption spectra the 3-0 band taken three different pressures are 
shown Fig. was found that this large range pressures was necessary 
order bring out the essential features the band. Fig. shows, most 
the strong lines are near the center the band and this region almost all 
the are blends two more lines. 

The wave numbers the lines the bands are given Table The accuracy 
measurement varies good deal) from line line depending upon the 
sharpness the line and degree overlapping. probable that the 
average error about 0.07 but for some lines the error may great 
0.2 

The analysis the bands proved difficult task. Unlike the usual 
progressions bands, was found here that every band appeared quite 
different from all other bands. There did exist series about six strong lines 
each band which appeared similar all the bands but when these 
strong lines were correlated obvious correlations appeared exist between 
the other lines. 

Though the lines near the center the band show regularity, weak 
lines which occur the high-frequency end the bands form regular series 
and proved suitable point from which start the analysis. These 
regular series can seen Fig. Starting from these series, with the aid 
the excellent lower state combination differences Benedict, Plyler, and 
Tidwell (1958), was possible obtain analysis. The regular series proved 
branches subbands and the corresponding and branches could 
found. When all the lines the high-frequency series had been assigned 
subbands, together with their corresponding and lines there still remained 
many strong lines unassigned. this stage was apparent, both from the 
relative intensities the and branches the subbands and from the 
number subbands, that the band was the perpendicular type with two 
subbands for every value The subbands giving rise the regular 
structure described above were those with —1. extrapolating from 
these series was possible find the subbands and 
complete the analysis. 

The analysis shows that the major features the bands are just those 
expected for perpendicular bands symmetric top and are fact very 
similar the infrared bands NH; which have been analyzed Benedict, 
Plyler, and Tidwell. primarily the poor resolution, caused the large 
line width, together with the difficulty registering lines with wide differences 
intensities vacuum ultraviolet plates, which caused the difficulty the 
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analysis. will discussed later, the principal reason for every band appear- 
ing different from all others, lies the rapid decrease the upper state 
value with increasing The few strong lines which appear similar 
subbands. 

order specify band electronic transition ammonia such 
the one have photographed here necessary give all the vibrational 
quantum numbers both the upper and lower states together with some 
symbol denote the symmetry the lower state with respect reflection 
(inversion). the band system under investigation, however, only the 
vibration excited both the upper and lower states and shall designate 
bands the symbols where shall understood that all 
other vibrational quantum numbers are zero. Very often the symbol 
which differentiates between the lower state symmetric and antisymmetric 
levels, importance and will omitted. 

From the analysis the bands, two unexpected facts emerged once. 
Firstly, the long wavelength side the 1663 band which Duncan (1936), 
the basis the isotope effect, had assigned the 0-0 band, very weak 
band was found. discussed previously continuous absorption prevented 
this band from being investigated with long absorption paths. was first 
thought that this was the 0-1 band. The position the band with respect 
the band can calculated exactly from the infrared band which has 
been measured Garing, Nielsen, and Rao (1959). The positions the 
strongest lines the and the bands were calculated and was 
found that these did not correspond the positions the few lines which 
could measured the weak band. The weak band did, however, fall 
position such that could considered additional band the main 
progression. Thus was concluded that the weak band was the band and 
that the vibrational numbering given Duncan must increased one 
unit. 

definite proof that the weak band not the band was found when the 
spectrum ammonia 300° was examined. The spectrum showed number 
new lines the edge the weak band and their displacement from the 
1663 band was such that there was doubt that they arose from the 
1(s) level the lower state. Thus there appears doubt that Duncan’s 
vibrational numbering must increased least one unit. The question 
whether not can advanced more than one unit will dealt with 
later. the remainder this paper, the new numbering, which the 1663 
band the 1-0 band, will used. 

The second unexpected result arising from the analysis was the displacement 
branches these subbands had been found the expected positions the 
branch lines could determined. Since these branch lines should some 
the strongest lines the bands, they should have been found without 
difficulty but none the lines could found the expected positions. After 
the analysis all other subbands was completed, several strong lines were 
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TABLE 
Vacuum wave numbers and assignments for the lines the NH; bands 
band 
K’ = 7-K" =6 K’ = 6-K” =5 K’ = §5-K" =4 K’ = 4-K" =3 K’ = 1-K” =0 
1 59215.00 
2 
3 59194 .42 248.60 59180.97 
4 59190.11 213.06 59114.09 
5 59183.77 278.43 186.31 
6 59175.96 
7 305.04 
J” OO) P(J) P(J) 
4 | 59122.80 
6 59135 .46 
band 
K’ = 7-K"” =6 K’ = 4-K” =3 K’ = 1-K” =2 K’ = 2-K"” =3 K' = 3-K”" =4 
59205.18 
4 254.14 181.4 59190.24 
5 280.04 
6 59207 .42 
P(J) P(J) 
5 59178 .62 
59164.49 
band 
K’ = 7-K”" =6 K’ = 6-K” = 5 K’ = §-K” =4 K’ = 4-K” =3 
60126.08 
4 60120.75 145.73 60046 .46 59947 .34 
5 60113.40 140.33 60021 .39 165.00 046.46 
60104.18 159.62 021.39 183.82 
152.34 994.26 178.61 021.39 202.12 
142.91 965.24 171.13 994.26 197.16 219.87 
9 161.91 965.24 
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TABLE (Continued) 
Vacuum wave numbers and assignments for the lines the NH; bands 


60130.80 
187.18 59991.02 214.28 004.00 
221.87 970.16 
P(J) P(J) R(J) P(J) 
082.90 
059.79 60155.70 076.20 
156.17 
018.58 
000.56 
028.87 037.46 60048.51 
009.60 019.24 024.42 028.95 
001.85 005.90 012.00 
P(J) 
59992.91 
972.23 
band 
R(J) R(J) R(J) P(J) R(J) R(J) P(J) 
61055.92 
021.31 111.10 774.29 
064.47 746.48 


| 
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TABLE (Continued) 


Vacuum wave numbers and assignments for the lines the NH; bands 


0 61061 .08 
1 61061 .08 
2 | 61061.08 077 .52 61019.68 092.05 61038 .95 61000 .33 
3 075.23 61089 .50 093.11 017.86 
4 089 .50 60995 .40 107.18 013.80 117.63 034,10 60952 .85 
5 103.28 990.63  60896.65 119.52 008.80 
6 115.41 984.58 871.85 130.56 000 .33 137 .67 026.65 899 .93 
7 126.29 846.46 140.05 60990 .63 
8 135.85 819.74 148.10 982.11 153.10 017.86 841.99 
9 144.35 792.09 971.54 
10 763 .06 164.53 780.89 
K’ = 0-K” =1 K’ = 1-K” =2 K’ = 2-K" =3 
1 | 61095.35 61055.92  61038.95 
2 112.85 055.92 017.86 | 61128.39 61070.67 61034.10 
3 129.53 053.23 60995.40 144.27 —_ 011.79 | 61158.80 61083 .33 61026.65 
4 146.03 050.46 973.71 160.19 064.37 60989.25 174.24 079.43 003.81 
5 162.80 047.00 951.19 175.78 061.08 965.87 189.05 075.23 60980.25 
6 179.93 043.30 928.44 191.44 055.92 942.09 203 .99 070.67 956.42 
7 197 .56 906 .35 207 .51 053 .23 919.01 218.81 064.37 932.08 
8 884.46 896.62 234.47 061.08 908.21 
9 884.46 
K’' = 3-K” =4 K’ = 4-K” =5 K' = 5-K” =6 
4 | 61187.49 61093.11 61017.86 
5 202 .00 089.50 60994.21 | 61214.20 61101.11 61006 .70 
6 216.04 083 .33 969 .30 227.84 095.35  60982.11 | 61238.93 61107 .18 60994.21 
7 230.22 077 .52 944.87 241.08 089 . 50 956.42 251.45 100.22 968.25 
8 920.07 931.34 264 .09 092.05 942.09 
9 896.62 906.35 276.69 084.50 915.17 
10 880.82 887.61 
11 864.46 
K’ = 6-K” =7 K’ = 7-K” =8 K’ = 8-K" =9 K’ = 9-K” = 10 K’ = 10-K” = 11 
J” P(J) P(J) P(J) P(J) P(J) 
7 60980 . 25 
8 952.85 60963 
9 924.59 934.81 60944 .87 
10 899 .93 909.97 917.38 60924 .59 
11 60902 .62 
K’ = 11-K” = 12 K’ = 12-K” = 13 
P(J) P(J) 
12 60878 .98 
13 60852 .89 


| 
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Vacuum wave numbers and assignments for the lines the NH; bands 
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band 
K’ = =7 K’ = 7-K" = 6 K’ = 6-K" = 5 
5 62002.19 
6 61993 .42 008 .30 
7 61984.51 997 .37 015.99 
8 62000 . 58 61838 .86 61695 .98 015.99 61869.76 61711 .67 
9 001.44 822.57 661.45 013.01 854.63 677.50 
10 804.76 625.72 838.86 642.46 
ll 588.99 603 .99 
K' = 5-K”" =4 K’ = 4-K"” =3 K’ = 3-K” =2 
2 62018.73 
3 62015.14 030.74 61959. 23 
4 | 62009.59 025 .03 61935.75 040.30 951.22 
5 017.51 61910.32 033 .09 925.85 61836 .53 048.63 941.46 61852 .22 
6 023.85 898.65  61791.66 039 .73 914.38 806.77 054 .83 929 .92 822.57 
7 028.58 885.59 760.18 044.19 901.31 776.08 058.76 791.66 
8 870.19 727 .67 047.14 886 .69 743.42 060.90 758.95 
9 693 .89 048.63 870.19 709.70 724.76 
10 658.70 674.43 688 .82 
11 622.52 637 .97 651.64 
12 600.40 
= 2-K"” =1 K’ = 1-K” =0 = 0-K” =1 
1 | 62020.62 62035 .01 62001 .44 62053 . 16 62015 .99 62000 .58 
2 | 034 .63 61980 .86 067 .39 013.01 61977 .45 
3 046.18 974.85 057 .73 61994 .21 61935 .75 079.79 007 .61 953.90 
055.79 966.78  61895.47 091.40 000 .58 928.27 
5 062.99 956.52 867 .67 | 071.25 982.64 879.28 102.06 61993 .42 901.31 
6 068.14 944.37 838.11 111.89 982.64 873 .36 
7 071.25 929.93 806.41 076.88 955.74 814.52 845.06 
8 071.25 772.48 816.36 
9 737 .04 076.88 742.90 785.98 
10 700.05 
K’ = 1-K” 2 K’ = 2-K” =3 K’ = 3-K" =4 
R(J) P(J) R(J) P(J) R(J) 
2 | 62083 .39 62029.62 61994.21 
3 095.64 023 .85 969.97 | 62111.89 62040 .30 61986 .57 
4 106.48 015.99 944 .37 122.19 032.52 960.69 | 62138.89 62048 .63 61977 .45 
5 116.14 007 .61 917.12 131.42 023.18 933 .17 146.95 039 .73 949.61 
6 125.06 61997 .37 888.75 138.89 013.01 904.34 154.25 028.58 920.23 
7 133.45 986 .57 859.21 000.58 874.22 161.12 889 .50 
8 141.85 976.13 829 .34 61989 .36 843.43 166.91 858.23 
9 965.12 799 .07 | 812.45 172.99 826.28 
10 769.13 | 781.09 794.17 
11 750.31 761.83 
12 719.87 
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TABLE (Continued) 
Vacuum wave numbers and assignments for the lines the NH; bands 
K'=4-K" =5 K’ = 5-K” =6 K’ = 6-K"” =7 K' =7-K" =8 
175.61 905.35 190.65 047.14 921.93 
180.74 873.36 194.83 032.52 888.75 905.35 
840.43 198.16 017.51 211.83 870.19 
806.77 201.52 001.44 820.40 834.65 849.76 
785.98 
752.03 
P(J) P(J) 
866.43 
band 
62973.21 
62962.77 969.54 
941.32 
62920.37 923.83 
981.96 992.18 873.17 884.37 782.43 
980.43 000.32 865.06 745.98 
62822.74 831.50 
63001 .02 
012.18 025.36 940.83 873.17 034.80 960.29 885.29 
015.31 913.35 828.53 926.23 841.71 
012.18 896.57 794.66 034.80 935.72 817.19 
758.52 771.07 
719.88 731.86 905.53 738.72 
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TABLE (Continued) 
Vacuum wave numbers and assignments for the lines the NH; bands 


837.75 
798.75 
772.98 814.16 
P(J) 
62920.07 
band 
965.18 
P(J) R(J) R(J) P(J) 


63994 .75 


63981 .98 987.21 63889 .93 896.44 63816 
975.38 63863 .15 980.14 874.67 
962.85 841.61 973.28 847 .93 

678.68 952.93 


63989 995.98 63915.20 
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TABLE (Continued) 
Vacuum wave numbers and assignments for the lines the NH; bands 
745.22 755.67 895.72 762.98 
668.14 
045.75 964.25 899.45 981.26 916.51 000.16 
761.23 910.16 778.67 
882.00 
685.81 
956.08 867 .00 006.46 891.78 
752.02 804.13 
P(J) P(J) P(J) 
874.67 
9 | | 63887 .25 
831.86 63867 .00 
band 
015.99 010.48 


66977 .93 
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TABLE (Concluded) 
Vacuum wave numbers and assignments for the lines the NH; bands 
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67080.93 
67076.84 
030.91 026.38 930.03 847.88 930.03 817.88 
002.59 888.34 791.65 791.65 
965.53 840.77 730.77 
788.48 
67085.48 
002.59 918.20 820.99 
845.02 864.80 740.16 760.79 
778.39 910.73 800.77 
733.99 
795.48 840.77 892.47 
810.59 
66951.14 
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found every band, which could only accounted for assuming they 
were lines the missing branch, displaced from the expected positions. The 
displacement the lines the branches varies with the value and must 
levels shown Fig. The explanation this observation will be.discussed 
later. 

order determine the rotational constants the excited electronic state, 
the energy levels the excited state were determined. This was done first 
calculating the energy levels the lower state from the data Benedict, 
Plyler, and Tidwell (1958) and then adding these lower state energies the 
frequencies appropriate observed lines. The constants the equations 
describing the rotational levels were then adjusted give best fit these 
upper state energy levels. This method evaluating molecular constants 
more suitable than one using combination differences the spectrum has many 
overlapped lines. 

attempt was made fit the rotational energy levels the upper state 
the equations usually applicable the degenerate state symmetric top, 
that 


—D yx. J(J+ 1)K?—Dr, 


Two major difficulties were encountered. Firstly, mentioned above, the 
level the observed levels which give rise branch lines are not 
the positions calculated from the levels giving rise and branch lines. 
The displacements from the expected positions were great This 
effect can best thought splitting between the and levels for 
although fact all levels are missing ammonia. The second 
major difficulty trying fit equation (1) the observed levels was that 
the energy difference between the and levels for given which 
the equation should and independent was found very 

Both the above difficulties can resolved assumed that the 
rotational levels are perturbed Coriolis interaction the type described 
giant /-type doubling Garing, Nielsen, and Rao (1959). They have shown 
that interaction between and type vibrational state can modify 
the rotational energies the latter such that for the and levels with 
K=1 


: 
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and for all other levels 


16(K+1)[(1—¢,)C,—B,] 


Here the equations Garing, Nielsen, and Rao have been modified letting 
and dropping small term 6’. The centrifugal stretching terms have 
been omitted but assumed that these will have the same form 
equation (1). 

Equations (2) and (3) introduce only one new constant beyond those given 
equation (1). This constant can determined the splitting 
the levels which according equation (2) equal Thus 
test whether not equations (2) and (8) are applicable the observed 
rotational levels determine from equation (2) and insert this value 
into equation (3) and see the last term will give the observed dependence 
the splitting between the and levels. This calculation was carried out 
and was indeed found that most this dependence could accounted for 
equation (3). This result shown Fig. where the splitting between the 


Fic. The open circles show the separation the and levels the same value 
for the state. The triangles show the same separation after adding the last 
term equation (3) the value each energy level. 


the correction has been applied. That there still should remain small depen- 
dence this splitting not unexpected (see Benedict, Plyler, and Tidwell 
1958, equation 3a) and must conclude that the equations (2) and 
account for the positions the observed levels. 

the final determination the rotational constants, the value for 
each vibrational level was determined from the level. With this value 
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and approximate values and the last term equation (3) was 
calculated and added the observed energies the rotational levels. With 
these corrections added, the rotational levels were represented the equation 


and were adjusted give the best fit the observed levels. The final 
constants vo, and are given Table II. The accuracy 


the centrifugal distortion constants given this table quite low and force 
constants determined from them might far from correct. 


TABLE 
Constants for the upper state NH; 

B, G $v vo q X1074 Dg 
(10.29)* (5.17)* 59224.7 

9.73 0.873 60123.8 0.758 —58 +91 

9.21 0.859 61054.7 0.550 —28 +49 

8.73 5.37 0.842 62014.9 0.529 +58 

8.27 0.839 0.428 —44 +42 
7.87 5.46 0.832 63997 0.258 +67 


*Values obtained by extrapolation: 
Bv = 10.29—.5777+.018v? 
Cv = 5.21+.0430 


The agreement between the values the rotational energy levels calculated 
from equations (2) and (3) with the constants given Table and the 
observed values not good. Differences between the observed and 
calculated values are common and few differences great are 
found. The largest difference seems occur terms with low and high 
but equations (2) and (3) there are eight constants which can adjusted 
fit the rotational levels any one vibrational state and difficult 
determine which levels actually fail fit the equations and which merely 
appear fail because the particular choice constants. The constant 
particularly troublesome since seems vary greatly with and single 
meaningful value could obtained. Though better fit the observed levels 
could obtained introducing additional terms, the relatively low accuracy 
the observations would make the new constants doubtful value. 

The vibrational energy levels the upper state can represented the 
equation 


The positions the band origins differ from the values measured Duncan 


One the unusual features the upper state the large value the 
interaction constant a3. The value 0.58 found for 


| 
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this constant can compared with value 0.015 for the lower 
state. Though true that these values are not exactly comparable 
difficult believe that the change the geometry and the force constants 
could change factor about 40. 

The value corresponds N—H bond distance 1.027 instead 
apparent that the bond distance little different from the ground 
state value 1.011 (Benedict and Plyler 1957). 

Section III the structure band system resulting from particular 
transition was discussed and three essential characteristics were given. 
shall now consider whether not these characteristics have been found 
exist the observed system. 

(1) There doubt that the bands are perpendicular bands. Both the 
number subbands and the intensity the branches the subbands can 
explained only assuming that the bands are the perpendicular type. Thus 
the upper state levels must degenerate (£) states. 

The question whether this degeneracy the upper state arises from 
degeneracy the electronic wave function the vibrational wave function 
requires some consideration. Since the upper state planar and 
symmetry (see the following section) the possible symmetry types the 
electronic wave function are and Since the electronic 
selection rules are determined the symmetry the lower state 
state would give allowed transition but the bands would the 
parallel type. The fact that perpendicular bands are observed therefore rules 
out this possibility. state can not combine with the ground state. 
forbidden transition Aj’ state can made allowed and can give 
rise perpendicular bands one quantum type vibration excited. 
Such band system would weak one whereas the 1600 system strong 
and according Walsh and Warsop (1960) the first member Rydberg 
series. therefore appears most probable that the band system allowed and 
has electronically degenerate excited state. 

(2) Table the subbands have every other line missing 
all their branches. Since these subbands are strong there doubt about 
this conclusion. Furthermore was found that the 0-0, 2-0, 4-0, and 8-0 
bands the even lines are present while the 3-0, and levels the odd 
lines are present. This alternation, according the argument given previously, 
arises from the fact that the upper state levels combine alternately with the 
symmetric and antisymmetric lower states and provides with proof that 
the upper state progression the vibration. 

Having established that the upper state progression the out-of-plane 
vibration can decide the question the planarity the molecule 
this state. Here the vibrational spectrum gives the best evidence. molecule 
such ammonia, the structure could changed progressively from 
planar pyramidal type, the energy the lowest level out-of-plane vibra- 
tion would drop very rapidly and finally become the antisymmetric com- 
ponent the ground state. the ground electronic state ammonia, for 
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example, this vibrational level has dropped until only 0.66 above 
the ground state. the upper state the 1600 band system the molecule 
were far from planar then successive pairs the and levels would 
coalesce and the observed alternation band structure could not occur. 
Therefore there remains only the question whether not the molecule 
slightly non-planar. the molecule were slightly non-planar would not 
expect successive vibrational quanta form regular series and, particular, 
would expect the second vibrational quantum larger than the first 
third. The fact that the first few vibrational levels form regular pro- 
gression must taken proof that the NH; molecule planar the excited 
state. 

(3) order show that the upper electronic state rather than 
symmetry, necessary show that the level combines with the 
antisymmetric level the lower state. Though there doubt that the 
level have numbered does combine with the level, mentioned 
previously, there can some question the correct numbering the 
vibrational levels. From the vibrational spectrum alone there proof that 
the numbering have used the correct numbering, comes from the rota- 
tional constants. For planar symmetric top molecule 2C, and expect 
very nearly equal With the vibrational numbering have 
used, 10.29 and 10.34. this numbering increased one then, 
extrapolating, find that 10.89 and 10.26. Thus appears 
that the rotational constants give strong evidence support the vibrational 
numbering which have used. must admitted, however, that have 
observed very large and possible that the other values are equally 
large. this so, then possible that the agreement between and 
which have observed fortuitous. 

From the above discussion appears that, though some the character- 
istics the band system have not been established with certainty, most 
the major characteristics can deduced with confidence and that these are 
agreement with the transition discussed Section II. 

addition the progression described above, weaker progression with 
about the same spacing has been observed. One these bands shown 
Fig. Bands have been observed 64687, 65639, 66607, 67602, and 68598 
The first observed member this progression the low-frequency 
side the band and further members occur between succeeding bands 
the main progression. Each band characterized very diffuse absorption 
maximum the high-frequency side and about strong diffuse lines Jower 
frequencies. still lower frequencies there are number somewhat sharper 
lines which show regularity except the first band the progression where 
perhaps accidental coincidences series strong broad lines appears. One 
particularly strong line which occurs near the center each band and which 
appears similar all the bands was measured and given above the 
position the band. 
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attempt was made analyze the bands but without success.* The 
diffuseness the lines prevented accurate measurements and the number 
lines observed considerably less than expected for either parallel per- 
pendicular band. possible that the extent the predissociation depends 
the and value the levels and that certain the levels are diffuse 
that they not give rise recognizable lines. 

the weaker progression the main progression the bands fade out 
longer wavelengths. Hence not possible determine whether not the 
first observed band the first member the progression. Thus, without 
analysis the bands and without even true vibrational numbering the 
bands impossible determine whether they arise from excited electronic 
state different from that the main progression from progression 
the same electronic state. consider electronic—vibrational interaction 
discussed below, also possible that the bands arise from progres- 
sion where the degenerate N—H stretching vibration. 


DISCUSSION 


Walsh and Warsop have considered the various electronic states ammonia 
and have assigned the 1600 system transition which electron 
excited from non-bonding orbital Rydberg orbital. They have also 
pointed out that the latter orbital probably correlates with p-type atomic 
orbital and that the excited state the molecule should become planar (or 
nearly so) and that the N—H bond distance should change very little. The 
results the analysis the bands are largely agreement with this assign- 
ment and only the assignment the upper state the rather than the 
state resulting from the configuration must this interpretation 
altered. 

Jahn and Teller have pointed out that molecule such NH; excited 
degenerate electronic state the molecule may unstable with respect 
motions the nuclei which correspond the degenerate modes vibration 
the molecule. Under these conditions the molecule assumes some symmetry 
lower than and the degeneracy the electronic state will split. 
detailed analysis the spectrum any molecule showing effect 
has yet been carried out and the 1600 bands ammonia, which have 
degenerate upper state, were examined see they showed evidence this 
effect. has been discussed earlier, the rotational analysis showed that the 
energy levels the upper state are those symmetric top modified only 
Coriolis perturbation the type described Garing, Nielsen, and Rao. 
true that the rotational energy levels can not reproduced very precisely 
the equations but, nevertheless, they give clear evidence molecular 
distortion. has been pointed out Longuet-Higgins, Opik, Pryce, and 
Sack (1958) and Moffitt and Liehr (1957) that one the major effects 
Jahn-Teller interaction the appearance one more the degenerate 


*These bands have now been shown parallel bands. report the analysis this 
band system will published the near future. 


q 
q 
q 
| 
a 
q 
i 
4 4 


— 


500 CANADIAN JOURNAL PHYSICS. VOL. 39, 1961 


modes vibration the absorption spectrum. have seen there 
clear-cut evidence for this, though just possible that the weak progression 
described previously does involve the degenerate vibration. Thus must 
conclude that the 1600 bands yield information regarding the Jahn- 
Teller effect. 

The separation between the and levels with and the depen- 
dence the separation between the and levels given have been 
interpreted above ‘giant /-type doubling’ the type described Garing, 
Nielsen, and Rao. This effect the result interaction between and 
type vibrational level which are separated only small energy dif- 
ference. the excited state ammonia, have knowledge the 
levels but the only ones likely close the observed state are those 
resulting from the type in-plane bending vibration. possible that this 
in-plane vibration may have much lower frequency than the corresponding 
one the ground state where lies 1626 also appears possible that 
the observed effect electronic rather than vibrational origin and that the 
interaction between the observed state and unknown non-degenerate 
electronic state. The question whether the interaction electronic 
level, which should show the effects electronic but not vibrational 
interaction. Unfortunately the few lines the band which can measured 
and assigned with certainty have, their upper states, rotational levels which 
are not particularly sensitive the value Finally can not rule out the 
possibility that slight distortion’ the molecule rather 
than ‘giant doubling’ which responsible for the observed effect. 

There remain several aspects the upper state, discussed here, which could 
investigated further. may possible theoretically obtain some 
understanding the large value and the particular value which 
has been found. The spectrum has part been photographed and 
this remains analyzed. Also the spectra the mixed species NH2D and 
have not been analyzed. These two molecules are particularly interest- 
ing since they are symmetry the sense that the electronic wave 
function will have threefold symmetry but are very asymmetric tops that 
the three moments inertia are widely different from each other. study 
the interaction between electronic vibrational and rotational motion these 
molecules would considerable interest. 


The authors are indebted Dr. Herzberg, who read the manuscript 
and made numerous suggestions, and Mr. Alberti for taking the plates. 
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RADIO-STAR SCINTILLATIONS AND THE AURORAL ZONE! 


ABSTRACT 


continuous series observations scintillations the radio star, Cassiopeia 
carried out Saskatoon frequency near Mc/s over period nearly 
years has been analyzed. The altitude-angle dependence the scintillations 
was very strong 1955 but weak 1958. This behavior suggests that the 
scintillations are not produced uniform layer the atmosphere. seems 
more likely that the scintillations arise most strongly regions the atmos- 
phere closely associated with the auroral zone and that these regions migrate 
southward during years intense sunspot activity. 


INTRODUCTION 


When radio star observed frequencies the very high frequency 
range the signal strength fluctuates random manner. These scintillations 
have been shown arise the ionosphere (Little 1951; Hewish 1952a), but 
the nature and precise location the ionic inhomogeneities that cause the 
scintillations are uncertain. The earlier investigations (see, for example, 
Hewish 19526) found that the scintillations showed strong solar-diurnal 
variation, with maximum occurrence during the night hours. Later, when 
regular observations were begun Canada 1953 (Hartz 1955) was 
noticed immediately that the solar-diurnal variation was not particularly 
strong. This fact was attributed the proximity the auroral zone. One 
possibility that new mechanism becomes dominant the production 
scintillations near the auroral zone, but convincing evidence this new 
mechanism has still presented. The usual assumption has been that the 
scintillations are generated the ionosphere heights 300 400 and 
that this height range applies equally well the scintillations observed near 
the auroral zone and those observed lower latitudes. 

The principal radio source for scintillation studies Canada has been the 
very intense radio star, Cassiopeia which circumpolar for all Canadian 
stations (Dec. +58.5°). With very simple equipment possible observe 
this source continuously, even though greatly changes its position the 
sky during the course day. For example, Saskatoon (latitude 
52.1° N.) the source about 6.4° north the zenith upper transit and 
about 20.6° above the northern horizon lower transit. The necessary equip- 
ment has been described Hartz (1955) and Reid (1957). 

The variation occurrence scintillations with altitude angle (or with 
sidereal time) has been studied Hartz (1955) for Ottawa and Reid 
(1957) for Ottawa and Saskatoon. Hartz observed that scintillations occurred 
much more often when the source was low altitude angles. interpreted 
this effect being due the longer transmission path through the ionosphere 
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for the smaller altitude angles. Using this interpretation, should possible, 
least principle, determine the height the active region studying 
detail the dependence the scintillations upon altitude angle. Hartz 
attempted this determination but his results were not conclusive. Reid noticed 
asymmetry about the time lower transit the scintillation rates observed 
Saskatoon, and subsequently detected similar asymmetry (in the opposite 
sense) some records obtained Ottawa. deduced that the irregularities 
responsible for the scintillations were greatly elongated along the magnetic 
field direction and tended drift along circular paths centered the magnetic 
dip pole. 

Still more recently, Hartz (1959) has pointed out that many features 
the observed sidereal-time variations might explained equally well 
postulating that the ionosphere becomes progressively more disturbed the 
auroral zone approached, and that the occurrence scintillations depends 
only the degree disturbance that exists that part the ionosphere 
traversed the radio path from the radio star the receiver. Since this 
path typically traverses the ionosphere points separated several hundreds 
kilometers for upper and lower transits, quite reasonable gradient 
ionospheric disturbance with distance from the center the auroral zone 
would required explain the observations. support this suggestion 
Hartz refers unpublished results obtained Ryan and Harrower (at 
Kingston, Ontario) which apparently showed that the scintillation behavior 
radio sources the southern sky was quite different from that sources 
the northern sky. Hvatum and Orhaug (1957) advanced somewhat 
similar suggestion explain rate that they observed 
Gothenburg, Sweden. their case, the change appeared suddenly the 
transmission path entered left region the ionosphere which the authors 
associated with the auroral zone. Also, the observations Tuominen, Kataja, 
and Riihimaa (1960), made Finland, indicate relatively small variation 
scintillation rate with altitude angle the source. Although they suggest 
that this effect due systematic motions can also ascribed the 
presence disturbed (auroral) regions high altitude angles. The purpose 
the present paper draw attention the fact that the nature the 
sidereal-time variation scintillations, observed Saskatoon, changed very 
considerably over period years. This change such rule out the 
variation ionospheric path-length the dominant factor the sidereal-time 
variation. seems more likely that disturbed region the ionosphere 
the neighborhood the auroral zone gives rise the more intense scintilla- 
tions. According this interpretation the Saskatoon data indicate movement 
the auroral zone toward lower latitudes during years high sunspot 
number, has been deduced previously Davies (1948) from optical 
evidence. 


ANALYSIS AND RESULTS 


For studies radio-star scintillations customary record the received 
signal strength means moving chart recorder. the analysis the 
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charts scintillation amplitude index and scintillation rate index are assigned 
each hour. Ideally, the amplitude index proportional the fractional 
fluctuation the signal amplitude, and the rate index the number 
fluctuation peaks that occur given time interval. practice, because 
the limited resolution the record, the two indices are assigned each 
hour the basis the visual appearance the record. relatively small 
number classes used for each index. Two analysts may differ the class 
assignment particular record, but there reason suspect that the 
use absolute scale would greatly improve the quality the data. Typical 
examples chart records corresponding different amplitude indices have 
been given Hartz (1959) and different rate indices Reid (1957). 

Saskatoon, the scintillation recorder was put into operation frequency 
about Mc/s, 1954 and consistent set records available covering 
the period from January 1955 through the autumn 1958. These records 
were scaled routine basis they were produced. Recently, the course 
study some the later data was noticed that the sidereal-time 
behavior was quite different from that found earlier Costain (1955) and 
that reported Reid (1957). Since several different analysts had carried 
out the scaling the records for various parts the 4-year period was 
first necessary verify the consistency the chart analysis. This was done 
having single analyst back the original records and scale sample 
periods overlapping each the original changes analyst. Also, randomly 
picked samples record were scaled and the indices compared with those 
derived earlier. While occasional differences (of one unit) occurred between 
the new and old indices there were systematic differences that could affect 
the averaged results. Also, was apparent from the records that while changes 
equipment sensitivity had occurred, the reliability the derived indices 
had not been affected thereby. 

The mean monthly index for each hour sidereal time after upper transit 
the radio source was derived. These indices were used compute 6-month 
running means. This computation resulted group monthly sets, each 
set consisting individual numbers, and each number representing the 
index for single hour sidereal time averaged over 6-month period. Each 
set was then subjected harmonic analysis determine the phase and 
amplitude the harmonic component representing the 24-hour (sidereal) com- 
ponent the variation. The results are shown Figs. and harmonic 
dials the type described Chapman and Bartels (1940). These dials show 
the time maximum angular position measured clockwise from the 
vertical, and amplitude distance from the origin, the 24-hour harmonic 
component the time variation. Figure shows the values for the rate 
index variation plotted the harmonic dial with the successive values joined 
straight line. Figure gives the same information for the amplitude 
index variation. 

Figures and give information concerning the fractional variation 
the index values within each monthly set. order indicate the long term 
trend the mean monthly indices relation the harmonic variations shown 
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Harmonic dial showing the phase and amplitude the 24-hour sidereal component 
the variation scintillation rate index. Time measured from upper transit the source. 


Fic. dial for the 24-hour variation amplitude index covering the same 
period that Fig. 
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Figs. and the data were treated the manner indicated Figs. and 
each monthly set the values corresponding the hours before and 
hours after upper transit were averaged produce mean 
The remaining values were averaged produce mean 
transit These two values are plotted against time (in months) Fig. 
for the rate indices and Fig. for the amplitude indices. The smoothed 
sunspot numbers for the same period are plotted also each these figures. 
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Fic. the 6-month means the scintillation rate index for the hours 
near upper transit, and for the hours near lower transit, together with the smoothed 


observed sunspot numbers. 
Fic, Graph showing the same information that Fig. for the amplitude indices. 


DISCUSSION 


examination Fig. shows interesting trend the sidereal-time 
variation the rate indices over the 46-month period observation. The 
amplitude the variation decreases from one index unit the beginning 
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the period small value the end. The annual loops formed the harmonic 
dial are obviously due solar-diurnal variation. The phase this variation 
indicates that the maximum occurs near midnight (in general agreement with 
the British observations) but the amplitude very much less than the 
sidereal-time variation. 

The phase the main (sidereal) variation, shown Fig. indicates 
that the maximum occurs about hours after upper transit about hour 
after lower transit. This displacement the maximum from the time lower 
transit just the asymmetry found Reid different method, using 
some these same data. interesting note that the period chosen 
Reid (starting April 1955) happened the year when the sidereal 
variation was strongest the Saskatoon observations. 1958 there was 
little variation the scintillation rate with position the sky. Figure 
shows that the effect can not ascribed general decrease scintillation 
activity. The mean level activity changed very little over the 4-year period, 
but the transit index’’ decreased the value the ‘‘upper transit 

The behavior the amplitude index, indicated Figs. and only 
slightly different from that the rate index. From Fig. may seen 
that the sidereal-time variation disappeared more rapidly for the amplitude 
index than for the rate index, being virtually complete about midway through 
1956. From Fig. the predominant effect this change seems have been 

seems most likely that the observed behavior the scintillation indices 
linked some way solar activity. may noted that the sidereal 
(altitude angle) variation was very prominent for periods when the sunspot 
number was appreciably less than 100, and very small for periods when the 
sunspot number was much above 100. Unfortunately there seems 
way which this link can established firmly from the present data. much 
longer series observations, extending over about years, would required. 

the other hand, the present data seem rule out the possibility that 
the sidereal variation due the variation path-length relatively 
uniform layer. This mechanism, considered some detail Hartz (1955) 
and Reid (1957), would lead decrease sidereal variation only for periods 
when the mean scintillation activity was low. 

The recurring suggestion that the scintillation activity Saskatoon 
directly connected with the auroral zone gains some support from the present 
measurements. The observations could explained southward migration 
the auroral zone during periods intense auroral (and solar) activity. Such 
migration was reported Davies (1948) based large number visual 
observations aurora obtained various times the period 1932 1947. 
deduced that the southern limit the auroral zone nearly overhead 
Saskatoon during sunspot minima and well the south Saskatoon during 
sunspot maxima. Although not clearly established, the northern edge 
the auroral zone seems subject similar migration. The variation 
auroral activity with geomagnetic latitude reasonably well known, and 
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there some evidence that the variation radio-auroral activity similar 
(Forsyth, Green, and Mah 1960). The radio-aurora (detected radar 
radio-scatter measurements) generally located heights near the 100-kilo- 
meter level and usually assumed that the scintillations are produced 
heights near above the 300-kilometer level, although the evidence supporting 
this assumption, least for the Canadian observations, scanty 
suspect. can argued that, whatever the nature other sources involved, 
radio-aurora should contribute the production scintillations. 

any case, the assumption geographic distribution scintillation- 
producing regions, similar that aurora, leads acceptable explanation 
the trends shown Figs. principle, possible use the details 
Davies’ description the auroral zone migration and the data Figs. 
and deduce the height the active region. However, the details the 
migration and the causative link between aurora and scintillations, not seem 
well enough established justify such analysis this time. 

Scintillation data for frequency Mc/s obtained Jodrell Bank 
England (latitude 53° N.) and covering the same period that the present 
study have recently appeared the literature. The detailed analysis these 
data, reported Chivers (1960), shows many interesting contrasts between the 
Jodrell Bank and the Saskatoon scintillations. particular, the solar-diurnal 
variation both scintillation rate and scintillation amplitude Jodrell Bank 
was more pronounced than the sidereal-time variation, for the whole period 
observation. The sidereal-time variation, which was weak 1955, became 
stronger 1958. This change occurred earlier the period and more abruptly 
for the scintillation amplitudes than for the scintillation rates. appears 
that similar changes were taking place simultaneously Jodrell Bank and 
Saskatoon, but opposite senses. The Jodrell Bank results are consistent 
with the suggested migration scintillation-producing auroral zone, 
assumed that the zone too far north observed from Jodrell Bank during 
years low sunspot activity and comes into view low elevation angles 
during years high sunspot activity. Further support for this picture pro- 
vided the fact that, contrast the Saskatoon results, the scintillation 
activity Jodrell Bank increased markedly with increasing sunspot number. 
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REPRESENTATIONS THE UNITARY GROUP 
AND WAVE FUNCTIONS! 


VENABLES 


ABSTRACT 
number wave functions besides the spherical harmonics are obtainable 


from the irreducible representations the two-dimensional unitary group. 
INTRODUCTION 


The irreducible representations the rotational group obtained from 
the two-dimensional unitary group are well known. Their use the 
quantum theory angular momentum has been extensively treated (Rose 
1957; Sharp 1957; Wigner 1959). Where are Euler angles 


(¢, 6, v) = 
Here 
2j+m—m’'—2q  2g—m+m’ 6 

spherical harmonics. The relation involves constants and perhaps taking 
the conjugate function, depending the convention use the particular 
author (Rose 1957; Sharp 1957; Wigner 1959). 

Besides the spherical harmonics seems worth while point out that 
there are wave functions, which can obtained from the unitary repre- 
sentation. These functions are among the simplest and most successful 
quantum They are, where constants are omitted (Bethe and 
Salpeter 1957; Schiff 1949), 

(i) hydrogen radial wave function 
(ii) hydrogen parabolic wave function 


(iii) harmonic oscillator wave function 


(x). 


Each function can derived from entirety, except for minor 
adjustment required the radial case. 


1Manuscript received December 14, 1960. 
Contribution from the University Alberta, Edmonton, Alberta. 


Can, J. Phys. Vol. 39 (1961) 


510 


a 
— 


VENABLES: UNITARY GROUP 511 


expressible terms Jacobi polynomials (Rose 1957; Sharp 1957) 


(0) = 


(1) dium’ (x) = 
LIMITING FORMS 


The Jacobi polynomials yield the Laguerre and Hermite polynomials under 
special conditions. (1959) and Tricomi (1955) give the relation between 
the Jacobi and Laguerre polynomials 


and similar relation holds for the Hermite polynomials 


These may obtained number ways. One method employ 
the differential equation (Poole 1960) (x), 


For the Laguerre equation change the variable and let 


For the Hermite equation change the variable and let 
followed 
The three wave functions will now considered turn. For the first case 
choose for the indices 


+2nH 


and let 


| 
t 
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This gives (1) 


which the limit (2) becomes 
The function the brackets constitutes the radial only the 


resence the factor which prevents precise agreement between 


and (i). 
order deal with the parabolic case the following indices are chosen 


and 
Taking the limit before gives 


2 —é/2 


similar function this can found with the variable, and the 
product the two functions will give (ii). 
For the harmonic oscillator (iii) only necessary take the limit (3) 


give 


Here have used (1) 
=n, 


and changed the variable x/+/m. 


DISCUSSION 


Instead the Jacobi equation (4) the equation (Whit- 
taker and Watson 1927) may found and used. procedures 
similar those already employed the wave equations the hydrogen atom 
and the harmonic oscillator can readily obtained. interesting observe 
that the Coulomb potential and the quantum differential operator are derived 
together simultaneously unified way. This seems suggest that further 
examination may interest. 
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ROOM TEMPERATURE STABILIZATION 
RADIATION-PRODUCED FREE RADICALS 
BARBITURIC ACIDS! 


ABSTRACT 


Experimental evidence has shown that boric acid glass may used 
stabilize room temperature organic free radicals produced radiation 
number barbituric acid derivatives. The method and the experimental results 
are discussed. 


INTRODUCTION 


The immediate effect radiation most organic compounds generally 
believed the production free radicals which initiate several unexpected 
reactions. Due the tremendous energy outbalance, difficult predict 
where and how bond ruptures will effected. Fortunately, experimental 
evidence indicates that radiation damage many organic compounds (Smaller 
and Matheson 1958; Gordy, Ard, and Shields 1955a; Meara and Shaw 1957) 
and even complex biological systems (Conger and Randolph 1959) follow 
general pattern where most important processes may observed, modified, 
and partially controlled (Patt 1953; Gordy and Miyagawa 1960). The under- 
standing radiation damage mechanism presupposes the knowledge the 
types free radicals being formed and study their properties. The method 
paramagnetic resonance has great potentialities for the identification 
radicals and has been used successfully number workers (Ingram 1955; 
Shoolery and Weaver 1955; Whiffen 1958). One the main difficulties that 
some the free radicals are very reactive character and must stabilized 
order obtain concentrations which are detectable with current instru- 
mental sensitivities. Low temperature offers immediate solution this 
problem. Free radicals are formed rigid medium which interferes with 
their recombination and allows their detection. However, when biological 
applications are considered, low temperature results may misleading. 
far from being proved that reactions which occur boiling helium nitrogen 
temperature, where molecular vibrational and rotational energies are much 
reduced and molecular resonance conditions altered, have some significance 
room temperature. Radiation damage mechanism may totally different 
more less exemplified Gordy, Ard, and Shields with irradiated 
acetic acid studied various temperatures. Room temperature irradiation 
and measurement under nitrogen helium atmosphere has also been used 
extensively various research groups. However, oxidation the unpaired 
electron pronounced with some free radical species that weak para- 
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magnetic resonances are detected only for short time. Further, when gaseous 
by-products are produced during irradiation, recombination the free radicals 
occurs independently the nature the atmosphere above the sample. 

Some two years ago, initial experiments were made study radiation 
effect series 5,5-disubstituted barbituric acids. resonance signals 
were observed when the samples were irradiated air; only two cases yielded 
some information when kept under nitrogen atmosphere. Since was intended 
carry out this research room temperature, search was made how 
stabilize the free radicals produced radiation these acids. The experi- 
mental material was dispersed number substances solid room 
temperature, irradiated, and tested for paramagnetic resonance absorption. 
Interesting results were obtained with boric acid glass. the purpose 
this paper present the experimental method and data obtained over 
15-month period. 


APPARATUS 


Figure shows block diagram the electron paramagnetic resonance 
spectrometer used this investigation. conventional design except 
that has been assembled using standard commercial units 


STRAND LAB. 
MODULATOR 
STRAND LAB. 


PRECISION GAUSSMETER 
MODEL M-Z 
NUCLEAR MAGNETICS CORP. 


VARIAN 
x-13 

KLYSTRON 

OiL BATH 


MODEL 


- 
CAVITY RESONATOR LOCK-IN AMPLIFIER 


MODULATION COILS 


MAGNET 


MAGNET POWER SUPPLY 
VARIAN MODEL V~-2200A 


PRECESSION FIELD SCANNING 
VARIAN MODEL V-4280A 


WAVEMETER 
PRO 559A 


OSCILLOSCOPE 
DUMONT TYPE 3044 


UNIVERSAL KLYSTRON 
POWER SUPPLY 
PRD 801A 


VARIAN 
MODEL G-I0 
RECORDER 


SORENSEN 
FRLO 750 
VOLTAGE REGULATOR 


Fic. diagram the electron paramagnetic resonance spectrometer. Nominal 
operating frequency 9500 megacycles per second. Modulation frequency kilocycles per 
second. 


blocks. operates nominal frequency 9500 megacycles per second. The 
transmission type microwave cavity resonator rectangular cell, 
silver-plated the inside, fitted with two coupling irises, 0.25 in. diameter. 
The Zeeman modulation coils are connected series. Each coil was made 
500 turns #32 H.F. wire mounted 0.50-in. spool, has resistance ohms, 
and can carry maximum current 0.32 amp. The modulation frequency 
kilocycles per second. The air gap between the two cylindrical pole pieces 
inch. The homogeneity the magnetic field the sample position 
the order 0.2%, determined with proton resonance meter. The klystron 


— 
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constant temperature oil bath. date, microwave frequency 
stabilization has been required, the frequency being stable enough for the type 
measurement being made. The spectrometer similar the one used by. 
and Shaw (1957) and its sensitivity the same order magnitude, 
unpaired electrons 100-mg sample. The information obtained 
the first derivative the microwave energy absorption the sample. 
Figure shows the signal observed from diphenyl picryl hydrazyl (D.P.P.H.) 
the dry state. The distance between points maximum slope slightly 
under oersteds, line width which compares quite well with published 
information (Gordy, Smith, and Trambarulo 1953). 


= 
=== 
= 


Fic. Paramagnetic resonance absorption signal from diphenyl picryl hydrazylin dry 
powder form. First derivative curve. 2.0036. Markers are oersteds apart. Distance 
between points maximum slope slightly under oersteds. 


Some the results reported here were obtained using earlier design 
the electron paramagnetic resonance spectrometer. The Zeeman modulation 
frequency was generated Krohn-Hite ultra-low frequency oscillator 
operating 700 cycles per second. The output the oscillator was amplified 
20-watt MacIntosh power amplifier and connected the modulation 
coils. The crystal detector current was fed Browning standing wave ratio 
amplifier, model TAA-16B, whose output rectifier system had been by-passed. 
This was followed Krohn-Hite model 310-AB variable band-pass filter 


| 
| 
| 
| 
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and synchronous demodulator driven the Zeeman modulation frequency 
oscillator. This earlier design was more flexible than the later one the sense 
that either the actual signal its first second derivative could selected. 
However, the noise level the detecting system was much higher and made 
necessary the reported modifications. 


MATERIAL 


The barbituric acids used this investigation are commercial products 
high purity. Samples from the same lot numbers have already been charac- 
terized members this laboratory (Levi and Hubley 1956; Chatten and 
Levi 1957). Pertinent published information (Chatten and Levi 1957) 
summarized Table for discussion purposes. Levi (1957) has reviewed the 


TABLE 
Melting points disubstituted barbituric acids 
Trade name Chemical name Melting point, 
Alurate 5-Allyl-5-isopropylbarbituric acid 
Amytal acid 
Cyclopal 5-Allyl-5-(1-cyclopenten-2-yl) barbituric acid 139.8-140.3 
Delvinal barbituric acid 
Dial 5,5-Diallylbarbituric acid 172.5-173.0 
Evipal acid 
Luminal 5-Ethyl-5-phenylbarbituric acid 174.1-174.9 
Mebaral acid 115.0-116.0 
Medomin barbituric acid 
Nembutal barbituric acid 128.8-129.4 
Neonal 5-Ethyl-5-n-butylbarbituric acid 
Nostal acid 180.3-181.2 
Phandorn barbituric acid 171.8-172.2 
Rutonal 5-Methyl-5-phenylbarbituric acid 225 
Seconal 5-Allyl-5-(1-methylbutyl) barbituric acid 96.1- 96.9 
Sigmodal barbituric acid 
Veronal 5,5-Diethylbarbituric acid 


chemical structure, synthesis, and nomenclature the barbituric acids. The 
compounds studied have all the same basic structure shown Fig. the 
malonylurea ring which are attached position two groups varying 


Fic. Structural formula 5,5-disubstituted barbituric acids. 


complexity. (In addition, evipal and mebaral have methyl group position 
3.) Barbituric acid itself non-toxic and has hypnotic proper- 
ties. only when both and are replaced alkyl, aryl, alicyclic 
group that the molecule possesses hypnotic, sedative, anaesthetic properties 
(Levi 1957). The heterocyclic system resonating structure; therefore, the 
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substituted barbituric acids are relatively stable molecules except when sub- 
mitted alkaline hydrolysis. Normal physiological processes destroy only 
the side chains and through direct oxidative attack (Maynert and van 
Dyke 1949). thus interest determine the effect radiation this 
type molecule since corresponding changes physiological properties will 
likely ensue. 


IRRADIATION THE SAMPLES 


The samples were irradiated using the facilities kindly offered this 
laboratory the Commercial Products Division Atomic Energy Canada 
Limited. These were 4000-curie Co-60 gamma cell 220, average dose rate 
350,000 rads/hour, and 500-curie Co-60 transfer case irradiator, average 
dose rate 220,000 rads/hour. Some the samples were also irradiated with 
X-ray therapy unit operated (peak) and ma. filter was 
used and the samples were simply taped the exit window frame the 
X-ray tube. this instance, the dose rate has not been measured; however, 
due the large proportion soft radiation, must the order 
few thousand rads per minute. Samples were irradiated either with the Co-60 
sources the X-ray therapy unit; differences were observed far the 
experimental data presented this paper are concerned. 


RIGID DISPERSIVE MEDIUM 


number substances were irradiated and studied for electron para- 
magnetic resonance (E.P.R.) absorption. The irradiation and the measure- 
ments were performed room temperature and air. The samples received 
dose which varied from megarads megarads. The E.P.R. determinations 
were done within one half hour after removal from the radiation source. 
Figure presents the results observed for commercial grade beeswax, sulphur, 
gelatine, paraffin, m.p. 55°-56° and boric acid. evident that, for this 
experiment, the spectrometer had operated maximum sensitivity and 
high amplification detect any possible resonance. 

Beeswax, Fig. 4A, gave signal which disappeared about two hours after 
irradiation. odd line spectrum width 100 oersteds, centered 
g-value equal that D.P.P.H. other attempts were made study 
this signal this observation eliminated beeswax possible dispersive 
medium use this experiment. Sulphur powder, Fig. 4B, was melted, 
allowed solidify rod form, and irradiated. Gelatine, Fig. 4C, was dis- 
solved little water, the solution was allowed set overnight mold 
and irradiated the following day. Paraffin, m.p. Fig. 4D, was 
simply heated, shaped into small cylinders, and irradiated. Boric acid, Fig. 4E, 
was heated slowly until viscous liquid formed; some this liquid was 
drawn out the end pyrex rod and rapidly shaped into cylindrical rod 
rolling between two microscope slides. E.P.R. signals could 
observed from these last four cases, even the highest radiation dose. 
understood that, possibly, with more sensitive spectrometer, resonances 
might have beenobserved room temperature from these irradiated substances. 
However, several experiments showed without doubt that the apparatus 
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Fic. Substances investigated for use rigid dispersive medium. beeswax; sulphur; 
gelatine; paraffin 55°-56° boric acid glass. First derivative curves. Markers are 
100 oersteds apart for and and oersteds for and Vertical arrows indicate the 
position resonance. 
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used this research was unable detect any absorption signal from these 
blanks. For practical reasons, only the paraffin and the boric acid methods 
were applied this study. 

Paraffin known easily form hydrogen bonding; therefore, has always 
been used with circumspection and mainly check the results obtained 
with boric acid. the other hand, boron forms variety complexes with 
many organic compounds (Lappert 1956). However, boron has two isotopes 
non-zero ground spin: B-10, 18.83% abundance, spin 81.17% 
abundance, spin 3/2 (Hollander, Perlman, and Seaborg 1953). The presence 
boron near resonating electron will then indicated further hyper- 
fine splitting the signal into equal intensity septets quartets. These 
types hyperfine splittings have not been observed this study. 


PREPARATION SAMPLES 


The microwave cavity resonator designed accept standard pyrex 
tubes i.d. The samples irradiated air were simply transferred 
one such tube after irradiation and studied. For irradiation and measurement 
under nitrogen atmosphere, the sample was introduced into long tube 
and flushed slowly for few minutes with nitrogen; the tube was then sealed 
during the flushing. Only the small portion the tube containing the sample 
was irradiated, the rest being carefully shielded from the main X-ray 
gamma-ray beam. After irradiation, the sample was shifted tapping 
the non-irradiated end the tube and the E.P.R. measurements were made. 
This technique was used O’Meara and Shaw (1957) and many others. For 
preparation paraffin, 30% mixture volume the dry sample powder 
was made with paraffin held 60° the mixture was molded into small 
cylinders in. long and 3.5mm diameter. The cylinders were irradiated 
and then gently pressed into position the end empty pyrex tube 
for their introduction into the spectrometer cavity. prepare the samples 
boric acid glass, the boric acid was rapidly heated hot plate 
155°+5° 5-ml beaker. Small amounts were used make sure that 
transparent liquid would obtained without the formation white crust. 
The powder sample was added and swiftly mixed with the viscous liquid. 
Twenty thirty per cent concentration was normally used. The small 
cylinders were shaped already outlined above. They were irradiated and 
mounted for study the end pyrex tube with drop household 
cement. 

Referring Table seen that almost half the barbituric acids 
studied were melted during the preparation the boric acid glass. Boric acid 
and the barbituric acids being polar compounds, molecular dispersion 
one into the other does take place. This represents ideal cases where individual 
barbituric acid molecules are trapped within spherical holes the rigid 
matrix. free radical likely diffuse from one site the other and yet 
may assumed that the energy levels the solute molecules are not 
markedly altered the presence the rigid surroundings. For those acids 
which remained solid during the preparation the glass, known that 
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barbituric acids are moderately soluble hot water. Therefore, the water 
resulting from the decomposition boric acid dissolves some the powder 
sample helps breaking the larger crystals into much smaller fragments. 
any case, molecular dispersion this second group compounds melting 
points higher than the boric acid glass also achieved but 
smaller extent. Incidentally, interesting note this time that E.P.R. 
signals recorded from glasses containing samples lower melting points than 
160° were noticeably stronger than those higher melting points. 

Glassy solvents were introduced spectroscopy early 1896 Schmidt 
for the study the phosphorescence ultraviolet irradiated molecules. This 
technique was then used extensively many others (Reid 1958) similar 
application. Lewis, Lipkin, and Magel (1941) showed, for instance, using 
fluorescein dissolved boric acid, that the phosphorescent state could 
characterized its absorption spectrum. Phosphorescence explained the 
presence metastable excited state, the triplet state, which delays the 
return excited molecule its normal ground state. The study the 
triplet state not too distant from the research reported this paper since 
both investigations are concerned with the study the behavior unpaired 
electrons: the former makes use the delayed light emission the excited 
molecule; the latter, the paramagnetic property the radicals which are 
formed. 


EXPERIMENTAL RESULTS AND DISCUSSION INDIVIDUAL CASES 


Alurate 

has already been shown that when boric acid glass paraffin rod 
was irradiated alone, paramagnetic resonance absorption signal could 
detected room temperature and the operating sensitivity the spectro- 
meter used for this research. now pertinent show that with the class 
compounds studied, the effect the solid matrix the shape the recorded 
signals was but negligible. Figure illustrates comparative study which could 
made with alurate. Photographs were obtained from irradiated 
alurate boric acid glass, photographs and from irradiated alurate 
paraffin, and photograph from irradiated alurate dry powder kept under 
nitrogen atmosphere. Concentrating the main structure the signal, 
which evidently odd line spectrum five resolved components, seen 
that the same type signals was obtained from the same compound irradiated 
under three completely different experimental conditions. The slight discre- 
pancies may easily accounted for the following consideration the 
individual signals. 

Figure was recorded minutes after irradiation. The symmetrical 
pattern has non-uniform spacing. Numbering the peaks from left right, 
the two satellites and stand out from the main pattern 
greater distance than the uniform spacing between the central lines. These 
two satellites may thought the normal doublet; however, their 
separation only 115 oersteds precludes this assignment. rather believed 
that Fig. complex pattern resulting from the superposition two 
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Fic. Paramagnetic resonance absorption signals from irradiated alurate. First deriva- 
tive curves. Photographs and are from boric acid glass; and from paraffin; 
and from the dry salt under nitrogen atmosphere. Times after irradiation: min; 
Markers are oersteds apart. Vertical arrows indicate the position D.P.P.H. resonance. 


equally spaced quintets: quintet (1, 7), average spacing oersteds, 
and quintet (2, 6), average spacing oersteds. The time effect 
study this signal indicates such possibility. Figure was obtained 
130 hours after irradiation. noticed that the relative heights peaks 
and have decreased with respect those peaks and Further, 
this decrease more and more pronounced Figs. and which were 
obtained and months after irradiation. the end months after 
irradiation, Fig. 5E, the two satellites are not resolvable from noise. The 
intensity the signal such evidently reduced; however, the contour 
the first derivative plot does not significantly differ from that months 
after irradiation. would therefore seem that during the first months after 
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irradiation, quintet corresponding some reactive free radical disappear 
even when the sample dispersed boric acid glass. Quintet II, associated 
with the more stable free radical which remains after months boric 
acid glass, also obtained from alurate irradiated paraffin, Figs. and 
5G, and the dry powder kept under nitrogen atmosphere, Fig. 5H. The 
signal width oersteds larger with the later two methods prepara- 
tion. However, effective reduction the spin-spin interaction may readily 
explain the narrowing the signal when alurate was dispersed through the 
boric acid glass. 

The presence quintets the resonance absorption pattern from irradiated 
alurate eliminates the possibility that the odd electron was located anywhere 
along the ring structure the molecule. Assuming the hyperfine splitting 
the observed signal was due proton interaction exclusively, the required 
free radicals with four nearly equivalent protons may obtained modifi- 
cation either one the alkyl groups position the malonylurea ring: 
such the removal one CH; the isopropyl group the removal the 
hydrogen atom carbon the allyl group. More complex mechanisms 
could formulated. Specific assignment the two observed quintets 
given free radical species would require more experimental work. never- 
theless useful note that the main effect radiation alurate was centered 
about the two side chains position Further, the boric acid glass preparation 
yielded the same information that obtained from the other two methods 
preparation. The boric acid glass, however, seemed stabilize additional 
free radical species which apparently was too reactive detected with 
the other two methods preparation. Other cases good agreement between 
the results obtained with the boric acid glass and the paraffin methods 
preparation were observed during this study. These will reported their 
appropriate sections. Alurate presented the definite advantage over the other 
cases studied that the results from the above two methods preparation 
could also compared those obtained from the dry powder itself kept 
under nitrogen atmosphere. 


Neonal, Amytal, and Dial 

Neonal, amytal, and dial, when prepared boric acid glass, all gave the 
same E.P.R. signal hour after irradiation. was, presented Figs. 6A, 
6E, and for these compounds the above order, well resolved triplet 
approximately equal intensity. The field separation between adjacent lines 
was oersteds. Three months after irradiation, the shape the signals 
from neonal and amytal, Figs. and 6F, was unaltered, while that from 
dial, Fig. was somewhat less resolved. The three signals were all narrower, 
the line separation being oersteds. This line separation remained constant 
over the next 12-month period shown Figs. 6C, 6G, and for neonal, 
amytal, and dial, respectively. The resolution the lines was preserved and 
their intensities were only slightly reduced. 

The appearance equal intensity triplet leaves doubt the 
location the odd electron. Nitrogen the only atom with spin these 
three molecules. Further, since nitrogen absent from the side structures 
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Fic. resonance absorption signals from irradiated neonal, from 
irradiated amytal, and irradiated dial, First derivative curves. The above 
curves were obtained from boric acid glass except photograph which was from paraffin. 
are oersteds apart for and oersteds for the others. Vertical arrows indicate the 
position D.P.P.H. resonance. 
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the molecules, follows that the odd electron located about one the 
two nitrogen atoms positions and the malonylurea ring. Referring 
Fig. obvious that once atom has been removed radiation, 
for example from position the odd electron would circulate around the 
ring structure and interact with not only the atom position but also 
that position along with its associated atom. other words, this 
were so, one would expect not single triplet but triplet triplets doublets, 
since these are non-equivalent atoms. However, has already been shown 
(Newman 1956) that the barbituric ring not planar. Position acts 
potential barrier which prohibits motion the around the ring. 
Therefore, quite possible that the odd electron compelled spend all 
its time close given atom, either position position 

When neonal was prepared paraffin and irradiated, complex signal was 
obtained, Fig. 6D. believed that this signal unresolved triplet 
doublets. The total width the signal and the separation each doublet 
are about the same those already reported for the boric acid glass prepara- 
tion. Therefore with the paraffin preparation, assumed that the same 
free radical produced; however, coupling hydrogen atom from the 
surrounding matrix promptly takes place. 

The resonance signals presented Fig. match perfectly well those given 
Cole al. (1957, 1958) and Foner al. (1958), who studied the paramag- 
netic resonance nitrogen atoms trapped molecular nitrogen matrix 
4.2° Figure even shows evidence side structure the main 
triplet, which was discussed Cole their second article (1958) 
terms zero field splittings atomic nitrogen. The only discrepancy 
the splitting constant which apparently factor larger when nitrogen 
bound the barbituric acid ring compared when free nitrogen 
atom simply trapped solid matrix. 


Cyclopal 

Irradiated cyclopal paraffin, Fig. 7A, under nitrogen atmosphere, Fig. 
7B, and prepared boric acid glass, Fig. 7C, yielded essentially the same 
E.P.R. signal right after irradiation. was complex pattern about 130 
oersteds, total spread. Four months after irradiation, the signal from the boric 
acid glass preparation, Fig. 7D, showed some changes the relative heights 
the first derivative peaks and definite narrowing. When preparing the 
boric acid glass, was noticed that, few instances, strong reaction 
would take place soon the cyclopal powder was added the molten 
boric acid. Some gas would liberated and the mixture would turn black. 
This would happen when the mixture was temperature slightly over 
165° These black glasses were irradiated. equal intensity triplet 
averaged spacing oersteds, Fig. 7E, was obtained shortly after irradiation 
and was still showing with minor contour changes months after irradiation, 
Fig. 7F. 

The signal from the black glasses undoubtedly due nitrogen nuclear 
coupling with the electron spin. similar free radical assignment that 
the previous group compounds just studied then suggested. believed 
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Fic. Paramagnetic resonance absorption signals from irradiated cyclopal. First deriva- 
tive curves. Photograph was obtained from paraffin, from the salt under nitrogen 
atmosphere, and the others from boric acid glass. Times after irradiation: 
hour; months. Markers are oersteds apart for and oersteds for the others. 
Vertical arrows indicate the position D.P.P.H. resonance. 


that the reaction temperature above 165° accompanied destruction 
the ring structure, 2-cyclopenten, position the cyclopal molecule. 
The mixed signal observed from the unaltered cyclopal molecule, Figs. 6A, 
would then taken include contribution coming from one the 
side chains position possibly the 2-cyclopenten group. Such contri- 
bution assumed removed preparation above 165° leaving 
strong equal intensity triplet due nitrogen interaction. The time effect 
observed the signal from the boric acid glass preparation seems favor 
such explanation. months after irradiation, Figs. and are 
almost alike, indicating that the whatever additional component responsible 
for the initial discrepancies slowly disappearing with time. 


Barbital and Medomin 

triplet total width oersteds was obtained from barbital, Fig. 8A, 
and medomin, Fig. 8B, both prepared boric acid glass. The line spacing 
was uniform about oersteds. appreciable time effect was observed 
the signal from barbital, shown Fig. 8C, months after irradiation, 
and Fig. 8D, months after irradiation. 

Both molecules have common ethyl group one the substituents 
position the barbituric acid ring. would therefore appear that radiation 
removes CH; from this ethyl group leaving radical which has the required 
two equivalent protons. 


Ipral, Sigmodal, Nembutal, and Seconal 

Ipral boric acid glass gave immediately after irradiation, Fig. 9A, 
quintet with uniform line spacing oersteds. Fifteen months after 
irradiation, Fig. 9B, the signal had decreased intensity; however, the 
original quintet with the same line spacing was still resolvable from noise. 
Ipral prepared paraffin, Fig. 9C, yielded complex signal which most 
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Fic. resonance absorption signals from irradiated barbital, and 
irradiated medomin, First derivative curves. The above curves were obtained from boric 
acid glass. Times after irradiation: hours; months; months; months. 
Markers are oersteds apart. Vertical arrows indicate the position D.P.P.H. resonance. 


likely due coupling with the surrounding matrix. interesting note 
that the molecule ipral has position two groups, and isopropyl 
group, which, removal single radiation, could contribute 
the observed signal. The isopropyl group with the loss CH; would give 
rise quintet, average spacing oersteds, already seen for alurate; 
the ethyl group without the end CH; would produce triplet, average spacing 
oersteds, the case barbital and medomin. These two signals being 
centered about the same g-value may added together form the observed 
regular and symmetrical quintet. 

Sigmodal, Fig. 9D, nembutal, Fig. 9G, and seconal, Fig. 91, gave also 
quintet average spacing oersteds immediately after irradiation 
boric acid glass. The shape the signal was not significantly changed after 
and months for sigmodal, Figs. and 9F, and after months for 
nembutal and seconal, Figs. and 9J. Sigmodal, nembutal, and seconal 
have butyl group common position Breakage radiation 
the C—C bond between carbon and would leave the required nearly 
equivalent atoms explain the observed quintets. For nembutal, contri- 
bution from the other substituent position ethyl group, should also 
considered already discussed for ipral. major contribution from such 
triplet would actually explain the relative weakness the two outer com- 
ponents the quintet reported for nembutal, Fig. 9G. 


Rutonal, Luminal, Mebaral, and Alphenal 

well-defined triplet triplets varying intensity was observed from 
rutonal, luminal, and mebaral, Figs. 10D, and 10K respectively, imme- 
diately after irradiation boric acid glass. The total width the signal 
was the average 140 oersteds. The distance between the components the 
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Fic. resonance absorption signals from irradiated ipral, from 
irradiated sigmodal, from irradiated nembutal, and from irradiated seconal, 
First derivative curves. The above curves were obtained from boric acid glass except 
photograph which was from paraffin. Times after irradiation: hour; 
months; months. Markers are oersteds apart for and oersteds for the 
others. Vertical arrows indicate the position D.P.P.H. resonance. 


large triplet was oersteds. The lines the small triplets were uniformly 
spaced oersteds. The signals were effectively unaltered after months, 
Figs. 10B, 10E, and 10L, and not significantly reduced after months, 
Figs. 10C, 10F, and 10M for rutonal, luminal, and mebaral respectively. 
Irradiation luminal sodium boric acid glass gave the same type 
signals right after irradiation, Fig. 10H; this signal did not change after 
months, Fig. and months, Fig. 10J. Luminal and mebaral were also 
prepared paraffin and irradiated, Figs. 10G and 10N. The same signal 
above with effectively the same component spacings was obtained. 
Assuming that the observed signals are not complex and result from given 
free radical species, the assignment triplet triplets varying intensity 
any given modified area the molecules rutonal, luminal, and mebaral 
not all evident. However, will observed that all the compounds 
reported this study, only rutonal, luminal, mebaral, and alphenal have 
phenyl group one the substituents position Alphenal prepared 
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Fic. resonance absorption signals from irradiated rutonal, from 
irradiated luminal, from irradiated luminal sodium, and from irradiated 
mebaral, First derivative curves. The above curves were obtained from boric 
acid glass except photographs and which were from paraffin. Times after irradiation: 
steds apart for and oersteds for the others. Vertical arrows indicate the position 
D.P.P.H. resonance. 


boric acid glass yielded also right after irradiation, Fig. 11N, signal which 
may taken partially resolved triplet triplets. The total width 
the signal the same that for the other three, 140 oersteds. The signal 
was not changed after months, Fig. 110. therefore very tempting 
relate the common signal observed from these four irradiated molecules 
their common structural feature, phenyl group position The proposed 
explanation would that radiation removes atom from the para 
position the phenyl ring. The ‘‘resonance condition the ring’’ assumed 
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modified the point that the remaining four atoms effectively 
interact with the free electron two groups two. 

If, the other hand, the observed signals are assumed complex, 
they may taken result from the superposition oersted 
wide odd line spectrum the center equal intensity, evenly spaced 
triplet triplets 140 oersteds, total width. The odd line spectrum would 
originate from side chain alterations, seen before; while the equal intensity 
triplet triplets would due the two atoms the malonylurea ring. 
has already been mentioned that the malonylurea ring not planar. There- 
fore, the odd electron not free move around the ring, position acting 
potential barrier. However, this odd electron may quite well encompass 
the remaining five positions the ring. hyperfine structure would then 
expected from the atoms position and 

evident that more experimental work required test the above 
two possibilities. 


Evipal, Delvinal, Phandorn, Nostal, Pernoston, and Barbituric Acid 

Signals varying complexities were obtained from this last group 
compounds. Interpretation most these signals will not possible. How- 
ever, believed that survey clinically important barbiturates would 
not complete without including the above compounds. 

The signal recorded from evipal right after irradiation boric acid glass, 
Fig. 11A, differed from that observed when prepared paraffin, Fig. 11B. 
Both were complex signals about the same total width, 130 oersteds. After 
months boric acid glass, the signal became asymmetric, Fig. 11C, and 
remained unchanged for the following year, Fig. 11D. However, its intensity 
was much reduced. 

Delvinal was prepared boric acid glass only. The signal, right after 
irradiation, Fig. 11E, was asymmetric and had total spread 100 oersteds. 
significant change shape intensity was observed after months, 
Fig. 11F. 

Phandorn, irradiated boric acid glass, gave symmetrical signal 
total width 120 oersteds, right after irradiation, Fig. 11G. After months, 
Fig. 11H, the signal was just intense and practically unchanged. 

The same complex signal, 100 oersteds wide, was obtained right after 
irradiation from nostal prepared boric acid glass, Fig. 111, and paraffin, 
Fig. 11J. This complex signal was reduced evenly spaced quintet after 
months boric acid glass, Fig. 11K. This quintet had line separation 
about oersteds. believed that the isopropyl substituent position 
the origin the quintet. The mechanism would similar that 
seen for ipral and alurate, that say, radiation would remove 
from the group leaving set four nearly equivalent protons. The 
component the signal which had disappeared after months boric 
acid glass, probably due the other substituent position the molecule 
nostal, 2-bromoallyl group. However, the two stable isotopes bromine, 
Br-79 and Br-81, have spin 3/2. therefore difficult speculate 
the possible contribution from this group. 
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Fic. 11. Paramagnetic resonance absorption signals from irradiated evipal, 
from irradiated delvinal, from irradiated phandorn, from irradiated nostal, 
from irradiated pernoston, from irradiated alphenal, and from irradiated 
barbituric acid, First derivative curves. The above curves were obtained from boric acid 
glass except photographs which were from paraffin. Times after irradiation: 
oersteds apart for and oersteds for the others. Vertical arrows indicate the position 


D.P.P.H. resonance. 


The signal obtained from pernoston irradiated boric acid glass was 
comparatively weak right after irradiation, Fig. 11L, and remained unchanged 
after months, Fig. The total width the signal was about 100 


oersteds. 
Barbituric acid was expected yield simple E.P.R. signal which could 
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have been used interpret those from some its derivatives. However, when 
prepared boric acid glass, asymmetric signal, Fig. 11P, oersteds 
wide, was obtained. 


GENERAL DISCUSSION 


The numerous experimental results presented this paper have undoubtedly 
shown that the boric acid glass method preparation may used stabilize 
room temperature the free radicals produced radiation barbituric acid 
derivatives. two cases, with alurate and cyclopal, was possible compare 
the results those obtained directly from the dry salt irradiated and analyzed 
under nitrogen atmosphere. major discrepancies the size and shape 
the signals were reported. Further, many the results were also compared 
those obtained when the compounds were dispersed paraffin and irra- 
diated. Neglecting the instances where hydrogen interaction with the surround- 
ing matrix was more less evident, the similarity the results was easily 
established. The boric acid glass method preparation will eventually 
applied other groups chemical compounds. hoped will useful 
certain cases study room temperature the effect radiation chemical 
and biochemical material. 

The barbituric acid derivatives studied are all the same basic formula. 
The main differences reside the types substituents position the 
malonylurea ring. is, therefore, interesting try correlate some the 
observed paramagnetic resonance patterns with the occurrence specific 
groups the side chains the molecules. quintet average spacing 
oersteds was obtained from molecules which contained isopropyl group, 
nostal, ipral, and alurate. Radiation was assumed remove single 
leaving set four nearly equivalent atoms. quintet with the same spacing 
was also observed from molecules which had methyl butyl group, sigmodal, 
nembutal, and seconal. Breakage radiation the bond between carbon 
and results the same free radical the above case and would therefore 
explain the similarities the results. Barbital and medomin which possessed 
ethyl group gave triplet average spacing oersteds. This triplet was 
taken originate from the two equivalent atoms obtained with the removal 
CH; radiation. also pertinent recall that triplet triplets 
total spread 140 oersteds and average spacings and oersteds was 
recorded only from all the derivatives studied which included phenyl 
group such luminal, rutonal, alphenal, and mebaral. One the tentative 
explanations given, the removal the atom the para position and the 
coupling pair for the remaining four atoms, will eventually have 
checked further experimental work. The remaining derivatives yielded 
signals which were either too complex not interesting from the correlation 
point view. This group included neonal, amytal, and dial which showed 
evidence the presence the free electron close one the atoms 
the malonylurea ring. 

The effect radiation the barbituric acid derivative molecules seem 
centered either one the side chains position The only exceptions 
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are neonal, amytal, and dial which show most important process the 
removal single atom from position the barbituric acid ring. 
reduction the length either one the side chains position might 
result decrease the duration the biological action these irradiated 
barbiturates. However, realized that the electron paramagnetic resonance 
method, combined with the boric acid glass method preparation, freezes 
initial processes only. Many secondary processes are importance under- 
stand radiation damage mechanism. Experiments elucidate the role 
secondary reactions were designed for nembutal. The results will published 
later date. 
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THE HALL EFFECT AND RESISTIVITY TELLURIUM! 


ABSTRACT 


The Hall effect and resistivity seven crystalline samples highly purified 
tellurium were investigated over the temperature range +350° 
The samples were grown from zone-refined tellurium with adaptation the 
method Tammann. The effective impurity concentrations were the range 
atomic per cent. The measurements were made with a-c. 
method high sensitivity. 

X-Ray investigations indicated that the purest sample was single crystal, 
but that the others were polycrystalline. However, differences the observed 
properties result crystal structure were not large. 

high temperatures pure tellurium exhibits anomalous behavior its 
Hall effect, and three proposed explanations this anomaly were investigated 
view the data this research. much lower temperatures, was concluded 
that tellurium behaved simple semiconductor with unique conduction and 
valence bands with non-degenerate carrier distribution. 


INTRODUCTION 


During the past twenty years the properties semiconductors have been 
studied many investigators. Most this effort has been concentrated 
those materials for which practical applications have been found, especially 
germanium, silicon and, more recently, various intermetallic compounds. The 
element tellurium, the other hand, has received very little attention 
although several anomalies its behavior are interest. 

All samples tellurium which have been studied far have been found 
p-type low temperatures regardless the nature the impurity 
content. the case with many materials which are p-type low tempera- 
tures, the effect tellurium changes sign the temperature raised. 
The temperature which this reversal occurs higher for samples with 
greater impurity content and reversal occurs enough impurities are 
present. For samples sufficient purity exhibit this reversal, second 
reversal the sign the Hall effect occurs higher temperature. This 
high temperature reversal unique property tellurium. Several theories 
have been suggested explain this anomalous behavior. The experimental 
evidence found the literature insufficient determine which, any, 
these theories correct. 

The purpose the present investigation provide more accurate 
experimental results range samples, including some with the highest 
degree purity, order throw more light the theories the anomalous 
Hall effect reversal. There disagreement also among the results various 
observers regarding the value the intrinsic activation energy, the variation 
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carrier mobility with temperature, and several other properties tellurium. 
further purpose this work improvement the accuracy measure- 
ment these properties. 

Crystalline samples high purity have been produced using the method 
zone refining which has not previously been used for tellurium. Measure- 
ments Hall effect and resistivity have been made a-c. method 
high sensitivity. This method free from errors due magneto- thermo- 
electric effects which are difficult eliminate from d-c. measurements. 
Although the results obtained improve the accuracy with which the properties 
tellurium are known, has not been found possible confirm without 
ambiguity any the theories the anomalous Hall effect reversal. 


PREVIOUS WORK 


The atoms tellurium crystals are arranged spiral chains 
close-packed structure. The work Scanlon and Lark-Horovitz (1947) and 
Bridgman (1938) indicates that tellurium exists only one crystalline form 
over the range temperatures and pressures normally available. 

The electrical properties tellurium have been studied using evaporated 
films, single crystals, and polycrystalline samples. Bottom (1949) studied the 
anisotropy the resistivity tellurium and found the resistivity 1.9 
times greater perpendicular the c-axis (chain-axis) than parallel it. 
found dependence the Hall effect orientation. 

All samples tellurium far studied have been p-type low tempera- 
tures. Fukuroi al. (1952) found that small additions antimony increased 
the acceptor activity and Moss (1952) showed that tellurium with more than 
0.3% antimony shows semimetallic properties. Kronmuller a/. (1956) studied 
the effects other impurities and found none which produced n-type behavior. 
This not surprising when the position tellurium the periodic table 
considered. 

Nussbaum (1953, 1954) investigated alloys tellurium and selenium. The 
electrical properties the alloys are similar those pure tellurium since 
the two elements and their alloys are isomorphous. 

Johnson (1948) stated that both the electron and hole mobilities vary 
However, later work Nussbaum (1953) contradicted results 
and showed proportionality hole mobility Calculations based 
Bottom’s data agree approximately with this. These results are compared 
with values obtained this work later section. 

Tellurium, common with most semiconductors which are p-type low 
temperatures, shows reversal sign Hall coefficient the temperature 
increased. still higher temperatures second reversal occurs tellurium 
which was first observed Wold (1916). Various attempts explain this 
anomalous reversal have been made Fukuroi and Tanuma (1952), Callen 
(1954), Epstein al. (1957), and Fritzsche (1952). Caldwell and Fan (1959), 
using the results absorption experiments conjunction with the calcula- 
tions Reitz (1957), give support the theory Epstein. These theories 
are discussed detail the section high temperature data. 
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The energy gap between the valence and conduction bands has been calcu- 
lated several observers from their experimental results. Values determined 
from the intrinsic resistivity data are the range 0.32 0.34 (Bottom 
1952; Nussbaum 1953; Kronmuller al. 1956). Fukuroi (1949) and Tanuma 
(1954) derived the energy gap using other data and found value 0.33 ev. 
Earlier workers had found higher values from 0.36 0.38 (Scanlon and 
Lark-Horovitz 1947; Johnson 1948). Measurements the energy gap have 
also been made photoelectric and optical absorption methods Moss 
(1952) and Loferski (1952, 1954). Their values are, general, still higher 
and range from 0.35 0.41 ev. 


EXPERIMENTAL PROCEDURE 


The Hall effect and resistivity crystalline tellurium were measured 
a-c. method which alternating current the specimen and constant 
magnetic field produced the Hall voltage. Several sources error which 
depend thermal gradients the specimen can superimpose spurious e.m.f.’s 
the Hall voltage. These are eliminated a-c. methods measuring the 
Hall effect. 

The samples were cylindrical form with probes for the electrical measure- 
ments welded place shown Fig. The Hall coefficient and resistivity 
such specimens can calculated from the following formulae: 


(1) 
(2) 


these equations, and are the diameter the rod, the magnetic 
induction, the Hall voltage, and the specimen current respectively; the 
distance between the resistance probes and the resistance between them. 
The derivation these equations assumes that the current density uniform 
all points the sample between the resistance probes and that the Hall 
probes are diameter normal the magnetic field. 

The voltages and where the voltage between the resistance 
probes, were measured with differential amplifier, and was measured with 
a-c. Kelvin bridge circuit. The current was found with the relation 

The temperature the specimen was controlled with cryostat low 
temperatures and with furnace high temperatures. The cryostat and 
furnace were designed fit between the pole pieces electromagnet. The 
temperatures were measured with thermocouple and potentiometer. 

The tellurium* before purification had purity about 99.99%. was 
then purified the method zone refining. Crystalline samples were pre- 
pared method based that Tammann (Buckley 1951). The mold 
for crystal growth was pyrex tube with closed capillary one end. This 
was inserted with the open end downward into long vertical pyrex tube 
containing some tellurium. tubular furnace surrounded the system. 


*The tellurium used was provided through the courtesy the International Nickel Co. 
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The system was evacuated and the tellurium was outgassed about 350° 
for hours. The furnace temperature was then increased melt the tellurium 
and inert gas was admitted into the system force the melt into the 
mold. The system was sealed and about minutes was allowed for the 
re-establishment thermal equilibrium. The temperature the furnace was 
then slowly reduced form sample. 

ensure that the cooling advanced along the mold, the furnace was 
constructed give large temperature gradient along its length. For this 
reason, the spacing the heater winding decreased and the amount insu- 
lation increased with distance from the top the furnace. The furnace had 
thermal gradient about C/cm over length about with the 
temperature increasing toward the bottom. The top the capillary was set 
near the top this range and the temperature this point was set slightly 
above the melting point before the inert gas was admitted. 

The furnace was slowly cooled reducing the heater current slightly 
frequent intervals. The large thermal mass the furnace prevented rapid 
changes the temperature for small changes the current. this way 
was possible keep the rate cooling below about Thus the 
velocity crystal growth could kept below about 0.4 cm/min. 

After the sample was removed from the furnace, was immersed hydro- 
fluoric acid remove the surrounding glass. 

jig was devised permit the accurate placement and the support 
the Hall and resistance probes. consists length teflon rod with 
axial hole accommodate the specimen and radial holes (all the same 
plane) accommodate the probes. One side the teflon cylinder was planed 
make flat and provide better thermal contact with the copper plate 
supporting the specimen chamber. 

The teflon jig containing the specimen was strapped copper strip, Fig. 
The Hall and resistance probes consisting short lengths 20-mil platinum 
wire with copper leads attached spot welding were mounted the jig and 
welded the specimen condenser discharge. Some trouble was experienced 
high temperatures due swelling the teflon, resulting breaking the 
welded contacts. 

The copper strip which the specimen was attached was joined long 
monel tube. This tube passes through short glass tube centered the top 
the removable joint the specimen chamber (Fig. 2). short length 
monel below the copper strip and teflon disks the monel tube ensure that 
the specimen centered the chamber. The monel with low thermal con- 
ductivity conducts little heat out into the chamber. The copper plate, 
however, helps give uniform temperature the specimen. 

The thermocouple leads were continuous from the junction the back 
the copper support external terminals. The specimen leads and the copper 
thermocouple lead were insulated with spun glass which was not impaired 
the highest temperatures reached (about 350° C). 

All leads passed through small side tube the removable joint (Fig. 2). 
Rubber tubing containing stopcock grease was placed over this side tube 
and clamped form vacuum seal. Vacuum wax sealed the monel tube. 
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Sealed with wax 


Rubber Clamp 


Specimen and 
Leads 


Joint 


Vacuum Pumps 


Specimen Chamber 
(Pyrex tube with heater 
wound outside) 


Monel Rod 


Teflon Spacers 


Copper Stri 
Specimen 


strapped copper 


Thermocouple junction 


Teflon Spacers 


Fic. The specimen holder. This holder inserted the cryostat furnace with the 
specimen coming between the poles the magnet. 
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The Hall probes were centered the magnetic field visual alignment 
the top the monel rod. The specimen was then rotated give maximum 
Hall voltage and thereby ensure that the Hall probes were normal the 
magnetic field. 

The simplified circuit diagram Fig. indicates the method measure- 
ment. The sensitivity the differential amplifier was approximately 0.04 
microvolt. This permitted accuracy measuring the Hall coefficient 
cm*/coulomb where the current through the sample milli- 
amperes. temperatures near the lower reversal temperature, the resistivity 
large and must less than prevent ohmic heating. high 
temperatures where the resistivity small larger currents were used. Thus 
the small values the Hall coefficient near the lower reversal point could 
measured +0.03 cm*/coulomb while near the upper reversal point the 
accuracy was somewhat better. For larger values the accuracy limited 
the calibration the attenuator the amplifier slightly better than 
1%. The accuracy resistivity measurements being made null method 
was between +0.1% for higher values and for lower values. 


EXPERIMENTAL RESULTS 


Measurements were made over the temperature range —190° +350° 
The samples were numbered serially but complete sets data were obtained 
over this temperature range for only four samples and 10. X-Ray 
investigations with Laue back-reflection camera were made five samples: 
and 10. Laue photographs taken intervals along the lengths 
the samples show that only sample single crystal. The c-axis this 
sample inclined about 60° the rod axis and nearly normal the 
diameter joining the Hall contacts. The other samples are polycrystals with 
various degrees preferred orientation and various grain sizes. Qualitative 
estimates show that the grain sizes and degrees orientation are largest for 
sample and smaller for samples and that order. Samples and 
have very small grains and almost preferred orientation. 

The experimental values the Hall coefficient and the resistivity 
were plotted the usual form: log vs. and log vs. The 
curves are given Figs. 3-7 inclusive. The calculations required the 
analysis these data were made with the values and taken from the 
experimental curves. Rewelding the electrical contacts produced changes 
the results which were small compared with differences between samples. 
Therefore only the average values for each sample are given. 

Some the values determined from the results are summarized Table 
The approximate effective impurity content was determined methods 
discussed the following section. The temperatures and are 
respectively the low temperature reversal Hall coefficient, the high tem- 
perature reversal, and the temperature corresponding the maximum the 
absolute value Hall coefficient occurring between and 
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TABLE 
Characteristic temperatures 


Effective 
impurity 
concentration 
7.4 253 
274 
270 498 
9.0 258 513 288 
1.4 213.5 513 238 
3.7 235 513 267 


(°K) 
Fic. Low temperature measurements 


DENSITIES AND ACTIVATION ENERGIES ACCEPTOR LEVELS 


the acceptors approach saturation temperatures below that which 
intrinsic conduction becomes appreciable, approaches constant value 
from which the number acceptors can found the equation 


(3) 
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3.0 3.5 4.0 4.5 5.0 
(°K) 


Fic. Intermediate temperature measurements 


For sample levels off 5400 cm*/coulomb giving acceptors 
per All the other samples have more impurities and the acceptors not 
reach saturation before intrinsic activation begins. For these samples only 
rough estimates the number acceptor centers are possible this method. 

method suggested Hutson (1957) but not previously used for analysis 
data tellurium can applied samples which not reach saturation. 
and are the energies the Fermi level, the top the valence 
band, and the acceptor levels respectively and the effective mass 
holes then the number holes approximately given 
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2 


Resistivity 


(°K) 


Fic. High temperature measurements Ru. 


and the number ionized acceptors 
(5) 


the extrinsic region where the number electrons the conduction band 
negligible From these equations function can defined. 


2 3/2 

The value can found from the measured Hall coefficient using 
and the function can calculated for number trial 
values the density acceptor levels. For the correct choice one 
would expect log nearly linear function 1/7 since the ratio 
the effective hole mass the electron mass does not vary strongly with 
temperature. 

set such curves for sample was plotted for various values The 
closest approach linearity was for The method very 
sensitive small changes the choice 

None the curves referred the previous paragraph remain linear 
the low temperature end. explanation this behavior that 
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(°K) 
Fic. Low temperature measurements resistivity. 


the acceptors may not all one kind have the same activation energy. 
Equation (6) can readily generalized cover the case more than one 
activation energy. However, since the approximation equation (4) valid 
only for spherical symmetry (scalar effective hole mass, and does not 
apply exactly tellurium, results obtained assuming more than one 
activation energy will not quoted. 

For estimate the acceptor density the remaining samples, third 


method was used. 
The relation for the Hall coefficient assuming isotropic, two-band model 


where the density free electrons the conduction band and the 
ratio electron mobility/hole mobility The condition that vanishes 
then: 
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Sample 


Resistivity (ohm cm) 


10?/T (°K) 
Fic. High temperature measurements resistivity. 
Also any temperature, 
(9) 
and 
(10) 
where 


and depends only the effective carrier masses and considered 


constant. 


Eliminating m2, and from equations (8), (9), (10), and (11) 


(12) 


§ 
1 ot x 
/ 
for 
l/g 
LY 
1 / off 
/ 
i y 
y 
d 
A 
t 
cy 
fof 
y J 
©) 
2. 4 


546 CANADIAN JOURNAL PHYSICS. VOL. 39, 1961 


Thus graph log vs. should linear, and inserting 
the known value for sample and the known energy gap, the values 
for the other samples were found from their reversal temperatures. 
These values are listed Table where the assumption made that 
the total acceptor density. 


RESISTIVITY HIGH TEMPERATURES 


From Fig. can seen that the high temperature intrinsic region 
the resistivity can expressed equation the form 


However, 
(13) 1/p 
Substituting for from (11) obtain 


For many semiconductors and constant, and under these 
conditions the value can determined from the constant The 
values determined this way for samples and are 0.310, 
0.312, 0.294, and 0.294 respectively, giving mean 0.30 ev. Previous 
observers have given values from 0.32 0.34 derived from resistivity curves. 
The discrepancy possibly due the fact, not previously reported, that 
there change slope the resistivity curves for all samples below 380° 
and our values are taken from the straight segment above 380° rather than 
from the mean slope over the intrinsic range. However, will shown, 
not constant for tellurium, nor does that there theoretical 
the above method must doubtful validity. 

The values for samples and were 7.1, 9.0, 8.93, and 10.3 
crystal orientation the various samples and within the range 1.9 
given Bottom (1949) for conductivity perpendicular and parallel the 

The high temperature Hall coefficient, the other hand, does not appear 
depend crystal orientation, being approximately the same for all samples. 
This again agreement with the work Bottom. 


VARIATION MOBILITIES WITH TEMPERATURES 


direct calculation from the experimental results possible only for the 
hole mobility the low temperature range where there are electrons 
the conduction band and the carriers are holes resulting from the ionization 
acceptor impurities. this temperature range 


(14) 
and 
(15) 
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where the number ionized acceptors and the hole mobility. From 
these equations one can obtain the hole mobility eliminating give 


The values log for sample found this way are plotted against 
equation having the form this curve the exponent varies 
from 0.8 for the lower temperature segment the curve 0.5 for the upper 


6 7 
log 


Fic. electron and hole mobility with temperature. 


temperature segment. These results not agree with the theoretical variation 
mobility with temperature for either lattice scattering scattering 
ionized impurities They are, however, agreement with the exper- 
mental results Nussbaum (1953). 

According Debye and Conwell (1954) the theoretical dependence 
mobility for the lattice scattering mechanism was derived Bardeen and 
Shockley the assumption that the conduction band edge located the 
center the Brillouin zone and that the dependence temperature con- 
siderably altered this not so. 

The estimation high temperatures less direct since necessary 
adjust the value one constant obtain fit between the curves for 
high and low temperatures. 

sufficiently high temperatures the impurity density negligible com- 
pared with the density intrinsic carriers. That is, m2, and 


(17) 
(18) 1/p 
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Eliminating from (17) and (18) gives 
(19) 


has been previously shown from the linearity the resistivity graphs 
that Eliminating from these equations gives 


calculated from this equation are shown for value chosen make the 
curve nearly possible continuous with the low temperature segment. 

The values corresponding the estimates were then calculated 
from (19) for the intrinsic range. the intermediate temperature range the 
formula 


(20) 


may used, taking the value given for sample Table 
This formula obtained from equations (7) and (9), together with 
1/p assuming that i.e. that intrinsic carriers still 
predominate. The correcting term less than 20% for any value 
plotted. 

The electron mobility can now obtained from the definition i.e. 
Estimates the electron mobility are also plotted Fig. 


THEORIES THE HALL EFFECT REVERSAL 


Three theories have been proposed for the anomalous reversal sign 
high temperatures. discussion these theories follows. 


(a) The Defect-Acceptor Theory 

this theory proposed Fritzsche (1952) and discussed Tanuma 
(1954) the reversion p-type behavior high temperatures supposed 
result from additional holes produced thermal-lattice defects acting 
electron traps. The authors this theory assume that the density 
thermal-lattice defects given 


where the density lattice sites, depends the ratio, and 
close unity, and the defect-activation energy. Equation (22) applies 
defects the Schottky type predominate. The authors also assume that 
all defect-acceptors are activated that high temperatures. 
produce predominance holes high temperatures, must exceed 
that, the temperature rises, begins appreciable somewhat 
below the higher reversal temperature With further rise temperature, 
quickly dominates the activity and causes positive above 
Thus, high temperatures somewhat above T2, should the dominant 
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source carriers and the graph log vs. should show slope propor- 
tional Thus the curve logpvs. high temperatures should 
exhibit slope which increases the rising temperature from the value corre- 
slope was not observed, consider this theory untenable. This does 
not imply that defects not occur but merely that they are not sufficiently 
numerous account for the reversal sign Hall effect this 
temperature. 


(b) The Dual Valence Band Theory (Epstein al. 1957) 

According this theory, lower valence band high mobility accommo- 
dates more holes the temperature rises that the average hole mobility 
increases relative the electron mobility. then becomes positive when 
the mobility ratio becomes less than unity (cf. equation (17) for 
the high temperature Hall coefficient). 


(c) The Dual Conduction Band Theory (Fukuroi and Tanuma 1952; Nussbaum 
1953; and Callen 1954) 

This resembles theory that the high temperature reversal 
due decrease with rising temperature. The proposed mechanism is, 
however, second conduction band low mobility lying somewhat above 
the edge (bottom) the main conduction band. The higher band accommo- 
dates more electrons the temperature rises. 

either theory (c) the mobility ratio varies with temperature 
such way that becomes zero the reversal temperature, agree- 
ment with the curves shown Fig. The curves for and Fig. 
suggest that the electron mobility varies with temperature nearly the 
normal way whereas the hole mobility does not. This would seem indicate 
that the dual valence band theory more likely explanation the 
anomaly than the dual conduction band theory (c). 

The work Caldwell and Fan (1959) also supports this view. Their 
measurements infrared absorption indicate the existence holes 
lower valence band and also that these holes have lower effective mass 


and higher mobility. 
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POWER SERIES EXPANSION THE MASTER 


vAN KAMPEN 


ABSTRACT 


order solve the master equation systematic approximation method, 
expansion powers some parameter needed. The appropriate parameter 
the reciprocal size the system, defined the ratio intensive and extensive 
variables. The lowest approximation yields the phenomenological law for the 
approach The next approximation determines the mean square 
the fluctuations about the phenomenological behavior. equilibrium this 
approximation has the form linear equation. The higher 
approximations describe the effect the non-linearity the fluctuations, 
particular their spectral density. The method applied three examples: 
density fluctuations, Alkemade’s diode, and Rayleigh’s piston. The relation 
the expansion recently given Siegel also discussed. 


THE PROBLEM 


Let fluctuating physical quantity, P(a,t) its probability distribution 
time possesses Markov character, P(a,t) obeys the equation 
where the transition probability per unit time from Equation 
(1) the differential form the equation and 
often called the The right-hand side can formally 
expanded series, 

where the nth moment the transition probability, 


The expansion (2) was given Kramers (1940) and Moyal (1949). The 
equation (Fokker 1913, 1914; Planck 1917; Wang and Uhlen- 
beck 1945) obtained from cutting off after 


Clearly, this cutting off not systematic approximation procedure, because 


(2) not series expansion powers some small parameter.* Our purpose 
provide systematic approximation method for For special case 


received November 28, 1960. 

Contribution from the Institute for Theoretical Physics, Rijksuniversiteit Utrecht, 
Nederland. 

*For instance, one might just well agree rearrange the terms (2) according suc- 
cessive derivatives and omit all terms with higher derivatives than the second. one 
might, like Lax (1960), expand “orders 

step from (1) (2) not entirely harmless. can shown for the example Section 
that certain terms are lost. These terms are not analytic our expansion parameter and 
therefore, treated power series expansion. Fortunately they are also extremely 
small. 
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(viz., the current fluctuations diode, see Section such expansion 
has been given previously (van Kampen 1959); now treat the general 


case. 
Note.—The name equation sometimes confined more 


restricted kind equation, wit 


with constants and shall call this equation—first studied 
Rayleigh (1891)—the linear equation. 


POWER SERIES EXPANSION (2) 


Rewrite the transition probability function the length 
the jump, and its starting point a’, 


(5) W(ala’) 


far the variable has not been specified; may extensive 
intensive quantity. now stipulate that extensive. The corresponding 


where measure for the size the system. The essential point that the 
size the jump, Aa, properly expressed terms the extensive quantity, 
whereas the dependence (5) properly expressed terms the intensive 
quantity More precisely, write 


the function longer involves implicitly, contrast with 
That this can seen from the examples Sections 

Since the dependence has now been made explicit, possible 
expand reciprocal powers The moments when expressed functions 


longer contain 

Next one has estimate how the successive derivatives P(a,t) depend 
For this necessary anticipate the kind solution one interested 
in. one studies equilibrium alone, the interesting values are order 
However, equation (2) can also used describe states that differ macro- 
scopically from equilibrium. such states will have mean value order 
and spread around this mean order Accordingly put 


the function time determined presently, the new variable. 


*This requirement may weakened: sufficient that can expanded power 
series Compare the example Section 
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After transforming from the variable the variable (2) becomes 
n=1 
The leading term the 


This can made cancel the second term the left subjecting 
the condition 


(6) 
This just the macroscopic, phenomenological law. The remaining equation 
(x,t 


n=2 n! 


convenient put that the new time scale increases with 
agreement with the increase the relaxation time. Then, expanding the 
right-hand side (7) powers and rearranging terms 


Or m=2 m! n=1 


this equation all factors are shown explicitly, that the relevance 
the several terms indicated their degree The positive powers 
have been removed means (6). 

For many purposes useful find the successive moments rather than 
the distribution function itself. The general formula (8) gives 


m=2 m! 


the upper limit being either whichever 
noteworthy that for many special cases exact results can extracted 
from (9). For instance, linear function, 


This can solved give const. will usually convenient 
choose equal that our equation then shows that 
for all later Similarly, happens linear, one finds 
rigorous equation for all a,(X) are linear, possible find succes- 
sive equations for all (x*); such case the example Section 


*More directly, from either (1) (2) one has rigorously linear, 
this equal which proves that (a) obeys the macroscopic equation. 
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THE ZEROTH ORDER APPROXIMATION 
Neglecting all terms that vanish for 


This would linear Fokker—Planck equation the type (4) the 
cients were not time dependent. Nevertheless the initial-value problem can 


solved. 


(11) P(x,0) 

define 

(12) s(r) log ’ 

that s(0) and ds/dr Then (10) reduces 


The solution this equation given Chandrasekhar (1943) 


where 


Thus finally find 
—s\2 


where 


order find the probability distribution with given initial dis- 
tribution 
one first has solve 


q 
q 
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This can always done quadrature, 
One then has determine and from (12) and (13). The final result 


This Gaussian peak, whose width the order the equilibrium fluc- 


tuations, and whose center moves according the phenomenological law.* 
Higher approximations will destroy the Gaussian shape. 


FIRST EXAMPLE: DENSITY FLUCTUATIONS 


box volume communicates through small hole with infinitely 
large volume, filled with dilute gas with density (Fig. 1). The number 
particles obeys the master equation (in suitable time units) 


(15) 


the extensive variable, the density the intensive one. The 
transition probability 


The moments are given 


*This result has been employed previously evaluating the entropy non-equilibrium 
states (van Kampen 1959). 
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The macroscopic equation (6) 


with the solution 


was expected. One finds furthermore 


s(r) 
Hence (14) takes the form 


This example sufficiently simple check the result (16) several ways. 
First one finds rigorously from either (1) (2) 


These equations can solved give 
This agrees exactly with the mean and variance the Gaussian (16). The 
higher terms therefore, while destroying its Gaussian shape, 


not affect the first two moments. 
second check one may compare the equilibrium value (16) 


with the correct expression 
N 
For this purpose expand log near its maximum, which 


This agrees with (17) when terms order are neglected. 
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final check possible solve (15) exactly with the aid the char- 
acteristic function, defined 


The solution 
Let this expanded second order 
This approximate expression just the characteristic function the Gaussian 
distribution (16). 
EQUILIBRIUM FLUCTUATIONS 
order describe the fluctuations equilibrium one has substitute 


for the general equation (8) its equilibrium value that the 
quantities become constants.* zeroth order equation (10) then 


reduces 


This the linear equation (4). The solution with initial dis- 
tribution (11) well known, 


(Of course, for all realistic cases <0, The difference between this 
solution and (14) consists the fact that the present one only valid for 
initial distributions which are confined range order the x-scale, 
i.e. order the a-scale, around the equilibrium value. 

The next approximation adds the right-hand side (18) 


(18) 


(Note that the resulting equation can also obtained directly from (2) 
breaking off after the second order after the first order, 
after the zeroth order, and omitting all higher terms.) The first order 
equation can solved the order treating the additional term (20) 
perturbation. That is, one substitutes the zeroth order solution (19) 
(20), and solves the resulting inhomogeneous equation for 

Alternatively one can find the successive moments, using (9), which for 
the region simplifies 


n=1 


*We write nth Aa, taken the equilibrium. denotes the kth deri- 
vative an(X) with respect also the equilibrium value 
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Explicitly, 


Now let the initial distribution given (11), that. 


then possible solve the set equations (22) successive orders and 
such units that —1, Then one has zeroth order 


omy’ 


Subsequently one finds order 
Finally order 


example the moment method compute the spectral density 
the fluctuations order The autocorrelation function obtained from 
(23) multiplying with and averaging over the equilibrium distribution 
For this one needs the equilibrium moments, which can found 
any desired order from (21) putting the left side equal zero 


With the aid these values one finds for the autocorrelation function 
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The last term may combined with the term give with 


The constant term exhibits the fact that does not vanish, that 
Hence (24) may also written 


From this follows the fluctuation spectrum 


42 — hai! | 


Bernard and Callen (1960) arrived the conclusion that the fluctuation 
spectrum entirely independent the non-linear terms. According (25), 
this only correct for (In fact can seen from our examples that 
both coefficients (25) not happen vanish.) 


SECOND EXAMPLE: DIODE 
(Alkemade 1958; van Kampen 1960) 
Consider the circuit Fig. The whole system kept constant tem- 
perature The two electrodes have different work functions, and electrode 


supposed operate under saturation conditions. The number electrons 
the left condenser plate obeys the master equation (in suitable time 


units) 


(26) 


Here and constant depending both work functions. The 
expansion parameter the capacity the condenser, rather 


For the extensive parameter take that eN. The 


The macroscopic equation 


(27) 


| 
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Equilibrium obtains for log which determines the contact potential 
Writing (27) the form 


one recognizes the usual characteristic the diode. 
The solution (27) 


(28) 
where integration constant given 


One subsequently finds 


s(r) = log 1+e” 
and 


(e" + e”? 


VAN KAMPEN: MASTER EQUATION 561 


This agreement with the result obtained previously (van Kampen 1960). 
The distribution P(N,t) the Gaussian (14) with taken from (28) and 
from (29). The equilibrium distribution is, this order, 


The spectral density the equilibrium fluctuations obtained from (25) 
substituting 


a3 = 0. 


The result again identical with the one obtained previously. 

According (30) the average number electrons equilibrium 
According (23), however, there higher order correction 
exact one, because the exact equilibrium distribution 


can easily found from (26).] means that the equilibrium distribution 
found from the zeroth order equation, viz. (30), differs from the exact one 
through higher order terms. this case these higher order terms not 
destroy the Gaussian shape, but merely cause shift. Accordingly the value 
for the contact potential, log not the exact one. Using the explicit 
value for one finds for the exact contact potential 


where and are the work functions both electrodes. The fact that the 
equilibrium value the potential lowest order, does not exactly corre- 
spond with the expected value, has given rise some discussion 
(MacDonald 1957; Alkemade 1958; Marek 1959). 


THIRD EXAMPLE: PISTON 
(Rayleigh 1891; Zernike 1929) 

one dimension heavy particle with mass subjected collisions 
with light gas molecules with mass and velocity distribution f(v) (Fig. 3). 
collision molecule with velocity changes the velocity the heavy 
particle into 


One thus obtains master equation for the velocity distribution with 
transition probability 


N QX, 2 i 
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where the number molecules per unit 

Although the concept intensive and extensive variables somewhat 
alien this system, similar distinction can made. Indeed, the individual 
jumps are due collisions with the gas molecules, and are therefore 
appropriately measured terms the momentum the other hand, 
the probability distribution these jumps depends the velocity itself. 
Thus one led the following identification 


The transition probability (32) then takes the form 
W(a;Aa) 


This expression not entirely independent However, can expanded 
powers that the previous treatment can still applied with 
minor modifications. This leads the usual expansion m/M. 

more elegant application the previous treatment, however, possible 
one puts 


The transition probability now takes the form 


entirely independent Hence the results the previous sections can now 
literally applied. 

The macroscopic equation, which describes the approach its 
equilibrium value rather unmanageable. therefore only study 
the equilibrium fluctuations. The zeroth approximation the linear 
Planck equation (18) with 
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the Maxwell distribution, 


a= RT Vv m e 


After transforming back the usual variables the Fokker- 
Planck equation becomes 


This the equation, but for the appearance 
rather than m/M. The lowest order the expansion m/M can, course, 
easily obtained from (35) simply with 

The next term the expansion (20), but all coefficients vanish because 
symmetry.* The term order can also found; one puts m/kT 
and m/M the result is, when written expansion 


2,8 


The terms the second line are order relative those the first line. 
They agree with the result Siegel (1960). 

advantage expanding rather than that the equilibrium 
distribution, when expressed the normalized co-ordinate independent 
Hence must satisfy each term the expansion separately. This can 
verified from (36). The equilibrium distribution, when expressed 
contains higher terms and hence does not satisfy the separate terms 
the expansion (see next section). That the reason why not 
exact solution (35). 


RELATIONS FOLLOWING FROM THE KNOWN EQUILIBRIUM DISTRIBUTION 


far the equilibrium distribution has been found special solution 
the master equation. often happens, however, that known priori. 
Then this known solution can used, substituting the general 
equation, find relations between the a’s. For instance, follows from (18) 
that lowest order the equilibrium distribution given 


*Recently MacDonald (private communication) has suggested asymmetric modification 
the Rayleigh piston model, for which these coefficients not vanish. 
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2 

The variance this Gaussian usually known from equilibrium 
statistical mechanics. Hence one obtains relation between and which 
corresponds the usual Einstein relation for Brownian movement. 

Additional relations can only obtained higher order terms the 
equilibrium distribution are also known. some cases one knows that (37) 


exact all orders. Thus, taking again —1, 


substituting (21) the separate powers yield 


whichever smaller. easily seen that for these equations are 
identically satisfied. For they lead 


For 


More generally one finds that all even derivatives vanish (at equilibrium) 
whereas for odd 


‘ —1 


The Rayleigh piston furnishes example equilibrium distribution 
that known all orders but not independent Indeed 


Putting m/kT one has instead (88), 


Substituting (21) and separating the various powers one finds the 
same way before that (39) still valid, but (40) must replaced with 


These relations can verified explicit computations. 
MacDonald (1957) has emphasized that (37) need not correct all 
orders. fact, our example Section certainly only correct zeroth 


order. Yet may asserted that the equilibrium distribution must 
Gaussian under the following conditions. 
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(i) The energy the total system quadratic the fluctuating quantity 
and 

(ii) the transition probability contains some physical parameter 
(other than factor, which permits scale down its magnitude without 
affecting its functional dependence and a’. 

The first condition ensures that, according equilibrium statistical 
mechanics, the probability distribution 


For instance, this condition rules out the case density fluctuations. The 
second condition necessary because (41) only exact for infinitely small 
interaction. The scaling parameter permits the limit infinitely small 
interaction, that (41) becomes exact. may then concluded that (41) 
must also exact for all other values the scaling parameter, because 
does not depend the magnitude but only its functional 
dependence and a’. 

Both conditions are fulfilled the case Rayleigh’s piston, the variable 
defined according (33) rather than (34), the scaling parameter being 
Hence the Gaussian distribution valid for all values M/m—which 
course well known. the case the diode, the electrical energy the 
circuit proportional the square that equilibrium 
statistical mechanics indeed leads (31). The exchange electrical energy 
with thermal energy the heat bath can scaled down decreasing the 
area the electrodes. (This area should occur factor the right-hand 
side (26), although for convenience was not written explicitly but absorbed 
the time unit.) This explains why (31) has exact. 


COMPARISON WITH SIEGEL’S EXPANSION 

Recently Siegel (1960) has given ingenious transformation the 

Kramers—Moyal expansion (2), namely 

The expansion only applies the neighborhood equilibrium (our Section 
5). The units are chosen such that and assumed that the 
equilibrium distribution rigorously The coefficients are related 


1 1 (m—k+1) +c 


2 


where denotes Hermite polynomial, 


Siegel suggests that (43) the successive values should constitute 
successive orders approximation, and shows for the case the Rayleigh 


2 
a 
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piston that, indeed, Aim order being given (33)). Aswe have 
obtained power series expansion for all cases, now possible check 
Siegel’s suggestion generally. From (43) follows with the aid (44), 

1 m 


This shows that Ajm order that the successive terms 
Siegel’s expansion are indeed increasing order 

Siegel’s aim rather give systematic expansion, which achieves two 
other purposes: 

(i) each approximation (i.e., (42) broken off arbitrary value 
the correct equilibrium distribution should exact solution; 

(ii) each approximation should meaning that all solutions 

The first purpose seems laudable one, albeit mainly for reasons 
elegance. Admittedly, awkward feature the systematic expansion 
that leads the denominator rather than (35), although 
not incorrect. For practical purposes, however, such the computation 
the fluctuation spectrum, systematic power series expansion more 
convenient. 

Siegel achieves the second purpose mathematical construction, which 
may roughly described follows. After cutting off the expansion the 
master equation mo, suitably chosen higher order terms are added such 
make the resulting equation This construction could also 
applied our expansion (8). However, not feel that one ought 
impose the requirement (ii) being ‘‘semidefinite” the successive approxi- 
mations. true that, for each for which the cutoff series not ‘‘semi- 
there exist solutions that not tend the equilibrium distribution 
Such behavior, however, only shows that the approximation used 
not good enough for these adding higher order terms 
make the equation one does not improve the approximation. 
One replaces the incorrect solutions other solutions, which tend 
but are otherwise equally incorrect. 


Aim = 
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THE FINITE BEAM SPACE-CHARGE LIMITED DIODE 
NOISE TRANSDUCER! 


McFARLANE? 


ABSTRACT 


method presented for calculating the noise properties cylindrical 
electron beam from space-charge limited diode electron gun. The multivelocity 
character the beam considered and correction made for finite beam 
diameter. The ratio anode noise current full shot noise found depend 
where the frequency and the anode voltage. 


INTRODUCTION 


Haus (1955) has demonstrated that four independent parameters are 
required completely specify the noise properties electron beam from 
thermionic cathode. shows that under small signal, single velocity 
assumptions, the electron beam can represented lossless passive trans- 
ducer, where output noise quantities are linear functions input quantities. 
There are found two noise invariants, and which remain unchanged 
under such conservative transformation. The non-conservative nature 
multivelocity electron beam was indicated Siegman (1957), who made 
detailed calculations the four noise parameters, and asa function 
the Fry-Langmuir distance variable the region space-charge 
limited diode immediately beyond the potential minimum. 

The beam noise properties beyond the electron gun anode are more 
directly accessible experiment than are those near the potential minimum. 
The present work describes method for calculating the anode noise para- 
meters terms input quantities plane slightly beyond the potential 
minimum. This essentially determination the transducer matrix ele- 
ments where the formal representation 


Vv, Vin 
®, Pin 
(1) Il, [715] Thin 
Ag Ain 


Input values II, and used the numerical calculations were those 
found Siegman. Included the analysis point-by-point correction for 
the finite beam diameter (Kornelsen 1957). This facilitated the use 
analog computer tosolve modified form the electronic equation (Hutter 
1952, 1953). particular electron gun examined three frequencies and 
over range anode voltages. 
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The results the computation for show dependence frequency and 
gun anode voltage not previously reported. companion paper presents the 
results measurements made this author and other workers and compares 
them with the present calculations. 


THEORY 


include the effect finite beam diameter the plasma frequency 
cylindrical beam, Parzen (1952) proposed that the d-c. current density 
replaced the electronic equation (Hutter 1952, 1953) effective value 


Joe, given 
(2) 


where the usual plasma frequency reduction factor. Under the assumption 
zero total current, the electronic equation becomes 


2 , 2 
(3) oy 1 Via 0 


3 
where the a-c. signal current modulation 


uo 
and the a-c. velocity modulation 


Kornelsen (1957) has assumed the region the gun between the potential 
minimum and the anode Llewellyn gap with zero initial acceleration. 
The d-c. velocity then given 


(6) 
where the transit time from the potential minimum and initial 
velocity used characterize the multivelocity beam. Equation (3) becomes 


where 


(8) 


and this the final form solved the computer. 

The plasma frequency reduction factor taken that the first 
radial mode cylindrical beam with surrounding conducting cylinder 
infinity. applying this form correction accelerated beam, 
assumed that the effective current density any point the gun can 
characterized the reduction factor for drifting beam the same radius 
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and d-c. velocity. Higher order modes and mode coupling are neglected. The 
result here obtained more accurate than the W.K.B. solution Rowe 
(1952) and the more restricted method Kosmahl (1959). 

The development the modified electronic equation (7) involves small 
signal, single velocity assumptions which are unquestionably violated the 
region the potential minimum. might also remarked that neither the 
modified electronic equation nor the Llewellyn—Peterson equations (Pierce 
1950) admit the existence development real kinetic power when the 
input plane selected the potential minimum, although the conservation 
maintained. overcome these difficulties, the present method 
solution restricts the application the electronic equation region 
the gun over which the assumptions made its derivation can reasonably 
presumed valid. The statistical calculations Siegman al. provide criterion 
for determining appropriate input plane and make possible specify 
the initial noise quantities used. 

input plane where the beam has developed real kinetic power, there 
exists correlation between the current and kinetic voltage fluctuations. 
shown Bloom (1955) and Beam (1955) the noise properties such 
cross section can represented three independent noise inputs. 
and are the total mean square fluctuations per unit bandwidth per unit 
beam area for the input kinetic noise voltage and current, the three inputs 
are specified 


with 

The correlation coefficients are such that 

<1. 


For convenience discussion, the following definitions are made, using 
quantities calculated Siegman. Their relation Haus’ noise parameters 
indicated. 
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Values for these quantities are published function the distance para- 


meter 
From (9) and (10) apparent that 


input plane having been specified, the values and the product 
are determined. will demonstrated, proper selection leads the 
conservation real kinetic power and with specified there results 
unambiguous determination all input quantities. 


THE ANALOG COMPUTER SOLUTION 
Kornelsen (1957) has discussed detail the use analog computer 
solve the finite beam electronic equation. Only short summary the method 
will presented here. schematic diagram the computer shown 
Fig. The procedure used Kornelsen briefly follows. 


Recorder Recorder 


Fic. Analog computer schematic. Two operational amplifiers are integrators. 
corresponds the quantity equation (7) and its time derivative. The circuit 
constants are determined equating coefficients the parallel equations. 

(1) The computer set solve equation corresponding exactly 


form equation (7). 
(2) Two time normalizations are made. the independent variable, 


real time the computer, and gun transit time, normalization such 
that point representing the gun anode 


and this serves define and 


n 
a 
i q 
f(t) 
> 
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(3) The quantity for given beam radius and frequency, 
determined function time using equation (6), and assumed initial 
velocity The function indicated Fig. the term 
normalized its maximum value and this set the computer loaded 
tapped potentiometer. 

(4) The necessary circuit constants are determined equating coefficients 
the parallel equations. 

(5) Kornelsen then introduces the initial values the noise current and 
velocity setting appropriate values for and the voltages repre- 
senting the initial values the function and its time derivative. The 
computer output voltages and V2, can then used evaluate the 
anode noise current and velocity. 

Considerable simplification possible step (5) and the results can 
conveniently summarized evaluating what are equivalent the transducer 
matrix elements equation (1) (McFarlane 1959). This follows from the 
observation that, for given gun voltage and frequency, the computer output 
voltages are linearly related the input voltages. 


The coefficients and are positive and their determination consti- 
tutes the complete solution. 

Noting that one the computer integrator time constants and that 
are the usual Fry-Langmuir normalized potential and distance variables 
respectively, define 


(13a) o= 2aR2C2ao— bo, 
(13d) 


gun anode voltage, 


2 é 


transit time from the potential minimum the plane determined 
from equation (6) and assumed value for 

After some lengthy algebra, summarized the Appendix, found that 
the noise quantities the electron gun anode are 


d 


; 


ti- 
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1 a 2 


2 


These have been obtained combining the noise outputs each the 
three independent inputs equation (9). the computation chosen 
conserve II’ equation The quantity represents the ratio the 
total mean square current fluctuations full shot noise. companion paper 
(McFarlane 1961) discusses experimental measurements this quantity 
the anode space-charge limited diode electron gun. 


CALCULATIONS 


The specific diode studied was the one which the indicated measurements 
were made. had the following properties: 


Type Pierce-parallel flow 
spacing 0.156 in. 
Effective beam diameter 0.023 in. 
Cathode temperature 1160° 


The evaluation the anode noise quantities required the specification 
the plane which the calculation was start, the four fluctuation terms 
(10a, this plane, and the transit time from the potential minimum 
the input plane. For the present work input plane was con- 
venient for computer operation. The potential this plane 0.359 volt above 
the minimum potential the gun and this plane and have attained 
their final values. 

The input noise quantities depend upon the frequency parameter 
Where wpm the plasma frequency the potential minimum 
defined 


2 1/4 
1/2 


Table shows the frequency parameter for the cases studied. 
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TABLE 
Frequency Frequency parameter 
(Mc/s) Gun voltage 
1,400 100 0.845 
200 0.502 
300 0.371 
500 0.253 
4,250 100 2.53 
300 
500 0.756 
1,000 0.450 
9,520 300 2.49 
1,000 1.01 


obtain the noise input terms for frequency parameters other than those 
calculated Siegman, the input plane was fixed and values 
II’, and were found interpolation. Siegman shows that the use 


(16) 


for the evaluation using the Llewellyn—Peterson equations, gives good 
agreement and 0.25 but predicts values for greater than 
those found from the statistical computations for frequency parameters 
0.5 and 1.0. The present calculations were based value equal 
twice that given equation (16). believed that more accurate 
determination was obtained larger values the frequency parameter. 

the present work, the value was obtained requiring invariance 
For all the cases studied, was less than its maximum allowable value 
unity. For frequency parameters greater than 0.5, was also less than 
unity. decreased from 0.5 0.25, however, value 1.035 was 
indicated. The discrepancy was not considered serious for the present work 
but was interpreted indicating that the validity the method calculation 
might open question for values the frequency parameter less than 
0.25. The difficulty can apparently removed for the lower frequencies 
advancing the input plane toward the anode. This possibly implies the 
necessity minimum transit angle between the potential minimum and 
the plane which the computation started. 

The results the calculation are summarized Table II. Both input 
and anode quantities are given. The theoretical dependence the smoothing 
below shot noise the anode noise current fluctuations can seen examin- 
ing the quantity the calculated values for which are shown 
Table 

CONCLUSIONS 


restricting the solution the electronic equation the conservative 
region space-charge limited diode, the noise transducer properties 
finite electron beam have been determined. The general form the equations 


i 


McFARLANE 


o 


apouy 


575 
v 
} © No 
ose 
use 
ual 
ate 
ice a| FBS 
lue 
ork 
ion 
lan 
the 
Ing . . . . . 
in- 
ive 
ons 


576 CANADIAN JOURNAL PHYSICS. VOL. 39, 1961 


TABLE III 
Frequency 
(Mc/s) Gun voltage 

1,400 100 
200 

300 

500 

4,250 100 
300 

500 

1,000 0.946 

9,520 300 
1,000 


shows similarity what might expected for lossy non-reciprocal 
transducer (Haus 1955). This due the presence unspecified corre- 
lation coefficient requiring that the transformation conserve real kinetic 
power, found that the parameter also conserved lossless trans- 
formation, remains within per cent its initial value. The dependence 
anode voltage and frequency variance with earlier theories. The 
present calculations indicate that, for the particular gun least, pro- 
portional the square root the frequency and inversely proportional 
the gun voltage. Theories not including the correction for finite beam diameter 
and starting the potential minimum show proportional The 
difference these two types variation sufficient that experiment should 
able distinguish between them. companion paper (McFarlane 1961) 
undertakes this investigation. 
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APPENDIX 
THE CALCULATION ANODE NOISE QUANTITIES 


The general case examined where the three independent noise inputs 
equation (9) are specified any plane, and the solution carried out 
determine the noise values the anode. conversion from velocity modulation 
kinetic voltage has been made, where the kinetic voltage 


MU Wy 


(A.1) 


The problem solved formally using the square roots quantities 
representing quadratic content. quadratic content that determined 
the gun anode, the results are physically meaningful. Unit bandwidth 
assumed. The algebra very long (McFarlane 1959) and only the results 
will presented. 


(1) Kinetic Voltage Input Only 
found that the anode quantities are 


kT, 
2, 1/2 c vila 
— p2\1/2 c 2a 
where 
aRsCow 
* — 
d = 2aRoCo. 


The computer input condition 


(2) Current Input Only 
found that the anode quantities are 


2, 1/2 1/2 Vo V; 
The computer input condition 


(3) Correlated Input 
The computer input for the correlated noise input 


| 

| 
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The calculation anode noise requires the determination real and imaginary 
components. 


wa 26 {| Yan | [ — 


where and refer the use the real and imaginary parts 

Using the appropriate computer input terms for each the three inde- 
pendent noise inputs, evaluating the required computer outputs equations 
(12a), and noting the relations equations and and the 


— 
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MICROWAVE NOISE ACCELERATED ELECTRON 


ABSTRACT 


Measurements have been made the noise current fluctuations the 
electron beam from space-charge limited diode electron gun, 1,400 Mc/s, 
4,250 Mc/s, and 9,520 Mc/s. Theories which not consider the finite beam 
diameter and the multivelocity nature the stream the region the potential 
minimum are poor agreement with experiment. The measurements here 
reported and those other workers are compared with the results theoretical 
calculations which these two effects are considered. Significant improvement 
agreement with experiment realized. 


INTRODUCTION 


Since the early measurements Cutler and Quate (1950) the noise 
current fluctuations beam electrons from space-charge limited diode 
gun, workers have been interested the noise figure longitudinal beam 
amplifiers, using such guns. Haus (1955) showed that two parameters, and 
II, are significant determining the minimum noise figure these devices. 
Subsequent investigations have centered the measurement these two 
quantities and the construction guns with increasingly complicated geo- 
metry and axial potential profiles, with view decreasing the noise figure 
the devices. Most the work has been done the region 3000 Mc/s. 
Although four independent parameters are required specify the noise 
properties electron beam (Haus 1955), experimental study these 
function frequency and gun voltage has received limited attention 
(Kornelsen 1957; Vessot 1957). Where agreement has been obtained between 
experiment and the more elementary theory, appears have been the 
consequence the proper selection frequency and gun operating conditions. 

The present work reports measurements space-charge limited diode 
gun the anode noise current fluctuations 1,400, 4,250, and 9,520 Mc/s, 
over range gun voltages from 100 1,000 volts. relating these fluctua- 
tions full shot noise, what has been obtained essentially measure 
Haus’ parameter i.e. the self-power density spectrum the noise current 
modulation. The smoothing below full shot noise calculated from theories 
which not take into account finite beam diameter and the multivelocity 
nature the electron beam the region immediately beyond the potential 
minimum are poor agreement with the experimental measurements. Con- 
sideration these effects has been given companion paper and the 
agreement between theory and experiment significantly improved. 
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THE EXPERIMENTAL VACUUM TUBES 


The electron gun used for the reported measurements was parallel flow 
Pierce type designed produce cylindrical beam electrons. employed 
oxide cathode front which wasa beam-forming electrode the required 
Pierce shape. 0.031-in. aperture the beam-forming electrode determined 
the diameter the electron beam leaving the cathode. The 
spacing was 0.156 in. which included 0.004-in. retraction the cathode 
behind the beam-forming electrode provide thermal isolation. simplify 
construction, the anode was approximated plane. The tubes used 
1,400 Mc/s and 4,250 Mc/s had anode aperture and inter- 
ception beam current the anode was small. improve coupling the 
electron beam the 9,520 Mc/s cavity, the anode aperture was reduced 
0.031 in. and this resulted increase interception for the two highest 
frequency tubes. 

The resonant cavities used couple the noise power resulting from beam 
current fluctuations the microwave receiver were located just beyond the 
electron gun anode. The tubes for measurements 1,400 Mc/s and 4,250 
Mc/s employed disk seal structure which external cavities were attached. 
9,520 Mc/s, the small cavity size required that placed entirely within 
the vacuum envelope. coaxial line coupled the noise power the receiver. 

The tube properties are summarized Table observed that for 


TABLE 
Electrical and mechanical properties the tubes studied 
Tube number D-1 D-2 I-1 
Cavity type External External Internal Internal 
Anode and cavity aperture (in.) 
Frequency (Mc/s) 4,030 4,250 9,520 9,520 
1,395 1,455 
Interception 3.8%-1.5% 0.2%-0.02% 2.8% 


200-500 200-500 


each type, the tube with the higher perveance value had the lower inter- 
ception. This interpreted mean that the cathodes tubes D-1 and I-2 
were positioned further behind the beam-forming electrode than for optimum 
operation. All tubes employed long cylindrical electron collectors biased 
volts above the anode potential order reduce the effect secondary 
electrons. Sealed each tube was ionization gauge provide continuous 
indication pressure within the vacuum envelope. During operation, pressures 


MEASURING APPARATUS 


The receiving system was based the radiometer technique Dicke 
(1946). block diagram the apparatus shown Fig. The electron 
gun was gated square wave and the beam noise coupled 
conventional microwave receiver. Following detection bolometer, the 
receiver output was synchronously detected with the modulating wave form 


— 
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ELECTRON RESONANT CRYSTAL LOW NOISE 
GUN CavITY MIXER PREAMP 


ELECTRON LOCAL auolo 
BEAM 
MODULATOR OSCILLATOR AMPLIFIER 


CALIBRATED 


SYNCHRONOUS 
DETECTOR 


SQUARE 
WAVE 
GENERATOR 


ATTENUATOR 


INDICATOR 
Fic. Block diagram noise measuring system. 


reference. operation, the receiving apparatus was linear 0.1 over 
greater than 30-db range input signals. The noise electron beams 
less than was readily detected. 


MEASUREMENT ANODE NOISE CURRENT FLUCTUATIONS 


The experiment consisted raising the cathode temperature and observing 
the magnitude the noise power coupled resonant cavity the gun 
anode, the beam current increased its space-charge limited value. The 
form the variation noise with d-c. beam current shown Fig. Two 


RELATIVE NOISE 


Fic. Electron beam noise the gun anode function current. 
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distinct regions the curve are apparent. low beam currents the diode 
gun was temperature-limited and the initial noise fluctuations the beam 
were unmodified the electron gun anode. the cathode temperature was 
raised, the beam current increased its space-charge limited value, and the 
greater charge density within the gun resulted decrease the anode 
noise current fluctuations. Extrapolation the linear portion the curve 
serves determine the full shot noise content the beam for any current 
value. The noise smoothing defined from the ratio measured noise power 
space-charge limitation, that value shot noise indicated the shot 
noise asymptote for the space-charge limited current. This form normaliza- 
tion obviates the necessity absolute noise power measurements. 


CALCULATION ANODE NOISE CURRENT SMOOTHING 


Three methods have been used compute the smoothing the anode 
noise current below full shot noise. For the latter two, values the current 
density and the beam diameter were required. space-charge limited 
planar diode, the current density given by: 


3/2 
(1) 


where the spacing and the gun voltage. effective 
beam diameter was found comparing the observed perveance with the 
theoretical current density. The mean perveance tubes D-2 and I-1 was 
and this was considered representative properly 
operating gun. This indicates that the effective beam diameter was 2.3 
in. which 73% the aperture the beam forming electrode. 


(a) The Llewellyn—Peterson Equations 

The computation was started the potential minimum. Full shot noise, 
uncorrelated Rack velocity fluctuations, and d-c. velocity equal the 
r.m.s. thermal velocity were the assumed input conditions. Using the co- 
efficients quoted Pierce (1950), can shown that the anode noise 
smoothing given 


(2) Van 8 Vie (4 m 


where the voltage equivalent the r.m.s. thermal velocity and 
given the cathode temperature assumed 1160° The 
first term (2) negligible high frequencies, since 0.1 volt. 


Computer Solution 

Kornelsen (1957), developing the finite beam solution, specified input 
conditions just those indicated above for the Llewellyn—Peterson equations. 
Using the general result for determined the previous paper (McFarlane 
1961), and noting the potential minimum, the smoothing 
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0a 
determining and the computer started point corresponding 
zero transit time. 


(c) Computer Solution 

This calculation has been discussed detail the companion paper 
(McFarlane 1961). The computation was started plane front the 
potential minimum using input quantities the noise parameters found 
Siegman al. (1957). The correlation coefficient was determined 
requiring the conservation real kinetic power. Its actual value only slightly 
affects 


COMPARISON THEORETICAL AND EXPERIMENTAL RESULTS 


The experimental values noise smoothing observed the gun anode 
are compared with the results the three theoretical determinations 
Figs. and The 1,400 Mc/s measurements tubes D-1 and D-2 agreed 
well within experimental error and they have been shown one measure- 
ment. 4,250 Mc/s only slight differences between the tubes were observed, 
while 9,250 Mc/s two distinct sets measurements appeared. 

The predictions the Llewellyn—Peterson equations and the first com- 
puter solution are strong functions frequency and, would anticipated, 
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Fic. Smoothing function gun voltage 1,400 Mc/s. 
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4,250 Mc/s 


Computer 


Llewellyn-Peterson 


Experiment Tube 


Tube No,0-2 


200 400 600 
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Fic. Smoothing function gun voltage 4,250 Mc/s. 
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Fic. Smoothing function gun voltage 9,520 Mc/s. 
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the agreement between them improves the frequency increased. Both 
indicate more rapid change smoothing with gun voltage than observed 
experimentally. 

The second computer solution predicts variation smoothing with gun 
voltage better agreement with that observed. all frequencies the theo- 
retical smoothing consistently greater than that observed, the discrepancy 
being somewhat larger high gun voltages. 

examine the experimental measurements for the theoretical dependence 
determined the companion paper, the observed smoothing 
values have been used calculate and the results are presented 
Table II. Considering that the measurements extend over frequency 
range 6.8:1, the data not show any serious deviation from the w!/? 
variation predicted. The increase the gun voltage was raised 
attributed the progressively higher cathode temperature required 
space-charge limit the gun. the cathode temperature was increased, the 
quadratic velocity input was increased, resulting higher value the 
anode current fluctuations. Zacharias and Smullin (1959) have examined the 
effect cathode resistance the noise parameters the emitted beam and 
find that results increase the quadratic velocity content and 
the presence negative real kinetic power. While not exactly clear how 
these effects modify the input conditions assumed Siegman, the net result 
further increase anode noise. 

The observed noise current fluctuations are higher than those pre- 
dicted theory. The difference cannot attributed interception noise. 
(1958, 1959) showed that aperture interception results increased 
beam noise only when relatively large fraction the beam intercepted. 
view the agreement between the zero and high magnetic field smoothing 
data also concluded that the excess noise does not result 
from the lens effect the anode aperture (Knechtli 1958). 

The experimental results Vessot (1957) and (1958, 1959), 
using higher perveance electron gun similar design, have been analyzed 
discover similar behavior the noise was present. The data pre- 
sented Table 

The first three columns are from Vessot’s measurements and the final 


TABLE III 
Values measured Vessot and McIntosh 


62 Voa/w'/? 


voltage 230 Mc/s 1,010 Mc/s 3,055 Mc/s 3,050 Mc/s 
700 


- 

“al 
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column based MclIntosh’s oxide cathode data. The author has made 
theoretical calculation for this gun. The dependence the 
noise current smoothing frequency and gun voltage reasonable agree- 
ment with the form predicted for the present electron gun. 


CONCLUSIONS 


The accurate description the transport fluctuation quantities current 
and velocity between the potential minimum and anode space-charge 
limited diode electron gun requires that account taken the non-con- 
servative nature the electron stream immediately beyond the potential 
minimum and the effect finite beam diameter the high velocity region. 
Using the results density function calculations specify the beam noise 
properties plane the gun beyond the potential minimum, possible 
determine the anode noise parameters means computer solution 
the finite beam electronic equation. 

The theory indicates dependence the anode noise current smoothing 
the one-half power the frequency and the inverse first power the 
gun voltage. The expected was found the measurements here 
reported and those other workers. 
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PHASE DIAGRAM VERY LOW 
THALLIUM CONCENTRATIONS! 


McLaren? AND WEINBERG? 


ABSTRACT 


The liquidus has been accurately determined from pure (419.505° 
the monotectic transition temperature (416.926° 0.42 at.% using 
precision resistance thermometry. The upper solidus was determined approxi- 
mately from measurements the distribution coefficient (~0.01) and the solid 
solubility limit (~0.004 at.% thallium zinc. value kcal/mole 
for the latent heat fusion pure zinc was calculated from the freezing point 
depressions. 


INTRODUCTION 


investigation has been made the phase diagram from pure zinc 
the monotectic composition part study solute segregation 
ingots very dilute zinc alloys. 

Heycock and Neville (1897), using resistance thermometry, determined the 
monotectic reaction temperature the system 416.5° and the 
monotectic composition 0.45 at.% They also determined two points 
the hypomonotectic liquidus 0.126 and 0.393 at.% TI. the present investi- 
gation, the liquidus and monotectic points were accurately determined 
resistance thermometry for series alloys the concentration range 
0.0002 0.88 at.% addition, the upper solidus was determined approxi- 
mately from measurements the distribution coefficient and the solid solu- 
bility limit thallium zinc. value kcal/mole for the latent 
heat fusion pure zinc was calculated from the freezing point depressions 
and was critically compared with Kubaschewski’s (1950) selected value 
kcal/mole. 


EXPERIMENTAL 


Alloy Preparation 

The alloys (1030-g samples) were prepared melting weighed amounts 
zinc* and 99.99% thallium under dry nitrogen atmospheres high 
purity graphite crucibles. The samples were held 650° for days com- 
plete solution. Certain the very dilute alloys contained 0.77 Mev; 
half life, years) tracer order locate the thallium for segregation 
studies and measure thallium concentrations for distribution coefficient 
determinations. 


received December 15, 1960. 
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*Both New Jersey Super Pure and Chemically Pure grades zinc were used solvent (see 
McLaren 1958). 
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Liquidus 

(a) Method 

The apparatus and techniques used for high precision thermal analyses 
several high purity metals platinum resistance thermometry have been 
fully described McLaren (1957, 1958, 1959). The alloys, contained 
graphite crucibles, were solidified under dry nitrogen atmospheres well- 
lagged metal block furnace. The platinum resistance thermometer was inserted 
into pyrex thermometer wells that were axially positioned the cylindrical 
melts. Freezing techniques have been developed eliminate greatly reduce 
the uncertainties caused supercooling liquidus point determinations 
very dilute alloys. 

During normal freeze, the alloy freezes radially from the outside the 
center the crucible rate controlled the over-all furnace cooling rate. 
very dilute zinc alloy supercools much 0.04° before solidification 
begins, then the temperature the alloy slowly rises toward the equilibrium 
liquidus temperature. This slow recalescence effectively masks the primary 
arrest the cooling curve and contributes uncertainty the liquidus 
point determination. overcome difficulty, freezing induced 
removing the thermometer from the well when the alloy becomes supercooled, 
allowing cool, and then replacing the well. This operation results 
the rapid nucleation and growth thin mantle solid metal the ther- 
mometer well which releases sufficient latent heat raise the melt immediately 
its equilibrium liquidus temperature. Thereafter, solidification proceeds 
normally from the crucible walls. 

The first solid freeze, immediately after overcoming the supercooling 
this induced freezing technique, has the solute concentration that closest 
the equilibrium concentration the solid the liquidus point the alloy. 
Therefore, the solute concentration gradient this first solid primarily 
determines the accuracy that attainable the temperature measurement 
the liquidus point. The most important factors setting this concen- 
tration gradient are the over-all solute concentration the alloy, solute 
segregation coefficient, rate freezing, and the slope the liquidus for the 
alloy system. 

Figure illustrates how p.p.m. thallium (containing some tracer) 
zinc segregates during slow induced freeze. This figure includes macro- 
photographs and autoradiographs horizontal and vertical sections cut from 
the ingot. The autoradiographs show that segregation freezing for this 
solute zinc great that most the appears swept into thin 
cylinder between the mantle and the main freezing interface. 

Since the solute concentration gradient the first frozen solid this alloy 
very small, the melt temperature slow induced freeze remains nearly 
constant for large fraction the total freezing time:* this enables precise 
temperature measurements readily made. The highest precision 
attained for very dilute alloys because these ingots solidify primarily 


*Typical cooling curves very dilute zinc alloys are shown McLaren (1957, 1958). 
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large grained outer shell advancing with smooth interface into the remaining 
liquid: this closely approximates the freezing single crystal with plane 
interface advancing into homogeneous liquid. 

For less dilute alloys and the same freezing rates the over-all grain size 
the solidified ingots decreases and zone equiaxed grains that nucleate 
the unfrozen liquid appears. result, significant amounts solute become 
trapped intergranular regions and the relation between melt temperatures 
and liquidus temperatures becomes less well defined. 

Table gives the reproducibilities the temperature measurements that 


TABLE 
Precision temperature measurement 


Conc. thallium Reproducibility 
(wt. liquidus temp. (°C) 
0.0001 0.0002 
0.001 0.0002 
0.01 0.0005 
0.1 0.002 
1.0 0.02 
1.3 (monotectic) <0.01 


were found repeated liquidus point determinations the dilute 
alloys. The experiments were carried out with interface growth rates between 
0.2 and 0.5 cm/hr. 

Results 

this investigation, liquidus point determination was made first the 
zinc solvent; then the thallium alloy was prepared and its liquidus point 
determined. Temperatures were related the liquidus point 419.505° (Int. 
1948) sample (S.4) New Jersey Super Pure Zinc having alloy melting 
range about 0.001° (see McLaren 1958, 1959): Table summarizes the 
freezing points that were measured. 


TABLE 


Liquidus and monotectic temperatures alloys related the 
liquidus point (Zn) sample (S.4) high purity zinc 419.505° 


(Int. 1948) 

liquidus monotectic Liquidus slope 
0.00075 0.00023 0.0012 
0.0059 0.0018 0.0130 
0.0071 0.0022 0.0139 2.6 —6.32 
0.0567 0.0177 0.109 2.6 
0.283 0.0885 0.571 2.6 
1.01 0.316 1.92 2.58 
1.36 0.425 2.56 2.58 
1.38 0.431 2.56 2.58 —5.94 
1.71 0.535 2.58 
2.06 0.644 2.58 
2.82 0.882 2.58 


Mean slope —6.18° C/at.% 
Mean deviation +0.37° C/at.% 


plane 


ize 
come 
tures 


that 


ween 


the 
(Int. 
the 


PLATE 


sections cut from ingot p.p.m. alloy following induced freezing (ingot diameter 
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(a) Micrograph polished and etched horizontal section cut from monotectic 
alloy showing the cross section from crucible wall, thermometer well, (radial 
length ~12 mm). 


(b) Micrograph showing fine spheroidal dispersion thallium-rich phase the fast- 
frozen induced mantle, 


PLATE 
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Solidus 

locate the solidus, measurements were made the distribution co- 
efficient, thallium zinc, defined the ratio where the 
concentration the solid that separates the liquidus temperature the 
alloy, and the concentration the initial liquid alloy. Two methods 
were tried. 

First, the ratio the relative concentration first-frozen solid 
the initial homogenized concentration was determined following controlled 
solidification rods alloy described Yue and Clark (1960). 
Turnings for beta counting were machined from the rods and pressed into 
standard samples for and counting. 

Secondly, the ratio was determined for alloys that were solidified 
the regular cylindrical ingots following induced freezing. Turnings were 
machined from the outside radius the ingots obtain samples the first- 
frozen material, C,, following recovery from supercool. this case, samples 
(15 the homogenized liquid, Co, were extracted and water-quenched 
pyrex ladles prior the ingot solidification; cuttings made across full 
section were used for standard counting samples. 

The samples used both methods were obtained from material solidified 
rates 0.2-0.5 cm/hr. 

For alloys initial concentrations the range 100 p.p.m., the measured 
values ranged from 0.006 0.02, with average value about 0.01. 
correlation between measured value and alloy concentration was not found. 
Large variations (up 30%) were found the counting rates selected 
samples and material from the same alloy. The presence micro- 
segregations thallium these samples was shown autoradiographs 
sample turnings. Smaller variations (10%) counting rate were found among 
samples material cut from the fast-frozen samples removed from the 
regular cylindrical ingots before solidification. attribute the uncertainties 
our determinations the distribution coefficient thallium zinc these 
microsegregations, beta counting being particularly sensitive self-absorption 
variations non-homogeneous samples. 

Thermal analyses the hypomonotectic alloys showed that the solid 
solubility limit thallium zinc must lie between 0.0071 and 0.057 wt.% TI. 
Prominent monotectic arrests (5% total time) were found melting and 
freezing (total time hours) 0.057 wt.% alloy. contrast, monotectic 
arrests were imperceptible thermal curves 0.0071 wt.% alloy these 
rates fusion and were barely detectable when cooling and heating rates were 
further reduced (total time hours). The monotectic arrests revealed the 
latter case were undoubtedly due non-equilibrium divorced monotectic. 
The solid solubility limit being this concentration range implies that the 
distribution coefficient thallium zinc the monotectic temperature 
the range 0.005 0.04. This general agreement with the results the 
concentration measurements made greater dilution. 

consideration the above measurements, have arbitrarily selected 
0.01 representative value the distribution coefficient thallium zinc. 
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Monotectic 

The monotectic composition, 0.42 at.% was calculated using the values 
determined for the mean slope the liquidus and for the depression the 
monotectic temperature below the freezing point pure zinc. This mono- 
tectic reaction temperature was determined 416.926+0.003° with 
pressure coefficient for atmosphere over the melt. This 
pressure coefficient identical with that previously found for pure zinc (see 
McLaren 1957). 

Figure shows micrograph sectioned monotectic alloy following 
induced freezing. The thallium-rich phase dispersed spheroidal state 
the zinc-rich matrix. Figure shows micrograph the region near the 
thermometer well: the thallium-rich phase appears finer dispersion the 
fast-frozen induced mantle than the slowly frozen main portion the 
ingot. 


Phase Diagram 

Figure shows the phase diagram for dilute determined this 
investigation. The liquidus was constructed using the freezing points given 
Table and the solidus was approximately located using distribution 
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coefficient 0.01 for thallium zinc. Freezing points determined the 
system Heycock and Neville (1897) and Honda and Ishigaki (1925) 
are included the diagram Fig. 


DISCUSSION 


The liquidus, within experimental error, appears vary linearly with 
thallium concentration from pure zinc the monotectic composition. The 
major uncertainty locating the liquidus arises assigning composition 
values for the series alloys. The amount thallium lost evaporation 
and incomplete solution during the alloy preparation not known. However, 
since the liquidus points the alloys have high reproducibility after pre- 
paration (Table and generally lie very close smooth curve drawn 
between the freezing temperature pure zinc and the monotectic temperature, 
probable that the concentration uncertainty the liquidus drawn 
Fig. inside +5% over most the range investigated. doubtful that 
whole ingot chemical analyses could made with higher precision, especially 
for the extremely dilute alloys. 


Latent Heat Fusion Pure Zinc 

Values for AH, the latent heat fusion pure zinc, can obtained from 
standard thermodynamic relations. 

1939, Hayes and Chipman derived form Van’t Hoff’s law relating 
the freezing point depressions very dilute metal solutions: 


(1) 


where the freezing point the pure metal, the solute distribution 
coefficient, the slope the liquidus, and the gas constant. For our 
experimental results this relation gives value 


1.53+0.1 kcal/mole, 


the uncertainty arising mainly from uncertainties solute concentrations. 
alternative expression containing and involving the pure metal 

only, that which Lumsden (1952) relates the freezing point, its 

pressure coefficient, and the decrease molar volume solidification: 


McLaren (1957) determined the pressure coefficient the freezing tempera- 
ture pure zinc +0.0043+0.0001° for atmosphere. The value 
commonly quoted (0.27 cc/mole (2.8% shrinkage)) that Pelzel and 
Sauerwald (1941), while Saeger and Ash (1932) obtained 0.36 cc/mole (3.8% 
shrinkage). These lead values 1.1 kcal/mole and 1.4 kcal/mole respec- 
tively. Because the uncertainty the value AV, equation (2) sub- 
stantially less precise than equation (1) means obtaining AH, but 
can combined with the result equation (1) give value 0.39 
cc/mole moderately close agreement with Saeger and result. 
Previous determinations the latent heat fusion pure zinc have been 
summarized Kubaschewski (1950, 1958) review paper the latent 
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heats fusion metals; selected value kcal/mole for for 
pure zinc. This value appears based partly his own calorimetry, for 
which quotes precision about 6%, and rather strongly heats 
fusion calculated Lumsden (1948) from Heycock and Neville’s (1897) 
freezing point depressions (solutes Sn, Bi, Tl, Pb, Sb, Mg), and Honda 
and Ishigaki (1925) from their own freezing point investigations (solutes Cd, 
Tl, Ag, Al). 

Both groups used alloys synthetic composition and the purity the zinc 
solvent not specified. certain that melt temperatures could not attain 
equilibrium liquidus temperatures during recalescence Heycock and 
Neville’s experimental arrangement which utilized massive resistance thermo- 
meters and small (250 ingots. Also, because the immersion effects* found 
Heycock and Neville (1895) earlier investigation the freezing 
point zinc, can clearly deduced that the 250-g ingots used their 
alloy investigations would provide insufficient overcome serious 
stem losses. Honda and Ishigaki’s thermal analyses were made only two 
alloys each solute system: they used thermocouple thermometry and cooling 
rates that estimate 10-100 times greater than those used this 
investigation. They report liquidus point depression 2.70° from pure 
zinc for alloy 0.50 at.% thallium: this depression 0.12° greater 
than the depression the monotectic reaction temperature found this 
investigation. 

Because these considerations, the better than precision attached 
Kubaschewski’s selected value not justified and the former 
1.6 kcal/mole for the heat fusion zinc selected from earlier literature 
survey Kelley (1936) should not have been discounted Kubaschewski’s 
review. surprising that these large uncertainties§ are still present the 
accurate evaluation simple thermochemical constant. 

The results this investigation dilute alloys support the lower 
value selected Kelley (1936) and the higher value the volume 
shrinkage solidification determined Saeger and Ash (1932). Additional 
work contemplated the establishment the binary phase diagrams 
extreme dilution for other solutes zinc; this will provide additional informa- 
tion for the determination the latent heat fusion pure zinc. 
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*Freezing temperature decreased 0.78° reducing ingot size from 700g (dia- 
meter mm). 

difficulty has been overcome our investigations the use large ingots and 
modern Meyers thermometers (see McLaren 1959). 

value for was also selected Rossini al. (1952) their compilation thermo- 
dynamic properties. 

§Some the uncertainty found the determination these thermodynamic properties 
may the result of, yet, unaccounted effects orientation, lattice defects, and porosity 
the solid phase. 
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EFFECTS EQUATORIAL SPREAD-F IRREGULARITIES 


ABSTRACT 


study made into the effects spread-F conditions continuous wave 
(c. w.) transmissions between Colombo and Waltair over distance about 
1300 kilometers. shown that under suitable conditions increased fading 
rate the received signals sufficient indication spread and gives 
genuine index its intensity. Pulse and methods are used simultaneously 
study the time variations spread-F intensity separated points. The 
results indicate that the extent spread-F occurrence the 
direction greater than 650 km. 

The probability distributions signal amplitudes the fading records 
were closer Rayleigh type under conditions than during spread-F 
conditions. suggested that the flutter fading caused the presence 
spread-F irregularities the appropriate zone reflection the F-region 
the ionosphere. 


INTRODUCTION 


The strong positive correlation reported between the radio star scintillations 
and spread (Ryle and Hewish 1950; Little and Maxwell 1951) suggested 
the occurrence the required irregularities extensive patches consider- 
able horizontal extent even (Briggs 1958). has been shown 
that the fluctuations sometimes observed the intensity the radio waves 
received from the direction some star constellations, are due their passage 
through the ionospheric irregularities. Effects similar this should also 
observable the case transmissions through the F-region spread-F 
irregularities are prevailing the relevant zone reflection. The present 
paper concerned with the study some these effects. With this end 
view signals from the regular broadcasting station Colombo have been 
picked and the signal fading observed studied relation the simul- 
taneously observed overhead spread-F conditions Waltair (geographic 
latitude 17°41’ N.; geomagnetic 7.4° N.). 


SCHEME OBSERVATIONS AND METHOD ANALYSIS 


Continuous wave signals 11.717 Mc/sec from Colombo (Ceylon) situated 
about roughly the south Waltair, the receiving site, were 
recorded using conventional equipment consisting communications type 
SX-42 receiver, d-c. amplifier, and Esterline-Angus pen recorder. About 
100 these fading records, each minutes’ duration, were taken 
during the months January and February, 1959, different times the 
night between 1800 hours and 2300 hours Indian Standard Time (I.S.T.), 
interval during which the transmissions are available. Simultaneously with 
each these recordings, spread-F conditions are observed Waltair using 


received November 1960. 
Contribution from the Ionospheric Research Laboratories, Andhra University, 
India. 
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the conventional vertical incidence ionospheric pulse sounding equipment. 
Each time the group height range (in km) over which the echo spreads 
has been recorded visual observation express the intensity local 
spread-F conditions terms the indices 0-10 accordance with the 
scheme defined Gopala Rao al. (1960). all these 100 occasions, pulse 
soundings were made only frequency 6.4 Mc/sec, which found 
nearly the equivalent vertical incidence frequency corresponding the 
first hop reflection the present transmissions. Exact evaluation 
this frequency will not necessary important since the equatorial spread-F 
characteristics would not differ much with the operating frequency. Examina- 
tion the and maximum usable frequency (m.u.f.) (3000 km) data 
Trichinopoly, and Trivandrum stations (All India Radio which are 
the vicinity the second and third hop reflection points the F-region, 
indicated that the frequency 11.717 Mc/sec was generally greater than 
the optimum working frequencies for the and propagation 
modes between Colombo and Waltair. Thus the mode expected 
mostly predominant these transmissions and gives information about the 
ionospheric F-region conditions location 650 south Waltair (near 
Madras). 

Fading rate cycles per minute each the records defined 
here the number maxima occurring per minute. order keep con- 
sistency only fluctuations greater than 10% the average amplitude 
fading were considered. 


RESULTS PRELIMINARY ANALYSIS 


striking feature most the records consistent increase the 
fading rate whenever spread occurred Waltair. This consistency itself 
suggested the occurrence spread-F conditions 650 south Waltair and 
also where the transmissions from Colombo encounter the F-region before 
reaching Waltair. Despite the simplicity equipment, well known that 
the method serious limitation that the waves received different 
group paths are not resolved. Observations spread-F echoes the 
method can therefore neither simple nor reliable the case with the 
pulsed radio echo method. The present series simultaneous observations 
have positively indicated that, under the present experimental conditions, the 
enhanced fading rate genuine representation spread-F conditions. 
Typical fading records taken January 18, 1959 2300 hours when 
there was spread Waltair and the one taken January 21, 1959 
2030 hours I.S.T. when spread-F index 6.4 Mc/sec Waltair are 
shown Figs. and respectively. 

All the fading rates the records connected with spread-F index 
Waltair are separately averaged and the rest are grouped into ranges 
with 1-2, 3-4, etc., spread-F indices. The average fading rates con- 
nected with each these groups are shown Table The average 
fading rate was only 10.7 cycles/minute under spread-F condition 
Waltair. Three- four-fold increase the average fading rate observed 
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TABLE 
Results preliminary analysis 


Spread-F fading rate 
index group cycles/minute 

10.7 
1-2 30.1 
3-4 33.8 
5-6 42.6 
7-8 34.8 
9-10 40.8 


whenever spread any degree occurred Waltair although systematic 
increase could observed the fading rate with increasing degree 
spread recorded Waltair. these circumstances fading rate cannot 
used measure the degree spread although general there 
positive indication that the enhanced fading rate caused spread-F 
irregularities. This discrepancy explained the following section con- 
sidering the possible time differences the occurrence certain intensity 
spread the overhead point Waltair and 650 south Waltair. 


TIME SHIFTS THE TEMPORAL VARIATIONS SPREAD 
SEPARATED POINTS 

second series observations during the months January and Febru- 
ary, 1960, continuous recordings together with overhead pulse observa- 
tions were closer time intervals between 1900 hours and 2300 hours 
I.S.T. Further, this series vertical incidence spread-F observations were 
made group spot frequencies which included 6.6 Mc/sec one 
the spots cover the frequency range over which spread-F echoes are 
observable. The variations, with time, the fading rate signals from 
Colombo and the spread-F index Waltair averaged all the spot 
frequencies well the index 6.6 Mc/sec alone obtained some 
the nights are illustrated Figs. and These figures readily show 
reasonable similarity the variations together with observable time shift 
expected. also clear that the time variations the average spread-F 
index and that 6.6 Mc/sec alone are very much similar indicating the 
uniformity spread-F index time variations different frequencies. The 
lack systematic increase found the values Table can now attributed 
the grouping the simultaneously observed parameters. the other 
hand shift either the curves such way that the nearest similarities 
the maxima and minima coincide, considerable correspondence evident 
the variations. plotting the observations the fading rates around the 
maxima and minima against the values the corresponding portions the 
maxima and minima the spread-F index variation curves ignoring the time 
shift, the scatter plot shown Fig. obtained. will seen from this 
figure that fairly good linear dependence exists between the fading 
rate and the average spread-F index. The existence systematic variation 
the fading rates with the changes the degree spread thus 
revealed. 
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2-3. Time variations fading rate signals from Colombo and spread-F 
index Waltair. 
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C.W. SIGNAL FADING RATE C/MIN 


SPREAD INDEX WALTAIR 


Fic. rate signals from Colombo versus spread-F index Waltair (time 
shifts ignored), 


Examination some the similarities, shown Figs. and arrows 
and several other similar variation curves, indicated earlier occurrence 
certain maximum minimum, the fading rates with respect those 
spread-F index variations Waltair. The distance and the 
observed time shifts may utilized obtain quantity with the dimensions 
velocity. Considering the average time shift corresponding the first 
major peaks, which about minutes, velocity 540 m/sec obtained. 
This represents only the apparent velocity the N.-S. direction. However, 
does not seem reasonable associate this with any movement drift 
because there some evidence (unpublished) believe that the timing and 
magnitude the first peak spread-F index closely related with local factors 
such height changes the evening hours. The first peak spread-F index 
each place seems later according the rate increase 
the early evening hours after sunset large small. this rate 
generally largest near the magnetic equator and smaller successively 
points farther from the equator, the peak spread-F index occurs generally 
earlier records when compared observations overhead point. 
Hence the observed time shifts are not interpreted due any genuine drift 
irregularities. 

this connection interesting and useful consider the results 
horizontal drifts Waltair reported Ramachandra Rao and 
Bhagiratha Rao (1959). about 2000 hours I.S.T. the drift reported 
mainly towards the north with average magnitude m/sec. 
assume that this drift pattern unaffected spread-F conditions and 
spread-F irregularities generated one location are carried these drifts, 
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will take least hours minutes for travel 650 km. The secondary 
maxima the spread-F variation Waltair denoted ‘S’ Figs. and 
occurring (approximately) around hours after the first major maximum 
the fading rate variations might have been caused the northward 
component the general F-region drifts. 


COMPARISON WITH FLUTTER PHENOMENON 


Sometimes along with the rapid fading under conditions quivering 
noise superposed the programme heard the loud-speaker giving 
impression that still higher fading rates might present than could 
recorded the instrument. likely that the amplitudes such fast- 
fading components are smaller than could recorded considered the 
present analysis which the fluctuations less than 10% the average 
amplitude fading were neglected. 

These experimental conditions remind one the ‘flutter phenomenon’ 
reported Subba Rao and Somayajulu (1949). They reported peculiarly 
rapid fading signals from Madras and Colombo transmissions received 
Waltair during the nights the months February June, 1948. This 
period corresponds roughly the same part the sunspot cycle the current 
investigation period. They reported that the quasi-flutter frequency was low 
(20-30 cycles/minute) before sunset and increased cycles/minute 
during night hours. mention was made them about any amplitude 
restrictions the evaluation the fading rates. However, from the fading 
records published them clear that this fast-fading flutter component 
much smaller amplitude than the average amplitude fading. The 
background the present investigation suggests that this flutter phenomenon 
might have been caused spread-F conditions prevailing the appropriate 
zone reflection. may noted that recently Wright (1959) also reported 
several similarities between the flutter the transmissions and spread-F 
phenomenon. 


AMPLITUDE DISTRIBUTION ANALYSIS FADING RECORDS 


The fading patterns the records this investigation have been 
generally random except few occasions when some kind quasi-periodic 
fading occurred. From observation large number these records, 
found that the depth fading somewhat reduced when transmissions 
are suspected affected spread-F conditions against normal con- 
ditions. order study this and other features, the probability distributions 
the signal amplitudes were studied for some the records. The deflection 
the pen recorder pointer varied linearly with the input signal voltage 
the receiver. The amplitude the signal terms the chart divisions was 
determined series equally spaced time intervals 0.625 second. The 
probability distribution the amplitudes was obtained. This was compared 
with the Rayleigh type distribution given 
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where ‘R’ any amplitude the observed range, and being the 
most probable amplitude. Substituting the experimental value for the 
theoretical values P(R) were calculated for various values ‘R’ and were 
multiplied suitable constant obtain best fit. comparison 
was found that distributions were nearer the Rayleigh type during 
normal conditions than when there was spread associated with rapid fading 
the records. Some the typical distribution curves are shown 
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Fic. Typical amplitude distributions fading transmissions from Colombo 
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McNicol (1949) has shown that Rayleigh distribution only special case 
which occurs when the downcoming waves consist only large number 
scattered waves approximately the same group path but with phase paths 
differing from one another random fractions wavelength. According 
his findings the deviation from the Rayleigh type observed this investi- 
gation indicates the presence steady component addition the ran- 
domly fading component. This leads the fact that under most severe 
conditions spread the depth fading reduced the proportion the 
steady component the signal increases. This consistent with the smaller 
amplitude the flutter fading referred earlier and supports the view that 
flutter phenomenon caused spread-F irregularities. interest 
note that similar situation namely the reduction the depth fading 
under fast-fading conditions has recently been reported Koch al. (1960) 
connection with transmissions through auroral regions. This leads 
some similarity between spread-F and auroral conditions regarding the nature 
the irregularities. 

Spread-F irregularities are therefore observable and can studied even 
with simple equipment under suitable conditions when periodic fading 
due other causes does not interfere vitiate the observations. Results 
the present investigation suggest that patches spread-F irregularities extend 
horizontally scale greater than the direction this 


location. 
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NOTES 


GAMMA DOSE HOLE UNIFORMLY CONTAMINATED PLANE: 
CONTRIBUTION GROUND PENETRATION* 


CLIFFORD 


assessing the protection against y-radiation provided hole 
contaminated field necessary consider the contribution the total 
dose produced y-radiation that penetrates the ground from contamination 
around the hole. the diameter the hole decreases this component becomes 
comparable the component due radiation backscattered from the sky 
skyshine. Usually the effect penetration ignored since can 
eliminated simply removing the contamination within few feet the 
edge the hole. However, when this contamination not easily removed— 
for example, the case foxhole was buried when the foxhole was 
dug—the ground penetration component could significant fraction 
the total dose. This ground penetration component could also important 
the case corner basement shelter particularly contamination was 
washed off the roof and piled against the outer basement wall. 

recent work the Defence Research Chemical Laboratories (Clifford 
1960) measurements were made the dose distributions particular hole 
ground penetration y’s from uniform contamination the surrounding 
plane. Some the more important results these measurements are given 
and compared with calculations from the work y-scattering carried out 
the National Bureau Standards Spencer. 

The cylindrical hole studied was meters diameter and meters deep. 
The wall, floor, and lip the hole were made sand and cement concrete. 
was chosen the contaminant since the y-ray energy (0.66 Mev) 
close the average energy the y-rays emitted fission few 
hours after detonation weapon. All dose measurements were made 
combining the results single source determinations give the values that 
would have been obtained with uniform contamination. During exposure 
the source was placed slight depression the concrete that 
its effective center was level, and covered with lead shield 
eliminate any skyshine dose the detectors the hole. 

The depth variations dose due ground penetration both along the 
axis the hole and along line halfway between the axis and the wall 
the hole are shown Fig. for uniform ground contamination 
The distributions down the center and halfway the wall are quite 
similar. The dashed curve gives the calculated dose along the central axis 
and was obtained from the work Spencer. lower depths agreement 
good while shallow depths the experimental values are lower than those 
calculated. 

*Issued DRCL Report No. 310A. 
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CENTRAL AXIS 
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EXPERIMENTAL 


DEPTH HOLE 


Fic. Dose rate due ground penetration y-radiation into hole concrete plane 
uniformly contaminated with Diameter hole meters. 


Figure shows the dependence the ground penetration dose the area 
contamination for two depths the central axis. Over 90% the dose 
due contamination within foot the edge the hole and was found 
that contamination the region beyond centimeters contributed less than 
1/10th the dose received from the first centimeters. For com- 
radius gives dose 12.1 mr/hr and over area 200 meters radius 11.5 
mr/hr height meter above the surface. the latter case skyshine 
doses 2-meter diameter hole meters deep are 0.74, 0.46, and 0.22 mr/hr 
depths 0.25, 0.50, and meter (Clifford 1959). depth meter 
the total dose (skyshine plus ground penetration) therefore approximately 
the dose meter above the surface and this roughly one-third 
due ground penetration which can eliminated clearing away the 
contamination within feet the edge the hole. For the case smaller 
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Fic. dose rate due ground penetration extent contamination. 
Dose rate measured axis 2-meter diameter hole plane contaminated with 
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diameter hole, elimination the lip penetration component would give 
even greater improvement shielding since the skyshine component more 
dependent than the lip penetration the diameter the hole. 

The individual dose measurements that had been integrated study the 
effect uniform contamination were analyzed order study the effect 
lip penetration detail and see whether not elementary method 
could used for dose calculations this and other shielding problems involv- 
ing lip geometry. Since infinite media data are available for many materials 
and energies (Goldstein 1954; Fano 1959) was felt that would useful 
relate the observed values this particular experiment corresponding 
known infinite media data. For the limiting case detector close the 
lip and below ground level with source well away from the lip, appeared 
reasonable expect that the detector would receive approximately one-half 
the scatter dose that would receive the corresponding lip penetration 
distance infinite concrete medium. Its dose build-up factor would then 
where B,, the dose build-up factor for infinite medium. 
Fig. summary curves the observed dose build-up factors are expressed 


~ 


LIP PENETRATION DISTANCE CONCRETE 


Fic. dose build-up vs. lip penetration distance, for lip pene- 
tration concrete hole using sources. 


fractions and are plotted against the concrete lip penetration 
distances mean free paths (m.f.p.). The data, total 230 points from all 
detector locations the foxhole, were arbitrarily divided into four approxi- 
mately equal groups depending the direction angle the source detector 
direction below the horizon. The lower curve summarizes the readings with 
20°, the others with between 20° and 30°, 30° and 40°, and 40°. For 
each group the average deviation point from the curve was 43%. From 
these curves would appear that the dose due ground penetration for 
shielding case involving lip geometry circular hole could approximated 
taking for distances less than m.f.p. shallow depths this 
elementary approach would increasingly error but more than factor 
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study lip geometry build-up factors further detail, laboratory measure- 
ments were made the dose transmitted between points intersecting sides 
cube homogeneous absorber, lucite+36% sucrose solution, 
using similar techniques the hole experiments. 

Source and detector positions were chosen that the dose fixed lip 
the angle penetration from the source plane. The observed dose build- 
factors were expressed fractions one-half the infinite medium build-up 
factor and are shown Fig. function The ratios are 
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Fic. dose build-up vs. penetration angle for lip penetration 
rectangular lucite tank using sources. 


grouped around value large values i.e. source close the lip 
and detector away from the lip, the ratios decreased for the and m.f.p. 
cases. these shorter lip penetration distances considerable amount the 
scattered radiation that reached the detector was generated material close 
the source. increases the amount this material becomes less and 
not unreasonable that the build-up should decrease. greater lip pene- 
tration distances, radiation scattered material close the source less 
likely reach the detector because the greater absorption lower energy 
radiation and decrease build-up with increasing less likely. 

The ground penetration dose hole uniformly contaminated plane 
can obtained integration over distribution point sources. For 
the foxhole studied, the average deviation, from the observed values the 
central axis, the ground penetration dose calculated using was 
about depths 50, 75, 100, and 125 centimeters. the most shallow 
depth measured, centimeters, the deviation was 14%. The use 
calculating the dose produced individual point source probably would 
error lip penetration distances greater than m.f.p. for rectangular 
lip lucite m.f.p. for circular hole concrete. However, discussed 
above, the case foxhole shielding the contribution from radiation that 
has penetrated more than m.f.p. would contribute only small error. 
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DIMENSIONS DES GRAINS COURS PRECIPITATION 
DES EMULSIONS IONOGRAPHIQUES 


Montréal, nous avons examiné avec microscope électronique 
RCA-EMU-2C des échantillons d’émulsion précipités diverses vitesses 
d’agitation (Demers 1958; Fournaux, Demers Demers 1960). L’agitation 
plus lente choisie permet peine d’éviter stagnation des nappes 
diverses densités; plus rapide évite peine formation bulles d’air. 
Pendant ces essais vitesse des réactifs dans méthode 
double jet appliquée, est restée comprise entre les limites choisies. 

distribution des diamétres granulaires lus sur les négatifs photographiques 
grandissement 143. distribution s’accorde sensiblement dans 
les trois cas avec une formule Gauss, avec 0.10 
0.037 0.033 0.18 Ces trois valeurs s’accordent avec une droite 
dans une représentation log 1/v, selon formule: 


0.088 


sorte une vitesse infinie assure diamétre plus petit possible, 0.088 
880 

aussi tracé distribution volume volume moyen est comme 

Nous avons ensuite prélevé des échantillons des époques successives 
cours précipitation d’un lot d’émulsions. Cette série mesures sert 
début les cristaux sont petits rares, fin ils sont plus gros plus nom- 
breux. volume moyen des grains augmente mesure que volume des 


*Adresse actuelle: stage laboratoire Professeur Max Morand, Sorbonne, Paris. 
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solutions versées augmente. nombre des grains reste fixe, est pro- 
portionnel sinon, est proportionnel V/n, suppose que varie 
avec selon une loi puissance constante: 


pour diamétre: 


Effectivement, entre tours 327 tours, diamétre augmente 
selon une puissance comprise entre Cela signifie tours 
versés, nombre des grains s’accroit peu, bromure précipité 
partir moment sert surtout accroitre dimension des grains déja 
présents, forme peu nouveaux germes. 

Des dénombrements exacts grains sur des portions aliquotes des échan- 
tillons permettraient mieux comprendre mécanisme question. Les 
dénombrements sur les grains trés fins, diamétre voisin 0.04 sont 
difficiles; ceux que nous avons faire s’accordent, leurs erreurs expéri- 
mentales prés, avec les mesures 


Nous remercions Conseil National des Recherches, des Recherches 
Jeunesse (Québec), pour leur aide financiére. 


Demers, 1958. les émulsions nucléaires, principes applications (Les 
Presses Universitaires Montréal, Montréal, Qué.). 
J., DEMERS, DEMERS, 1960. Can. Phys. 38, 1482. 


RECU NOVEMBRE, 1960. 
PHYSIQUE, 
UNIVERSITE MONTREAL, 


OBSERVATIONS UNUSUAL LOW-FREQUENCY PROPAGATION 
MADE NOVEMBER, 1960 


INTRODUCTION 


The intense solar eruptions 10, 12, and November, 1960, produced 
marked effects the propagation radio waves. Some these were recorded 
series experiments which have been undertaken the Radio Physics 
Laboratory for the purpose studying the lower ionosphere high latitudes 
means propagation long waves, particular the effects polar cap 
disturbance (PCD) events. High-frequency absorption measuring equipment 
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indicates that the pronounced polar cap absorption (PCA), occurring during 
such events, usually restricted geomagnetic latitudes greater than 60°, 
and only observed south about 60° geomagnetic latitude after the com- 
mencement the geomagnetic storm, i.e. the distortion the earth’s field 
can result equatorial extension the PCA. The class flare which 
began 1325 November was accompanied major sudden ionospheric 
disturbance (SID). few minutes later, marked increase cosmic- 
ray intensity was observed, which was accompanied marked disturbance 
the region. Since the PCD began coincident with the great geo- 
magnetic storm, which began 1348 November, presumably due 
the class flare 1009 November, the observed effects the PCD 
geomagnetic latitudes 60° and below was unprecedented magnitude. 
the purpose this note give outline the effects observed 
November, 1960, and tentative interpretation them. The detailed deduc- 
tions depend much accumulated information from our studies such 
events, but not intend discuss this evidence here, since will 
published separately the near future. 


THE OBSERVED EFFECTS 


The field strength 77.15-kc and 82.05-kc transmissions from the U.S.A.F. 
transmitters Thule and Goose Bay are being monitored Churchill, and 
the field strength 70.384-kc transmissions from the Decca Navigator trans- 
mitter located Comfort Cove recorded Ottawa. addition frequency 
stable 80-kc transmissions from transmitter operated the Radio Physics 
Laboratory Ottawa are recorded Churchill. this experiment both 
phase and amplitude variations are recorded, the phase variations are measured 
with reference very stable oscillator located the receiving site. The 
various transmission paths are shown Fig. 

Records the field strength the various transmissions made 12-13 
November are shown Fig. together with similar records made 8-9 
November for comparison. The sudden beginning the PCD November 
strikingly shown the record for transmissions from Comfort Cove 
Ottawa. The onset the PCD effect followed the large SID, which began 
1323, and reached maximum 1332. The SID effect clearly one great 
magnitude, since the I.f. field strength increased value greater than the 
usual quasi-peak nighttime value (SID effects result general increase 
the field strength for these transmissions). The enhanced signal began 
decrease after 1332, the usual way, until 1410, when new and continuing 
source electron production became apparent. the absence PCD, the 
excess flare-produced ionization would have decayed, resulting return 
the I.f. field strength normal sketched the broken curve. The solar 
flare proton stream apparently produced ionized layer the base the 
ionosphere whose gradient was sufficiently sharp support the enhanced 
propagation radio waves. The SID was not seen the northern records 
because occurred too near sunrise. The enhanced field strengths the early 
phases the PCD are evident all recordings. 


611 


Fic. Map showing the various propagation paths. 


The second observable feature began 1900 November and showed 
marked decrease all field strengths. diurnal variation was 
observed sunset and sunrise and November for the Comfort 
Cove Ottawa transmission, and diurnal variations all are evident for 
high latitude transmissions for period days. 

The deduced phase variations for 80-kc transmissions between Ottawa and 
Churchill for 12-13 November are also shown Fig. together with similar 
record made November for comparison. The usual quiet day variation 
shows decreased phase lag day (lower reflection heights), increased 
phase lag night, and slow phase fading night. The SID 1325 scarcely 
discernible, since occurred too soon after sunrise produce marked 
effect. The small change detected also confused small decrease which 
began somewhat earlier, 1318. The SID effect although very small quite 
different from normal since the phase lag continued decrease from 1330 
until 1354, which time more rapid decrease began. 1515 the phase lag 
was about 0.5 cycle less than normal. There was little diurnal variation 
the phase for period days. The PCD effect November was still 
progress when second PCD began 0300 November. 
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NOTES 


not possible from our recordings relative phase lag deduce 
apparent height reflection. The change phase from night-to-day 
cycles consistent with the interpretation that the predominant mode 
propagation that once reflected sky wave. For such model the change 
phase from night-to-day corresponds change apparent reflection 
height about 11-16 km, and the decreased reflection height during the most 
intense period the disturbance, 12-13 November, corresponds 
depression the reflection height 4-6 km. thought that these 
waves are reflected from heights below during the day, perhaps from 
heights near km, and hence during the peak the disturbance reflection 
probably occurred from heights near km. 


DISCUSSION 


The beginning observable effects propagation 1354 and 1900 
November and 0300 November coincide closely with pronounced 
enhancements ground neutron monitor count rates recorded Ottawa and 
Churchill. Low-frequency propagation apparently very sensitive ionization 
density and gradient changes the the ionosphere. Weak events 
evident the nighttime propagation the northern circuits 10-11 and 
11-12 November are interpreted being PCD event following the class 
flare 1009 November. This flare originated from the same sunspot 
group those and November. Measurements made Churchill 
(Montalbetti and McEwen 1961) hydrogen emission from the night sky 
indicate the presence considerably more hydrogen than normal the 
nights 10-11 and 11-12 November. The remarkable feature lack 
any diurnal variation during the PCD suggests production electron densities 
(and electron density gradients) near km, which are adequate reflect 
waves day night independent ionospheric processes, such attachment 
free electrons form negative ions during the night. The lower edge the 
ionization produced must determined mainly the characteristics the 
bombarding particles, such the steepness the energy spectrum. 

This lack diurnal variation has already been pointed out Ortner 
al. (1960), who recorded high latitude transmissions during previous 
events. These workers addition concluded that the pertinent reflecting 
strata was located heights km, but this not consistent with our 
observations. 

When the PCD event weak, i.e. weak proton event affecting propa- 
gation high latitudes, strong event affecting propagation near the low 
Collocation propagation recordings together with recordings Mc/s 
cosmic radio noise made Mont Joli polar-directed antenna. This latter record included 
for comparison starting times the PCA with the beginning the effects observed our 
Comfort Cove Ottawa transmission. The solid curves show the variations observed 
12-13 November, and the broken curves show variations observed normal day for com- 
parison (which for the signal amplitude records the 8-9 November, and for the phase record 
the 10-11 November). The short vertical lines represent the times when the sun’s zenith 
angle 102° and 90° 50’ (ground sunrise and sunset) the mid-point the path. The vertical 
arrows with captions SC, SECR, and represent the times the beginning the 


geomagnetic storm, the beginning sudden enhancements cosmic rays observed Ottawa 
and Churchill, and the beginning the visible solar flare. 
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latitude extent the disturbance, the diurnal variation field strength 
transmissions opposite that normally observed, i.e. the daytime field 
strengths are greater than normal, and the nighttime field strengths are less 
than normal. The transition from night-to-day occurs quite suddenly during 
twilight conditions the path mid-point. marked feature 
pattern the absence fading, indicating quasi equilibrium the ionization 
produced the base the ionosphere. When the proton event strong the 
diurnal variation disappears entirely, and the observed field strength about 
that normal day. This accounts for the apparent differences between the 
observed behavior 16-kc propagation from Rugby received Upsala 
(during the February 1956 event (Belrose a/. 1956)), and received 
Kiruna (during the May and June 1959 events (Ortner al. 1960)). the 
disturbance weakens, the nighttime propagation which predominantly 
affected, i.e. the excess ionization due the slowly dying disturbance lies 
between those heights where the waves are normally reflected during the day 
and during the night (between 65-85 km). This conclusion substantiated 
our measurement phase. 


BELROsE, H., and WEEKEs, 1956. Atmospheric and Terrest. 


Phys. 281. 
and McEwen, 1961. Can. Phys. 39. This issue. 
J., EGELAND, A., and 1960. IRE Trans. AP-8, 621. 
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THE IONOSPHERIC DISTURBANCE NOVEMBER AND 16, 1960 


The purpose this note summarize the absorption effects recorded 
ionospheric sounding stations Canada during the disturbed period between 
November and 16, 1960. The locations the stations and the characteristics 
the ionosondes are tabulated Table Ottawa, St. John’s, and Winnipeg 
are considered middle latitude stations; Churchill located the 
auroral zone, and Resolute Bay situated the polar cap region. The routine 
sounding program consists frequency sweeps through the h.f. spectrum 
15-minute intervals. 


TABLE 
Geomagnetic Peak pulse 
Station Latitude Longitude latitude power (kw) 
Churchill 94.2° 68.0° 
Winnipeg 97.4° 59.8° 
St. John’s 52.7° 58.0° 
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The best parameter for the qualitative measurement absorption that can 
obtained from ionograms the lowest radio wave frequency which 
echo recorded The fmm increases with increasing absorption. When 
the exceeds the highest radio wave frequency, which the ionosphere can 
reflect that time, echoes will observed, condition generally known 

since relative the over-all sensitivity the sounding system. 

good indication the presence abnormal absorption can obtained 
comparing the observed value with the value normally 
observed other days the same time. The parameter will used 
indicate the difference between the observed and its monthly median 
value for the same time. Fig. shown function universal 
time for all five stations, for the period November Times sunrise 
and sunset height and times occurrence sporadic are 
also shown. 

may noticed from Fig. that the durations black-out the polar 
and auroral zone stations (Resolute and Churchill) are quite similar. The two 
black-outs are separated period partial recovery lasting for few hours, 
the nights November and 15. The other three medium latitude stations 
all show similar characteristics but the period partial recovery begins one 
day earlier the night November The recovery period charac- 
terized strong absorption and sporadic which are typical the auroral 
ionosphere. The durations the main black-outs Winnipeg were some- 
what longer than Ottawa and St. John’s. 

The main features the black-outs are summarized Table II. 


TABLE 
Black-out starting Black-out starting 
Nov. Nov. 
Approx. time Approx. time Conditions First 
start Duration start Duration prior normal day 
black-out black-out the after the 
Station hours hours black-outs black-out 
Resolute Bay 1530 0430 Occasional hig 
25/11/60 
Churchill 1515 0530 Occasional high 
19/11/60 
Winnipeg 1800 1200 Quiet 18/11/60 
St. John’s 1800 1545 Quiet 18/11/60 
Ottawa 1815 1330 Quiet 17/11/60 


Acknowledgment given the Department Transport for operation 
the ionospheric sounding stations. 
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NOTES 


HYDROGEN EMISSIONS DURING THE PERIOD NOVEMBER 9-16, 1960 


Churchill, Manitoba (geomagnetic latitude photometer 
operated continuous basis monitor the emissions from the night 
sky. consists simply narrow band interference filter oscillated scan 
wavelength interval about 200 some seconds and photomultiplier 
tube detector. The photometer programmed view five positions 
sequence 40° steps, starting 10° above the northern horizon. complete 
survey along the meridional line made minutes. 

Fig. the intensity the emissions, grouped into intensity intervals 
given the legend for the figure, plotted function universal time 
and for the five directions. Figure also shows the maximum variation 
the component the magnetic field Churchill for each hourly interval. 
1007, and Any variation greater than 4007 indexed 

The strongest hydrogen emissions were observed the night November 
12/13 when the intensity reached some 600 According Chamberlain’s 
(1954) analysis this corresponds flux the order protons/cm? sec. 
Within the time resolution each meridional survey systematic motion 
evident. Strong bursts emissions appeared over the entire meridional 
line about the same time. These bursts follow strong radio-wave absorption 
peaks measured the Churchill riometer aimed towards the zenith. 
Table the times strong auroral-type absorption peaks and the times 
strong emissions are given. 


TABLE 


Time peak Time peak Max. value Max. value 


0700 Nov. 0740 270 >12 
0800 350 
1015 1030 650 >12 
1145 1210 330 
0045 Nov. sky) 
0115 
0530 0555 
0615 0625 


Whether there were strong emissions between 0600 hours and 0618 hours 
November cannot determined aurora occurred such intensity 
saturate the photometer. The time delay between the radio-wave absorption 
and the light emissions the order half hour. During the night 


*Defence Research Telecommunications Establishment, Ottawa, Ont. 
Research Northern Laboratory, Churchill, Man. 
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VERY DISTURBED 


Fic. Intensity emissions and magnetic disturbances Churchill, Manitoba, for 
the period November 16, 1960. 


NOTES 619 


November 12/13 strong rayed auroral forms and, more important, strong 
type aurorae were observed. Observations Churchill indicate that strong 
hydrogen emissions occur during type auroral displays. 

For the two nights prior the strong disturbance, hydrogen emissions 
were observed practically all night although they were not intense. 
November 10/11 the peak the emissions appeared move from the northern 
the southern horizon while November 11/12 rather weak intensities 
were measured covering the entire sky for many hours. would expected 
that the proton fluxes causing emissions are responsible, part, for the 
radio absorption phenomena, then absorption should have been observed 
November 10/11. This does not appear the case, though Belrose has 
reported evidence for enhanced low-level ionization this time. 

the basis absorption and magnetic observations, there were apparently 
two large events, namely one during November 12/13 and one during Novem- 
ber 15/16. Large proton fluxes and unusual aurorae were observed with the 
first event. unfortunate that the nights November 14/15 and November 
15/16 were completely overcast. Thus unknown whether large proton 
fluxes were associated with the second event. Observations from other stations 
whether type aurorae occurred would valuable. 

The 600 measured November 12/13 must taken unusual. The 
peak intensities observed with the photometer operation since December 
10, 1959, have been approximately 150 February 10/11, 1958, when 
similar geomagnetic storm occurred, intensities reached values 850 
(1959) measured with scanning spectrometer. 


BELROsE, and Ross, 1961. Can. Phys. 39. This issue. 
CHAMBERLAIN, 1954. Astrophys. 120, 360. 
1959. Atmospheric and Terrest. Phys. 14, 200. 
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GEOMAGNETIC VARIATIONS BETWEEN NOVEMBER AND 
NOVEMBER 16, 1960 


NIBLETT 


Exceptionally large geomagnetic disturbances have been reported from 
observatories many different parts the world during this interval. 
Canadian observatories the disturbances between the hours U.T. and 
U.T. November were comparable with the largest that have been 
recorded this country. 
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Figure shows mean hourly values horizontal force, Victoria, 
Agincourt, and Meanook observatories during the five Greenwich days. 
Baker Lake and Resolute the mean hourly values the principal horizontal 
component for Baker Lake and for Resolute) have been plotted, since 
not recorded these locations. The lowest plot Fig. shows the hourly 
range Meanook observatory. The times two sudden storm com- 
mencements, 13:48 November and 13:040n November 15, are indicated 
arrows. The locations the five observatories are given Table 
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Fic. Horizontal force variations Canadian magnetic observatories. 
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TABLE 
Geographic Geomagnetic 
Observatory Latitude Longitude Latitude Longitude 
Victoria 123.4° 54.3°N. 292.7° 
Agincourt 43.8 79.3 55.2 346.9 
Meanook 54.6 113.3 61.9 300.7 
Baker Lake 64.3 96.0 73.9 314.8 


The stations south the auroral zone, Victoria, Agincourt, and Meanook, 
show similar pattern disturbance throughout the five days. Two magnetic 
storms are evident, separated more than day relatively little activity. 
Significant features the first storm are decrease from about 16:30 
19:30 November 12, enhancement from then till the 13th, 
and extremely large depression between and the same day. 
Reasonably quiet conditions were not restored till approximately 02:30 
the 14th. The very large rapid fluctuations typical magnetic storms near 
the auroral zones developed hours after the SSC 13:48 and persisted 
till about 21:00 the 13th. 

the second and less intense storm the horizontal force went through 
moderate depression immediately after the SSC, strong enhancement which 
reached its peak about hours the 16th, and final moderate depres- 
sion. Moderate rapid fluctuations set after the SSC and large storm 
fluctuations occurred between 21:00 the 15th and 09:00 the 16th. 

All the variations described above were larger Meanook than Victoria 
Agincourt. This suggests that the succession events was controlled 
these locations auroral zone currents, not equatorial ring current. 
Classical patterns which are often obvious magnetograms from 
low and middle latitudes are seldom evident Canada during individual 
storm. the other hand the horizontal force the three observatories shows 
tendency depressed during the local a.m. hours and enhanced during 
the p.m. hours. This consistent with the usual behavior disturbance 
daily variation near the auroral zone. 

The two observatories north the auroral zone show intense but less 
regular disturbance patterns. The SSC the 15th almost obscured 
large background noise both stations. Baker Lake shows scarcely any 
contrast between the second storm and the irregular fluctuations preceding 
it. Systematic trends during the a.m. and p.m. hours are not obvious there. 
Resolute both storms are characterized enhanced during the p.m. hours 
and depressed during the a.m. hours. 

Some data the two sudden storm commencements and indices are 
given Table II. third SSC reported the IUGG Association Geo- 
magnetism and Aeronomy have occurred 10:21 November 13. Con- 
ditions the Canadian observatories were much too disturbed this time 
for identification possible. 
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TABLE 
Activity 
Amplitude 
Periods 
Resolute Nov.12 13:46 


Norte: Small differences in the times recorded for SSC’s at different observatories are probably due to clock 
errors on the magnetograms. The K indices for Resolute are not available. 2Kp is the daily sum of the eight plane- 
tary K indices. 2K is the daily sum of the eight indices at the specified observatory. “Periods of Kmax” refers 
the to three hourly intervals during which the K index was a maximum. 


RECEIVED FEBRUARY 1961. 

DoMINION OBSERVATORY, 

DEPARTMENT MINES AND TECHNICAL SURVEYS, 
CANADA. 


OBSERVATIONS AURORAL IONIZATION DURING THE PERIOD 
NOVEMBER 11-16, 1960 


The present note describes observations auroral ionization made the 
Defence Research Telecommunications Establishment during the period 
from November November 16, 1960, inclusive. Measurements both 
auroral luminosity and the amplitude echoes obtained radar systems 
operating ultra high frequencies were made. These measurements were 
supplemented observations the magnetic disturbance. The data are 
displayed graphically the accompanying diagrams. 
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Fic. Radar, photometric, and magnetic observations, November 13, 1960. 


The auroral radar observations were made two similar systems operating 
488 Mc/s and 944 Mc/s, each having peak power kw, pulse length 
500 usec, and beam width 2.5°. The transmitting and receiving sites are 
situated just west Ottawa. The antennas illuminate area the 100-km 
height level centered roughly over Moosonee, Ontario (geomagnetic latitude 
62.2° N.). from this region that echoes are obtained. 

The histograms give values the peak amplitude echoes for intervals 
1/100th day. Echo amplitude plotted 5-db steps above noise levels 
dbm and —123 dbm for the 488 Mc/s system and 944 Mc/s system 
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Fic. Radar, photometric, and magnetic observations, November 16, 1960. 


respectively. Comparison the relative amplitude echoes the two 
frequencies during periods when both systems were operating simultaneously 
suggests that the 488 Mc/s system would probably have been saturated 
large part the time from late November until noon November 13. 

two-channel recording photometer located the receiving site having 
beam width, bearing, and elevation similar that the radar receiving 
system. This photometer arranged look alternately through 5577 
and filter. The results are plotted below the radar observations 
the diagrams. The ordinate scale gives the peak value the intensity for 
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periods 1/100th day increasing decade steps. scale value unity 
corresponds intensity less than about kilorayleigh; scale value 
indicates intensities greater than 100 kilorayleighs. The photometric 
records are interrupted daylight and the presence cloud cover. 

noticed that the relative intensity the two wavelengths does not 
always remain constant. This made evident inspection the results 
the two wavelengths for the period just before noon November 13. 
significant that during the particular times when the intensity the 3914 
band exceeded that the line, Montalbetti and McEwen (1961) 
report unusually strong emissions Churchill. 

The peak value the change the magnetic X-component Moosonee, 
Ontario, measured three-component flux gate magnetometer, also 
given the figures. Portions this record are missing because failure 
the recording meter. Deviations from quiet day values greater than 1000 
gammas were frequently obtained during the period question. 

The period from November November 16, 1960, showed unusually 
high auroral activity particularly during the nights November 12/13 and 
November 15/16. Although there direct correlation between the radar 
and photometric observations aurora, strong features appear both records 
associated with large changes the magnetic field. 


and 1961. Can. Phys. 39. This issue. 


RECEIVED FEBRUARY 1961. 
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DEFENCE RESEARCH BOARD, 
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AURORAL RADAR OBSERVATIONS MC/S DURING THE PERIOD 
THE NOVEMBER, 1960, SOLAR EVENT* 


Simultaneous records from three auroral radars similar characteristics 
have been analyzed for the period November November, 1960. The 
radars are located Churchill, Saskatoon, and Ottawa, for which the geo- 
magnetic latitudes are 68° N., 60.5° N., and 57° The radars operate near 
48.5 Mc/s with pulse 1-kilowatt peak power. The data plotted 
are the occurrences echoes half-hourly intervals, starting the hour 
and half-hour. The solid bars Fig. indicate auroral echoes obtained 
Churchill, Saskatoon, and Ottawa. The histograms are the half-hourly meteor 
echo counts the Ottawa auroral radar. For reference other solar and 
geophysical events, the times the two major sudden commencements (SC) 
magnetic storms November and November are indicated the 
arrows. 

*Issued N.R.C. No. 6226. 
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NOV, 


$ 


Fic. occurrences Mc/s auroral radars. The solid bars are the occurrences 
auroral echoes Churchill (C), Saskatoon (S), and Ottawa (O). The histograms are the 


meteor echo rates the Ottawa auroral radar. indicates the time sudden commence- 
ment. 
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The primary event solar flare 13:25 U.T., November, followed 
cosmic-ray increase, and followed turn 13:50 U.T. sudden com- 
mencement. Within minutes the SC, echoes were observed for 7-minute 
period Ottawa. Strong signal absorption, indicated abrupt drop 
the meteor rate the Ottawa auroral radar, began the period 14:30 U.T. 
The meteor count reached minimum from 19:00 20:00 U.T. when 
meteors were detected. Thereafter there was slow return meteor echoes, 
the mean rate during the following hours (13 November) being about one- 
third normal. the other two radars, the meteor rates were also reduced 
fraction the normal rate for period approximately hours. Since 
the meteoric reflections occur heights around km, the meteor effect 
interpreted caused widespread low-level absorption, polar cap 
absorption (PCA) event. 

the hours following the November sudden commencement, the 
behavior the auroral echo occurrence the three stations noteworthy. 
Echoes started Ottawa 21:30 U.T., November, and were nearly con- 
tinuous for hours. The height reflection auroral echoes km. 
Activity the Saskatoon radar during the same time was present but weaker 
and more intermittent, while Churchill there were echoes observed 
the entire 24-hour period. This geographical pattern occurrence could 
attributed real variation the ionization distribution. However, when 
considered the light observations showing large values radio-frequency 
absorption, and reports intense visual and magnetic activity high lati- 
tudes, appears more likely that intense echoes 48.5 Mc/s would have been 
present all stations not for strong absorption low heights. 

November, another sequence solar flare, cosmic-ray increase, and 
sudden commencement occurred. The 13:05 U.T. was followed 
minutes echoes, this time only the Saskatoon radar. During the hours 
following the SC, the meteor rates fell low values, but recovered rapidly 
normal Ottawa with the coming darkness. feature peculiar the 
November sequence the short period nighttime absorption centered 
around 08:00 09:00 U.T. November. 

The combined results three widely separated v.h.f. auroral radars show 
interesting patterns activity. Short bursts echoes are sometimes asso- 
ciated time with sudden commencements. These echoes are not intense 
which may account for the fact that only the more favorably situated radars 
detect them. The evidence cited for this large PCA event suggests that 
such times the strong low-level absorption significantly alters the detection 
capabilities v.h.f. auroral radar, and may distort the apparent geographical 
distribution the reflecting ionization. 


RECEIVED FEBRUARY 1961. 

DIVISION RADIO AND ELECTRICAL ENGINEERING, 
NATIONAL RESEARCH COUNCIL, 

CANADA. 
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PROPAGATION METEOR BURST SIGNALS DURING THE 
POLAR DISTURBANCE NOVEMBER 12-16, 1960 


MAYNARD 


meteor burst communications system, short signal impulses are observed 
receiving terminal each time suitable ionized trail established the 
lower ionosphere. Most signal bursts arise from meteors originating the 
vicinity 100 km. 

D-region absorption present over large area near the center the 
transmission path, then meteor burst signals will absorbed, provided the 
absorption takes place below the meteoric region. 

Since meteoric signals are random amplitude and occurrence, the effect 
absorption the signals must considered statistically. This may 
done measuring the time that the signal level remains above certain gate 
level decision level. This time, expressed fraction the sampled period, 
called duty cycle. The dependence duty cycle decision level has been 
fairly well established. first approximation, found that for 10-db 
change decision level the receiver, the duty cycle changes approxi- 
mately can shown that 10-db change decision level 
essentially equivalent 10-db change transmitter power 10-db 
change transmission loss due absorption. rough estimate changes 
transmission loss due absorption can made observing deviations 
from the normal diurnal variation signal duty cycle. This note will describe 
the behavior meteor burst propagation link during the large disturbance 
November 12-16, 1960. 

Figure shows the Mc/s absorption mid-path meteor burst 
communications circuit between Goose Bay, Labrador, and Ottawa, Ontario, 
measured Mc/s riometer. Times shown these figures are U.T. 

Figure shows plot 104.1 Mc/s meteoric signal duty cycle over this 
path. The shaded region this graph represents the normal limits the 
variation duty cycle due the diurnal variation meteoric activity and 
the random nature the meteoric signal. Signal duty cycle seen drop 
below its normal range for three intervals during the four days plotted. These 
periods signal black-out coincide fairly well with peaks absorption 
measured the riometer mid-path the link. These three periods are 
the only times test period over months that 104.1 Mc/s meteor 
burst signals have been significantly affected absorption. Gaps duty 
cycle plot this figure and Fig. are due mainly cluttering the meteoric 
signals auroral-type signals. 

The 41.3 Mc/s signal duty cycle shown Fig. affected greater 
extent the absorption. Significant changes duty cycle occur early 
1400 U.T. November 12, about hour after the solar flare 1325 U.T. 

During the period from 18:00 20:30, November 12, the 104.1 Mc/s duty 
cycle deviated factor greater than 1000 from normal, causing 
signal black-out. 


Can. J. Phys. Vol. 39 (1961) 


ABSORPTION (db) 


6 12 18 6 2 6 12 06 12 


NOV. I2 1960 NOV. 13 1960 NOV. 14 1960 NOV. 15 


107! 


073 


METEORIC SIGNAL DUTY CYCLE 


1075 
7 6 12 6 12 18 6 12 18 6 12 18 
NOV. 12 1960 NOV. 13 1960 NOV. 14 1960 NOV. 15 1960 


METEORIC SIGNAL DUTY CYCLE 


18 6 18 


6 18 6 12 1s 12 
NOV. I3 1960 NOV. 1960 NOV. 15 1960 


Fic. 104.1 Mc/s duty cycle. 
Fic. 41.3 Mc/s duty cycle. 


6 12 
NOV. [2 I960 


NOTES 629 
12 
NO 
RECORDS 
INTERFERENCE 
i 


630 CANADIAN JOURNAL PHYSICS. VOL. 39, 1961 


This sudden decrease duty cycle coincides with the abrupt rise absorp- 
tion 1800 shown Fig. This change duty cycle could interpreted 
oblique signal absorption excess 104.1 Mc/s. 

approximate relationship between oblique signal absorption and vertical 
signal absorption might assumed 


2 
oblique frequency used over the transmission path, and the angle 
incidence the layer. 

one interprets the observed signal effects terms the absorption one 
might expect observe the mid-path riometer, appears that the signal 
effects are adequate assume the absorption was located mainly below the 
meteoric region 100 km. 


RECEIVED FEBRUARY 1961. 

COMMUNICATIONS LABORATORY, 

DEFENCE RESEARCH TELECOMMUNICATIONS ESTABLISHMENT, 
DEFENCE RESEARCH BOARD, 

SHIRLEY Bay, OTTAWA, ONTARIO. 


IONOSPHERIC ABSORPTION NOVEMBER 12, 1960 


active region appeared the sun early November, 1960. During its 
passage across the disk gave rise number large and outstanding flares, 
which turn produced series noteworthy terrestrial disturbances. Flares 
this region emitted relativistic particles least three occasions and 
caused number large geomagnetic and ionospheric storms. Table are 
listed the principal flares and the sudden commencement geomagnetic storms 
that, the basis the optical and radio importance the flares, are thought 
result from them. least one case, some ambiguity exists concerning 
the causative flare. 

Throughout the disturbed period, the Radio Physics Laboratory operated 
Mc/s cosmic noise receivers, riometers they have been called (Little 
and Leinbach 1959), Resolute Bay, Coral Harbour, Churchill, Cape Jones, 
Val d’Or, and Ottawa. Churchill and Ottawa, additional riometers 
Mc/s were also operation. Figure shows recording made one 
the Ottawa Mc/s riometers during the period 1200-1800 U.T. November 
12. The significant features this record are the rapid solar noise bursts start- 
ing 1328 and 1344, the sudden ionospheric disturbance beginning 1326 
and continuing until least 1430, the intense and spiky solar noise during 
the period 1420 least 1800, the intense (greater than db) polar cap 
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type absorption evident after 1800, but for which evidences can seen 
this figure early 1540. (It should noted that the receiving antenna 
was not directed the sun that absolute levels solar noise are not avail- 
able from the riometer record; nevertheless the character the solar noise 
the interval can seen Fig. 1.) 

Polar cap absorption (PCA), though usually confined geomagnetic 
latitudes greater than about 60°, has been known occur lower latitudes 
after some large perturbation the earth’s magnetic field. November 12, 
the sudden commencement magnetic storm began very soon after PCA 
producing flare, that this type absorption was observable that day 
all the riometers. Moreover, since the magnetic field was disturbed for 
number days, polar cap absorption was observable all stations least 
until November 18, and some stations until November 22. more detailed 
description these events will published forthcoming paper; the 
present contribution discusses some the distinctive features the data. 

Some features the early stages the PCA November are portrayed 
latitude-time profile Fig. The respective curves give the times the 
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profile showing start polar cap absorption Mc/s for Novem- 
ber 12, 1960. 


following features: the start the PCA, curve sharp change 
the rate increase absorption, curve and the attainment the 4-, 6-, 
and 8-db absorption levels. each case the horizontal bars indicate the un- 
certainty the timing that point, while the dots curve give the most 
probable time determined close inspection all available data the 
respective stations. the start the flare, the ionosphere was illuminated 
the sun all stations except Resolute: the apparent delay the start 
the PCA there should not taken evidence for delay the arrival 
the solar protons. Furthermore, the Coral Harbour records the SID, 
occurred, cannot distinguished from the start the PCA. Since the ground 
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level cosmic-ray increase did not start before 1345 U.T. (Rose, unpublished 
data), seems reasonable assume that the PCA also did not start earlier 
than this, even though increased absorption evident the records early 
1326 U.T. the southern stations the late start for the PCA must 
taken somewhat uncertain because the presence solar noise all the 
recordings. attempt ascertain the importance the solar noise the 
various parts the riometer records, the solar emissions were evaluated from 
independent recordings Mc/s made with lobe-switched interferometer. 
While these data cannot used remove the uncertainty the time the 
start the polar cap absorption, they permit better appraisal made 
the location time the higher absorption contours the respective 
riometer records. Furthermore, such prominent features the sharp change 
the rate increase absorption, curve remain unaffected the solar 
noise. this basis the shape the absorption contours shown Fig. must 
considered real, leading the conclusion that the build-up absorption 
the southern stations proceeds more slowly than the more northerly 
stations. 

approximately 1900 U.T. November 12, the ground level cosmic-ray 
monitors observed second pronounced increase counting rate (Rose, 
unpublished data). While the Mc/s riometers showed change this 
time—they were already beyond their dynamic range—the Mc/s riometers 
Ottawa and Churchill both showed sharp increase absorption. 

Yet another absorption feature observable some the records 
November 12. Coincident with the magnetic storm sudden commencement 
1348, the Cape Jones record and, lesser extent, the Churchill record showed 
the short duration absorption feature that has been identified elsewhere 
ray induced and feature the sudden commencement the auroral zone 
(Brown al. 1961). Similarly, November approximately 1304 U.T. 
the same feature evident the Churchill and Cape Jones records. 

The start the PCA November occurred before local sunrise the 
ionosphere each the observing stations: the steep increase absorption 
took place the ionosphere became illuminated the sun. The actual starting 
time the weak nighttime PCA somewhat uncertain because auroral 
absorption during the early morning hours, but there evidence from the 
various stations for starting time 0720 U.T. with uncertainty +30 
minutes. not possible identify any systematic variation starting 
time with latitude for this event. 

The PCA November did not exceed whereas the two earlier events 
showed absorption greater than db. 


The authors are indebted the various people who assisted with the 
operation the riometers and with the data analysis. They have also made 
use unpublished data from the Fraunhofer Institute, Freiburg, Germany, 
the High Altitude Observatory, Boulder, Colorado, and the National Bureau 
Standards, Boulder, Colorado. 
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Brown, R., Hartz, R., LANDMARK, B., LEINBACH, H., and 1961. 
Geophys. Research. press. 
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THE SOLAR NOISE BURST NOVEMBER 12, 1960* 


The intense burst solar radio noise which accompanied the now well- 
known flare November 12, 1960, has been observed the discrete fre- 
quencies 2800 Mc/s and Mc/s the Radio Observatory the Radio 
and Electrical Engineering Division the National Research Council. These 
observations will presented and the essential information noted. Figure 
shows reconstruction the high-frequency burst while Fig. shows portion 
the low-frequency record. 

2800 Mc/s event.—When the solar patrol observations were begun 40™ 
U.T., the solar flux was found increasing slightly, and 13" 20" U.T. 
increase units flux had been recorded excess the daily level 
147 flux units (where unit flux equal watts per square meter 
per cycle per second). this time, 13" 20™, the increase became more rapid and 
showed noticeable departure from the residual receiver noise; 13" 23™ the 
flux had reached flux units, 13" the flux had reached 100 units, and 
27™ level 3000 units. The intense part the burst shows number 
oscillations, with peak value flux 5500 units 13" The burst 
declined rapidly from this peak value level 300 units 14" 06". This 
was followed enhanced post burst increase and the occurrence 
number superimposed relatively smaller bursts. 00™ the quiet sun 
level was reached and further activity was recorded during the remaining 
two hours observations. 

Mc/s event.—A portion the Mc/s film record shown Fig. (note 
that the times recorded the film are E.S.T. and are minute fast). This 
record has the output the 2800 Mc/s radiometer simultaneously recorded 
single curve which moves upwards for increasing signal strength. The 
sloping curve across the upper track represents portion the 2800 Mc/s 
record the burst level increases from flux units 160 flux units. Fig. 
these levels are shown broken lines. The time, 13" which the 
85-flux level reached has been taken the high-frequency burst commence- 
ment. The Mc/s burst started abruptly and relation 
the somewhat arbitrarily chosen starting time for the higher frequency burst, 


*Issued N.R.C. No. 6234. 
Can. J. Phys. Vol. 39 (1961) 


636 CANADIAN JOURNAL PHYSICS. VOL. 39, 1961 


FLUX 
5600 
NOVEMBER 12, 1960 
RECONSTRUCTION 
UT. 
1330 1430 
FLUX 


180 


1500 1600 1700 
Fic. 2800 Mc/s solar burst. 


delay time 189 seconds obtained. The Mc/s storm continued until 

The two faint bands noise approximately 13" 23™ and 24™ 
are probably due other weaker type III bursts associated with the solar 
eruption. The other broad, faint bands noise 13" 26™ are due inter- 
ference. Other burst activity occurred before the principal event, from 30™ 
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LETTERS THE EDITOR 


Under this heading brief reports important discoveries physics may published. These 
reports should not exceed 600 words and, for any issue, should submitted not later than six weeks 
previous the first day the month proof will sent the authors. 


Frequency Measurement Standard Frequency 


Measurements are made Ottawa, Canada, using N.R.C. caesium-beam frequency resonator 
reference standard (with assumed frequency 192 631 770 c.p.s.). Frequency deviations 
from nominal are quoted parts per negative sign indicates that the frequency below 


GBR, 16 ke/s 


Date, 
January MSF, 8}-hour 24-hour WWVB, 
1961 60 kc/s average* averaget 60 kc/s 

1 N.M. — 133 — 136 N.M. 

2 —129 —134 — 133 N.M. 

3 —148 — 138 —139 —155 

4 —141 —140 —140 —154 

5 —143 —134 —139 N.M. 

7 —137 —138 —140 N.M. 

9 —150 — 138 —140 —149 

10 —157 —141 —144 —154 

1l —146 —140 —142 — 152 

12 —154 —141 —141 — 146 

13 — 153 —145 —148 —154 

14 —148 — 145 —147 N.M. 

15 — 146 — 146 —144 N.M. 

16 N.M. — 146 —147 —154 

17 N.M. —145 —147 —159 

19 —122 — 146 —147 — 142 

20 N.M. —144 — 147 —149 

N.M. N.M N.M 
N.M. N.M N.M 
23 — 167 N.M. N.M — 156 
24 —155 —147 — 150 —155 
25 —155 —150 — 149 — 155 
26 — 153 —144 —149 —152 
27 —151 — 142 —147 —151 
28 —149 —149 —150 N.M 
29 — 152 —149 — 152 N.M 
30 — 158 — 153 — 153 — 155 
31 — 162 — 153 —153 —154 
Average —150 —143 —145 —153 

Midmonthly 
mean of WWV —150 


Note: N.M. no measurement. 
*Time of observations: 22.45 to 07.15 hours U.T. 
tTime of observations: 15 to 15 hours U.T. 


RECEIVED FEBRUARY 13, 1961. S. N. KALRA 
Division OF APPLIED PHYSICS, 

NATIONAL RESEARCH COUNCIL, 

Ottawa, CANADA. 


1[ssued as N.R.C. No. 6220. 
2Cf, Kalra, S. N. 1959. Can. J. Phys. 37, 1328. 
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ANNOUNCEMENTS 


The New Unit Atomic Mass* 


the 10th General Assembly the International Union Pure and Applied Physics, 
held Ottawa September, 1960, the following resolution was adopted unanimously: 
“The 10th General Assembly the I.U.P.A.P. recommends the adoption the 
exact number the relative nuclidic mass the carbon isotope mass number 
12. This action will effect unification the physical scale relative nuclidic masses 
and the chemical scale atomic 


This resolution was recommended the Study Committee (of Nuclidic 
Masses (President, Mattauch), the Commission (of Symbols, Units, and 
Nomenclature, and the International Union Pure and Applied Chemistry. 

The new unit mass designated the symbol thus, Readers are 

the instruction the General Assembly the Commission Nuclidic 
Masses has arranged, public service, distribute reprints new mass table (Everling, 
Mattauch, and Wapstra, Nuclear Phys. 18, 529 (1960)) which nuclidic masses are 
expressed terms both the new unit and the old unit atomic mass. These reprints may 
obtained from the undersigned cents per copy copies for $4.00. 


H. E. DucKwortH 
SECRETARY 
COMMISSION ON NUCLIDIC MASSES 


RECEIVED MARCH 2, 1961. 
McMASTER UNIVERSITY, 
HAMILTON, ONTARIO. 


*This work was supported in part by the U.S. Air Force through the Air Force Office of Scientific Research 
of the Air Research and Development Command, and by the National Research Council of Canada. 


The New Standard for the Metre! 


event took place October, 1960, the 11th International Conference Weights 
and Measures which great interest every scientist concerned with the measurement 
length. that date the physical prototype for the International Metre, platinum iridium 
bar bearing two fine lines, which, since 1889, has preserved the world’s basic unit length, 
was relegated the position honored scientific museum piece the adoption the 
following resolution base henceforth the International Metre wavelength light: 


Métre est longueur égale 650 763, longueurs d’onde dans vide 
krypton 86. 

Définition Métre vigueur depuis 1889, fondée sur Prototype Inter- 
national platine iridié, est abrogée. 

Prototype International Métre sanctionné par Premitre Conférence 
Générale des Poids Mesures 1889 sera conservé Bureau International des 
Poids Mesures dans les mémes conditions que celles qui ont été fixées 1889.” 


This new definition the metre was drawn intentionally broad general terms that 
would not complicated practical instructions for its realization which will changed 
from time time take account new spectroscopic, interferometric, electronic, and other 
techniques that will permit better realization the ideal definition. The International 
Committee Weights and Measures was invited second resolution the 11th Con- 
ference issue such instructions from time time and its October meeting issued its 
first, which were follows: 

“Conformément paragraphe Résolution adoptée par Onzitme Con- 
férence Générale des Poids Mesures (octobre 1960), Comité International des 
Poids Mesures recommande que radiation krypton adoptée comme étalon 
fondamental longueur soit réalisée moyen d’une lampe décharge cathode 
chaude contenant krypton d’une pureté non inférieure pour cent, 
quantité suffisante pour assurer présence krypton solide température 
64° cette lampe étant munie d’un capillaire ayant les caractéristiques suivantes: 
diamétre intérieur millimétres, épaisseur des parois millimétre environ. 

estime que longueur d’onde radiation émise par colonne positive 
transition entre les niveaux non perturbés, lorsque les conditions suivantes sont 


satisfaites: 
1Issued as N.R.C. No. 6239. 
Can. J. Phys. Vol. 39 (1961) 
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capillaire est observé bout que les rayons lumineux utilisés chemi- 
nent cathodique vers anodique; 

partie inférieure lampe, compris capillaire, est immergée dans 
bain réfrigérant maintenu température point triple degré prés; 

densité courant dans capillaire est 0,3+0,1 ampére par centimétre 


Physicists will probably interested the way which the transition from the prototype 
standard the wavelength standard was effected. 1927 the 7th International Conference 
Weights and Measures authorized relationship between the wavelength the cadmium 
red line and the International Metre which made the wavelength normal air 6438.4696 
metres. this value all spectroscopic wavelengths have since been based well all inter- 
ferometric determinations length for science and industry. The value chosen was the average 
number experiments conducted 1927 relate this wavelength the International 
Metre. should pointed out that this action risked the possibility that, some later date 
when the metre was defined officially wavelength basis, there might arise the necessity 
accepting either discontinuity the length the metre revision all spectroscopic 
tables and interferometric determinations length. The 1927 relationship between the metre 
and the cadmium line was defined for normal air instead vacuum. This introduced measure 
uncertainty succeeding years the constitution normal air. When the consultive 
committee for the definition the metre met 1953 consider the procedure used 
the transition from the physical standard the metre wavelength standard recognized 
the necessity that the new definition terms vacuum wavelength. therefore decided 
that the 1927 wavelength air for the red line should reduced vacuum means the 
formula—often called the Edlén formula—for the dispersion normal air adopted Rome 
1952 the Mixed Commission for Spectroscopy. Having thus established the vacuum 
wavelength for the red line cadmium was possible relate 
the vacuum wavelength any other spectral line deemed potential candidate for use 
the new base for the metre. The spectral lines three elements particular were intensely 
studied for their general suitability and particular for their sensitivity changes wave- 
length and line symmetry result perturbations imposed the atoms such factors 
pressure, temperature, electric current, etc. The elements were isotopes cadmium, mercury, 
and krypton. Finally, krypton was chosen the most suitable. Details the research work 
which this decision was based and the discussion them can found the Procés Ver- 
baux the International Committee and the Consultive Committee for the Definition 
the Metre for the years 1953 1960 well elsewhere the scientific literature. The 
International Committee accepted the recommendation the Consultive Committee and 
presented the resolution cited above the 11th General Conference. 

The International Committee and the General Conference could reassured that 
duality the metre had been introduced the action 1927 defining the cadmium red 
line terms the international prototype metre some recent Canadian work reported 
the Conference the time its meeting. Dr. Baird the Division Applied 
Physics, National Research Council, made observations four 1-metre bars compared directly, 
recently, and different times, with the international prototype Sévres and produced 
evidence show that the new definition the International Metre preserved the continuity 
this unit which the order accuracy and precision with which the International 
Metre has been reproduced since 1889. 

perhaps worth noting some facts about the International Conference Weights and 
Measures which unfortunately not too familiar the average physicist. The Conference 
probably the oldest international scientific organization the world. has operated 
since 1875 under the authority treaty called the Convention Métre and generally 
responsible for the whole field fundamental physical standards measurement and the 
orderly development our international measurement system. The Conference itself meets 
least once six years. International Committee Weights and Measures, elected 
the Conference, consisting members chosen for their individual qualifications the 
field fundamental measurement, responsible for the actual carrying out the work under 
the authority the Conference and for the direction the laboratories the International 
Bureau Weights and Measures Sévres, France. aid the International Committee 
come decisions technical matters there are number consultive committees consisting 
distinguished experts representing the major national physical laboratories and others 
appointed for their individual distinction the field concerned. Each under the chair- 
manship member the International Committee. These advisory committees are 
present six number and cover the following fields: electricity, photometry, thermometry, 
the definition the second (of time), the definition the metre, standards measure for 
ionizing radiation. 

earlier decision taken the 10th Conference 1954 was great importance and 
might also noted. This was resolution which authorized the International System 
Units (abbreviation S.I.), which give the basic dimensions and the corresponding units for 
our scientific measurement system. These are: 


3 


| 


| 
| 


ANNOUNCEMENTS 641 


length, metre; mass, kilogramme; time, second; electric current, ampére; tempera- 

ture, thermodynamic degree; brightness, candela. 
system evolution the original metric system consisting only length, metre; 
mass, kilogramme. addition being the basis all scientific measurement also 
the metric countries, now constituting substantial world majority, the basis all legal 
RECEIVED MARCH 2, 1961. L. E. HowLettT 
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